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Abstract 



Bound states of heavy quarks can serve as a laboratory for inquiry into the behaviour of 
the fundamental strong and electroweak interactions. This thesis examines observations of 
i?^, -B", and mesons produced in proton-antiproton collisions at centre-of-mass energies 
of -s/s = 1.8 TeV. The S-meson decay products are recorded using the Collider Detector 
at Fermilab (CDF) located on the Tevatron collider at the Fermi National Accelerator 
Laboratory in Batavia, Illinois, USA. 

Four S-meson decays and their charge conjugates are studied: J/ipK'^, 
J/^K*(892)°, B+ il){2S)K+, and B^ V(25) ii:*(892)°. Using a data sample corre- 
sponding to jCdt = (109 zb 7) pb~-^, statistically significant signals are observed in all the 
channels. Topological similarities between the B decays are exploited to measure the six 
relative branching fractions {B) of each channel with respect to the other channels. The 
ratios involving the B^ J /ij^ mode are 

B{B^ J/ipK*{892f^ 



B{B+ V>(2g) K+) 
B{B+ J/iIjK+) 

B{B°^ip{2S) i^:* (892)0) 



= 1.76 ± 0.14[stat] ± 0.15[syst] 
= 0.558 ± 0.082[stat] ± 0.056[syst] 
= 0.908 ± 0.194[stat] ± 0.100[syst]. 



B{B+ J/V' K+) 

The indicated uncertainties are statistical and systematic, respectively. In ratios involving 
unlike B-meson species, equal production rates for and B^ mesons have been assumed. 

Absolute branching fractions are extracted by normalizing the above ratio measurements 
to the world-average value, B{B+ J/i^K^) = (1.01 ± 0.14) x lO'^: 

B{B^ J/V'K*(892)°) = (1.78 ± 0.14[stat] ± 0.29[syst]) x 10"^ 

B{B+ 'ip{2S) K+) = (0.56 ± 0.08[stat] ± 0.10[syst]) x lO'^ 

B{B^ ^l)(2S) K*{892f) = (0.92 ± 0.20[stat] ± 0.16[syst]) x 10"=^. 

The B+ 'il^{2S)K+ and 5° iIj{2S) K* {892f reconstructions are the first observa- 
tions of these processes at a hadron collider. The branching-fraction ratio measurements 
are consistent with phenomenological predictions that employ the factorization Ansatz, and 

iii 



the absolute branching-fraction measurements are consistent with previous world-average 
limits and values from e"*" e~ colliders. The measured branching fractions that involve ip{2S) 
final states constitute the world's most precise measurements of these quantities. 
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Chapter 1 

Introduction 



At the high-energy frontier of experimental particle physics, the last quarter of this cen- 
tury has seen the discovery of a new generation of heavy and (apparently) fundamental 
constituents of matter, the bottom and top quarks. This dissertation describes a study of 
some of the properties of the bottom-quark member of this new generation, properties that 
can aid in the understanding of matter and the fundamental forces that act upon it. The 
present chapter opens with a chronology of experimental results that puts into context and 
motivates the study of exclusive nonleptonic decays of B mesons, composite particles that 
consist of a bottom quark and a light antiquark. 

An understanding of the properties of all types of matter, be it particulate or astronom- 
ical in scale, appears to require the existence of agents, or forces, to mediate its interactions. 
Unlike the infinite-range forces of quotidian experience, namely electromagnetism and grav- 
ity, the 'weak' and 'strong' nuclear interactions did not receive recognition until the decades 
near the turn of this century. Fifteen years before Rutherford reported evidence for the 
atomic nucleus in 1911 iQ, Becquerel demonstrated the existence of spontaneous radioac- 
tivity emitted by phosphorescent substances Q, a harbinger of a hitherto unknown 'weak' 
force. The 'strong' force, the first evidence of which was reported by Chadwick and Bieler 
in studies of a-particle collisions with hydrogen nuclei Q, was postulated to bind together 
the known baryons (protons Q and neutrons [||) into atomic nuclei. The concept of the 
'meson' was introduced by Yukawa in 1935 as a hypothetical mediator of the strong nuclear 
force experienced between proximate nucleons [^. Yukawa's meson, called the pion (vr), was 
eventually discovered in cosmic rays in 1947 after having been mistaken for the muon, 
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the first evidence of which was reported a decade earher [^. By the 1950s, the evolution 
of cosmic-ray experimental techniques and the advent of the modern particle accelerator 
made possible the study of a more massive type of matter in the form of 'strange' mesons 
and baryons. Whereas the production of the new strange matter was thought to involve 
the strong interaction, the long lifetime of the strange particles ('^10^^'^ s) indicated that 
the weak interaction played a primary role in their decay Q. The subsequent "explosion" 
of new particle discoveries during the 1960s and 1970s hinted that two large classes of par- 
ticles, the mesons and baryons, were composite, just as the multitude of different atomic 
elements foretold the existence of subatomic structure in the early years of this century. 



1.1 The Fifth Quark 

The notion of quarks found its origins in the early 1960s in the course of searches for an 
organizing principle to describe the proliferation of hadronic particles and resonances ob- 
served by the experiments of that time |jT^. To this end, Gell-Mann [ pi] ] and, independently, 
Ne'eman |12| refined an application of the SU(3) representation|^, which was originally for- 



mulated in terms of "fundamental" p, n, and A baryons [|T3|, to introduce an organizational 
framework of the known baryons and mesons. The SU(3) approximate-symmetry inter- 



pretation was further extended by Gell-Mann |14] and, independently, Zweig |15| | with a 
hypothesis that hadrons consisted of 'quarks'. Gell-Mann and Zweig hypothesized three 
'flavours' of quarks: up (u), down (d), and strange (s). 

Although the quark idea met with immediate success by explaining the taxonomy of the 
experimentally observed particles and resonances, evidence for quarks as dynamical objects 
was to come from future experiments. Studies of the deep inelastic scattering of electrons 
by protons, where the incoming electron scatters off the target proton to produce a massive 
hadronic recoil system, were able to probe the structure of nucleons with measurements of 
the differential scattering cross section as a function of the recoiling hadronic invariant mass 
for different values of four-momentum transfer between the electron and proton. In the re- 
coil invariant mass region beyond the resonances, the cross section exhibited only a weak 



dependence on the momentum transfer between the electron and the hadronic state |16|. 



^ SU(3) denotes the 'special unitary' group in three dimensions, where the matrix operator that effects 
transitions between members of this group is unitary and has determinant +1 (special). 
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This observation, which was cahed 'scale invariance' because the cross section appeared to 
be independent of the momentum transfer and to depend on a dimensionless quantity that 
related the momentum of the recoiling hadronic system to the incident proton momentum, 
gave credence to quark parton models of nucleon structure advanced by Bjorken |17] and 



Feynman |18], which predicted this scaling behaviour. Further confirmation of this par- 
ton picture was supplied six years later by observations of azimuthal asymmetries in the 
production of hadrons using electron-positron annihilation . 

In spite of the parton model's successful description of deep inelastic scattering results, 
the quark idea suffered from a theoretical deficiency when used to classify, for example, the 
baryonQ resonance [pO|]. The baryon, thought in the quark model to be made up 
of three u quarks, each in the same state, appeared to be symmetric under the interchange 



of any of the constituent quarks, thereby violating the Pauli principle |21, E3]. An exact 



SU(3) symmetry, that of a new quantum degree of freedom for quarks, termed 'colour' 



was hypothesized |23| to resolve the conflict with Fermi statistics by rendering the A"*""*" 
wave function antisymmetric. With this new colour symmetry, there also emerged a theory 
of the strong force called quantum chromo dynamics, or QCD, which described the strong 



interactions as being between spin-1/2 quarks and mediated by spin-1 gluons |24, 25, 26|. 
It was postulated that all flavours of quarks and antiquarks were each endowed with one 
of three colours and anticolours, respectively. Unlike photons, which do not carry electric 
charge, the gluon mediators carried colour and were thought to form a colour SU(3) octet. 

QCD also postulated that all naturally occurring particles are colour SU(3) singlets. 
This was motivated largely by the inability of the experiments to produce isolated quarks. 
In QCD, since gluons are themselves colour sources, they are self-interacting; this property 
incites the QCD coupling to grow in strength as the separation between two colour sources 
increases. If two quarks within a hadron or a pair of hadrons are made to recede from each 
other in an energetic process {e.g., in a particle accelerator), the potential energy accrued by 
the increased interquark separation will make it energetically favourable for pairs of quarks 
to be produced from the vacuum and to interact with the receding particles and with each 
other. This process continues until all the quarks are once again confined inside hadrons 



^In this dissertation, references to specific diarge states imply the additional charge-conjugate state, 
unless obviated by context or noted otherwise. 
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producing 'jets' of particles. Evidence of such jets arising from energetic quarks was first 
reported in 1975 in e+ annihilation studies [27|. As predicted by QCD, jets due to gluon 



Bremsstrahlung in processes like e 'J* ^ qqg were subsequently discovered at higher 



centre-of-mass energies [28 



An explanation of the revolutionary and largely unexpected 1974 discovery of the J/ip 
meson in e~ annihilation and p-Be fixed-target experiments [^9| proved to be one of the 
quark parton model's greatest achievements. The unusually high mass and long lifetime 
of the J/V' meson indicated the presence of fundamentally new physics; in turn, the quark 
parton model established the observation as a manifestation of a fourth quark, charm (c), 
in a bound state with its antiquark to form the J/V' meson. 

The interpretation of the J/ ^p meson as a cc bound state was buttressed by the discovery 
of the i/y , or i/)(2S'), meson in its e"*" e~ decay channel |^l|, a resonance that was identified 
immediately as a radial excitation of the J/ip state. Quite analogously to electrodynamic 
explanations of the states of positronium observed two decades previous, QCD was able 
to predict the charmonium cc bound states and their narrow widths. The subsequent 



experimental observation of the decay ip{2S) J/ipir'^ 7r~ [32] served to complement the 
dilepton channels in clarifying the spectroscopy of the charmonium system. 

The discovery of charm brought the count of known fundamental fermions to four quarks 
and four leptons (the electron, e~ , and muon, fi~ , and their associated neutrinos) in two 



generations, thereby vindicating theoretical prejudice toward a lepton-quark symmetry |33|. 



This symmetry was temporarily broken by the discovery of the r lepton in e"*" e colli- 



sions [34 1, setting the stage for a third generation of fundamental fermions. 



1.1.1 Discovery 

Evidence for the b quark, often referred to as the 'bottom' or 'beauty' quark, was initially 
obtained with techniques analogous to those used in the discovery of charmonium. In 1977, 
a significant excess in the rate of dimuon production was observed in 400-GeV proton- 
nucleus collisions by a fixed-target experiment at Fermilab [^5|. The original enhancement, 
observed near 9.5 GeV/c^, was interpreted as arising due to decays of bottomonium, a bb 
bound state (T), and was rapidly confirmed and resolved into two resonances, namely the 
T{1S) and T(25) mesons M. 
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Extensions of nonrelativistic potential models tuned on the cc system were used to 
calculate the dielectronic partial width of the T state, r(T — > e"*" e~), which depended on the 
electric charge carried by the b and b quarks. Comparisons of the calculated r(T — > e"*" e~) 
with the area under the observed T line shape suggested that the b quark would join its d 
and s-quark counterparts in possessing an electric charge of —1/3. The emerging pattern 
of generations implied by the properties of the observed quarks^ helped presage the 1995 



discovery of a sixth quark, top (t), observed in 1.8-TeV pp collisions |^, |40|| . 
1.1.2 Electroweak Interactions 

Whereas the weak interactions of leptons have been observed to be strictly intragenera- 
tional [41 1, quarks in their mass-eigenstate generations. 



XI) 



may interact weakly with quarks in generations other than their own. The typical notation 
used to describe the degree of this 'mixing' is the Cabibbo-Kobayashi-Maskawa (CKM) 



matrix [||, ||, 




Vud Vus Vub 

Vcd Vcs Vcb \ s \ , (1.2) 




Vtd Vts Vtb 

which, by convention, leaves the +2/3-charged quarks unmixed; the states (d', s', b') are the 
weak eigenstates. Under the constraints that there be three quark generations and that 
the CKM matrix be unitary, the mixing can be parameterized with three angles and one 
complex phase. The pursuit of measurements to constrain the CKM matrix and observe 
the charge-parityQ (CP) violation in the 6-quark system (a result of a non-zero complex 
phase) constitutes a major component of many experimental programmes at modern particle 
accelerators [^] . 

Contemporary understanding of the weak interactions has its foundations in Fermi's 
field theory of (5 decay, introduced in 1934 EH]. The idea of the four-fermion interaction 



'^Analogous to the need for the charm quark in the second quark generation to cancel out unobserved 
strangeness-changing neutral currents in the electroweak theory (GIM mechanism ^^), the lack of observed 
flavour-changing neutral currents {e.g., B'' fi^ fi~) in the 6-quark system inferred the existence of a weak 
isospin partner to the third-generation b quark. Moreover, the weak isospin of the b quark was measured, 
via angular asymmetries in e"*" e~ bb production, to be Ts = —1/2 [B7|, suggesting a doublet structure 
similar to that for the lighter generations. 

^'Parity' refers to the quantum mechanical operator that inverts a spatial displacement vector, i.e., 
r — > — r. Parity violation was first observed in studies of /3 decay in polarized Co*'" atoms |44[. 
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was retained for several years before Sakurai introduced the universal V — A (vector and 
axial-vector current) modification [^] to accommodate the experimentally observed parity 
violation. Unfortunately, the V — A Fermi theory had the shortcomings of unitarity viola- 
tion (the nonconservation of probability in predicted cross sections that grew quadratically 
as a function of centre-of-mass momentum) and nonrenormalizabilityQ in its predictions of 
cross sections at high energies. The 'standard model' of electroweak interactions, developed 
primarily by Glashow, Weinberg, and Salam in the 1960s ||4^ and based on the gauge group 
SU{2) X U{1), hypothesized four intermediate gauge fields to avoid these difficulties: the 
W^, , Z^, and 7 bosons. The first three of these were thought to be endowed with 



mass via the Higgs mechanism [50|; the fourth (7) boson was the massless photon of the 
electromagnetic interaction. The weak and electromagnetic forces were thence consolidated 
into a single theory. Soon after the Glashow- Weinberg-Salam standard model had been 
proven to be renormalizable |l5l[| , neutral weak-current interactions, predicted by the exis- 
tence of the boson in the theory, were discovered in a neutrino-antineutrino experiment 
using the Gargamelle bubble chamber at CERN in 1973 [^]. Complete confirmation of 
the electroweak theory, however, came a decade later with the discovery of the and Z^ 
intermediate vector bosons p^]. 

The term in the standard-model electroweak Lagrangian that plays a large role in the 
6-quark decays investigated in this dissertation represents the charged-current weak inter- 



action between the fermion fields |3C]: 

^CC = E?^ - + T-WP Ci, (1.3) 

where g is the SU{2) gauge coupling constant, 7'^(1 — 7^) are the Dirac 7 matrices repre- 
senting the V — A current, T~^ and T" are the weak-isospin raising and lowering operators, 
respectively, are the massive weak charged boson fields, and the index i represents the 
fermion families. In charged-current weak interactions of the b quark, the fermion fields, 
^i, are either left-handed SU{2) doublets, ^3 = ( ^, ) , or right-handed SU{2) singlets, 



L 



(^3 = (&')^, where b' is the weak eigenstate defined in Equation 1.2 



theory is renormalizable if the predicted amplitudes of physical processes remain finite at all energies 
and for all powers of the coupling constant, often at the expense of the introduction of a finite number of 
arbitrary experimentally-determined parameters ]4m]. 
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1.1.3 Hadroproduction 

In the present study, b and b quarks were produced m 1.8-TeV cohisions of protons (p) and 
antiprotons (p) by way of the inclusive process 

p{kp) + P{kp) ^ b{kb) + b{ki) + X, (1.4) 

where X denotes the 'underlying event' and kp and kp (kh and /cj) are the momenta of the 
baryons (6 quarks). Note that the p [p) baryons each comprise several partons: the uud 
(uud) 'valence' quarks, gluons, and many 'sea' quark-antiquark pairsQ. 

A perturbative QCD formula for the invariant differential hadroproduction cross section 
of a & quark with energy Ef, and mass rrii, can be expressed by convolving the partonic 
cross section (a) with the parton distribution functions of the two hadron reactants in the 



form [54 1 
Ebd^a 



V / dxi [ dx2 
77 "'0 Jo 



d^kb 



d^kh 



F[{xuQ^)Ff{x2,Q 



^^^^^^-^{xikp,X2kp,ki,;mb,iJ.,A) i^^^-c^pi^^ /^2^ 

(1.5) 



where xikp and X2kp are the momenta of the incoming partons, Ff'J' are the parton distri- 
bution functions for the i^^ and j^^ parton in the p and p baryons, respectively, and is 
the square of the four-momentum transfer. The parameter ^ represents the energy scale of 
the process and, by assumption, \Q\ = /x. The quantity A is an experimentally-determined 
parameter used in the description of the dependence of the strong coupling constant, a^, 
on the energy scale, /i. 



Integrating Equation 1.5 over the momentum k^ yields the total cross section for the 
production of a 6 quark, 

^(^) = 5Z / '^^^ i dx2aij{xiX2s;mb,n,A)Ff{xi,f/)Ff{x2,n'^), (1.6) 
^Jo Jo 

where s is the square of the centre-of-mass energy of the colliding proton and antiproton. 
The threshold condition for bb production is met when the square of the parton-parton 
centre-of-mass energy, s = xiX2S, satisfies the condition s = 4m^. 

The heavy mass of the b quark makes possible QCD calculations of perturbation 
series in powers of the running strong coupling constant, a^. The first terms in the series 
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The 'underlying event' refers to the aggregate product of lower-energy interactions between those partons 



not directly involved in the "hard" scattering part of Equation 



1.4. 
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that contribute to the cross section are 0{a'^) quark-antiquark annihilation or gluon-gluon 
fusion processes: 

q + q b + b (1.7) 

g + g ^ b + b. 



Figure LI depicts Feynman diagrams of these lowest-order bb production mechanisms |55|. 
The next-to-leading 0{al) terms in the perturbative series arise from processes ||5^ like 

q + q b + b + g 

g + g ^ b + b + g (1.8) 
g + q b + b + q 
g + q b + b + q, 

some examples of which are illustrated with Feynman diagrams in Figure Due to 

interference with diagrams containing virtual gluons, the two processes in Equation |1.7| 
can also contribute at 0(0^). For high energies, i.e., when ki{b) S> tti;,, where A;t(6) is 
the momentum of the b quark projected onto a plane perpendicular to the axis of the two 
incoming partons, some of the next-to-leading order 0{a^) mechanisms can contribute to 
the cross section by amounts comparable to the ©(a^) contributions [54|. 



In Equations and 1.6, there is a degree of arbitrariness in the value of the renor- 
malization scale, /i, that contributes a relatively large uncertainty to QCD predictions of 
6-quark production because they are not calculated to all orders in Og. The value of fi is 



typically assigned to be near a physical scale, such as irib or ym'^ + kq^{b); however, these 
choices of /i are "bootstrapped" because the fact that b quarks are confined inside hadrons 
requires that extractions of mi, depend on the renormalization scheme used, model-specific 



definitions of rrn,, and the value of n itself [30, 5€]. 



(Pa 

Nason, Dawson, and Ellis (NDE) have calculated — „ , the inclusive differential 

dyb dkq^[b) 

6-quark production cross section, as a function of rapidity, yi,, and A;t(&) fsj]. Rapidity is 



a measure of the polar angle of a particle's trajectory, usually with respect to the collision 
axis, and is defined for a b quark as 



Eb + k,{b) 
Eb-k,{b) 
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Figure 1.2: Examples of some of the Feynman diagrams for the next-to-leading order 
{0{ag)) mechanisms of bb hadroproduction. 
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Collision Type 


PP 


PP 


e+ e 


Accelerator 


Tevatron 


LHC 


CESR 


LEP, SLC 


[GeV] 


1800 


14 000 


10.58 


92 




20-40 


~200 


~ 1.1 X 10-3 


~ 9 X 10^3 




4 X lO^'' 


2 X 10-^ 


0.25-0.33 


~ 0.2 



Table 1.1: Comparison of hadroproduction and leptoproduction cross sections for bb pro- 
duction at various centre-of-mass energies |57|. The parameters of the five accelerators are 



summarized in Reference |30| 



where kz{b) is the projection of the 5-quark momentum onto the beam axis. The rapidity 
variable is useful to descriptions of high-energy particle production because the shape of 
the particle- multiplicity distribution, dN/dyb, is Lorentz-invariant; reference frame trans- 
formations amount to linear shifts in the origin of [50] . 



Table 1.1 compares the maximum 6-quark production cross section {o'^^^) applicable 
to the present study (Tevatron) with that in other collision types at different energies. 
Although the rates of bb hadroproduction exceed those of leptoproduction by three orders 
of magnitude, they constitute a significantly smaller fraction (a^g/cTtot) of the total cross 
section (cxtot) than do the analogous bb leptoproduction fractions. The consequentially low 
signal-to-background ratios for 6-quark production at hadron colliders pose experimental 
challenges that must be overcome. 



1.2 The B Meson 

Analogously to studies of the charm sector, the discovery of the bottomonium states lead 
naturally to the expectation that, at masses somewhat greater than those of the lowest-lying 
T resonances, "open" bottom meson production would occur in the form of bare B'^ (bu) 
and (bd) mesons. First evidence for i?-meson production, reported in 1981 by the CLEO 
collaboration ||58|] , was obtained through observations of increases in the single electron and 
muon inclusive cross sections in e~ collisions. The enhancements were attributed to 
inclusive semileptonic decays of B mesons via the processes B X ii'i, where i denotes 
either the e or fj, lepton flavour and X represents the remaining hadronic system. 
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Prior to the first evidence for the existence of B mesons, however, Pritzsch argued that 
". . . the only reahstic method to discover the B mesons" was through the reconstruction of 
their decays to charmonium states [59|. Fritzsch's assertion, which launched a considerable 
amount of theoretical work on the subject (see, for example, References [Q), was moti- 
vated by the following points: cc mesons are readily produced in weak i?-meson decays; 
the /i^ and e+ e~ decay modes of the cc states can be identified easily, most notably in 
hadronic collision environments with their sizeable backgrounds; and the relative heaviness 
of the charmonium states forces the remaining hadronic system in each B decay to be rel- 



atively simple because of the attendant restricted multiparticle phase space [59|. In spite 
of theoretical expectations that B mesons would first be reconstructed in their charmo- 
nium final states, it is interesting to note that the first full reconstruction of B decays was 
achieved using final states containing and Z?*(2010)^ mesons, such as B^ vr^ and 

B^ D*(2010)^ vr^ vr^ |^l|. More mention of the mechanism of charmonium production 



in 5-meson decays, which is the subject of this investigation, is given in Section 1.2.2 and 
thereafter. 



1.2.1 Hadronization 

The 'fragmentation' of a 6 quark into a colour-singlet hadron, in this case a B meson, is 

a long-distance, nonperturbative QCD process. Models of fragmentation typically employ 

Eb+P\\{B) _ 

a parameter Zf, where Zf = ^ ^ ^^^^ |62| ], which represents the fraction of available 
energy-momentum carried by the B meson. The symbols Eb and p^^{B), respectively, 
represent the hadronized i?-meson energy and the momentum component parallel to the 
direction of fragmentation. Simple kinematical arguments support the claim that functions 
describing b ^ B q processes, where q is the light-quark counterpart to the q quark in 
the B meson, peak at high values of due to the expectation that the majority of the 
6-quark momentum is transferred to the B meson |63, 64]. Peterson et al. noted that the 



quantum mechanical amplitude describing fragmentation would be inversely proportional 



to the energy transfer of the process, AE = Eb + Eg — Ef, |65]. Assuming nib ~ ^B, then 



^ml + z}k{b)^ + + (1 - zj)2A:(6)2 - + A;(6)2. (1.10) 
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The fragmentation function, Dff{zf), was estimated from Equation I.IC with the expression 

1 1 



1 1 ^ 

2/ 1-^/ 



where the factor of in the denominator arose from longitudinal phase space and the 
quantity ef, ~ nig /ml was the "Peterson parameter" p5| . 

The Peterson parameterization benefits from its dependence on a single experimentally 
determined parameter (e^). Chrin has estimated a value of ef, = 0.006 it 0.002, which is 
based on a survey of several experimental e~ observations [|6^ ] . The search continues for 
a better understanding of 6-quark fragmentation and its sensitivity to the type of collision 
environment in which the b quark was produced and the flavour(s) of the non-6 quark(s) 
constituting the final-state hadrons. 

1.2.2 Nonleptonic i?-Meson Decays 

Nonleptonic decays of B mesons feature a rich interplay of the weak and strong interactions. 
Knowledge of i?-meson decay rates may be used to obtain angle and phase information in the 



CKM matrix (Equation 1.2) by way of the charged-current electroweak interaction described 
in Equation |1.3| . The practical extraction of information about the weak processes, however, 
is confounded by the fact that quarks are necessarily confined inside bound colour-singlet 
states, requiring the invocation of quantum chromodynamics to complete the description. 
The relatively heavy mass of the b quark, nominally 4.1 to 4.5 GeV/c^ or about 5 times 
the mass of a proton, is such that QCD descriptions can be decoupled into 'short' and 
'long-distance' dynamical regimes |Q (see Section 1.2.3| ). Short-distance effects exploit 



the asymptotic-freedom property of QCD, which enables the perturbative calculation of 
corrections to the electroweak Hamiltonian due to exchanges of hard gluons []67|| . In the 
context of nonleptonic decays of B mesons, long-distance QCD effects generally involve the 
hadronization of the decay products and include the exchange of soft gluons, the creation 
of qq pairs from the vacuum, and interactions of the hadrons in the final state. That long- 
distance QCD processes are nonperturbative has so far thwarted calculations of these effects 
from first principles; nevertheless, the separation of the decay mechanism into short and 
long-distance components is currently the most successful approach to understanding many 
S-meson decay processes. Experimental measurements of S-meson decay properties are 
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C 
C 



u,d * * u,d 



Figure 1.3: Feynman diagram of the 'spectator-internal' weak decay mechanism for a B 
meson (here either a bu {B~) or bd {B'^) state) decaying to charmonium {J/ip or ip^lS)) 
and a strange meson {K~ or K (892)'^). The decay is colour-suppressed because it only 
occurs when the cs pair, itself a colour singlet, conspires with the c quark and the u or d 
quark to form colour-singlet cc and su or sd mesons. In this mechanism, the u or d quark 
is assumed to be a 'spectator' of the weak process. 

clearly essential to the pursuit of an understanding of the dynamical complexities of the 
underlying phenomena. 

This dissertation concerns the experimental study of exclusive nonleptonic decays of B 
mesons into J/ip or tp{2S) vector- meson and K'^ pseudoscalar-meson or i^*(892)'^ vector- 
meson[] final states. Figure exhibits the weak-interaction aspect of these decays, which 
are categorized as 'spectator-internal' processes by virtue of two assumptions: (i) the light- 
quark (spectator) constituent of the B meson fails to participate in any process before 
hadronization and final-state interactions occur and (ii) the quark with which the spectator 
forms a colour singlet in the final state derives from an "internally produced" W boson. 
These decays are considered to be colour-suppressed because they can only occur when 
the W bosons' hadronic decay products, themselves constituting colour singlets, conspire 
with the charm quarks from the flavour-changing decays and the spectator quarks to form 
colour-singlet charmonium and strange mesons, respectively. Finally, these processes are 
CKM-favoured, as their rates depend upon \Vcb\'^ (refer to Equation |l.2|) , where \Vcb\ = 



0.041 ±0.003 |0| 



Gourdin, Kamal, and Pham state that "it is generally believed that the best place 
to study the importance of colour-suppressed processes in 5-meson decay is to look at 



final states involving a charmonium and a strange meson" |70|. Precise measurements 



^The K*{892)^ vector meson was the first meson resonance to be observed [pSl p9 
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of properties of this class of decays can improve knowledge of wave functions in the B 
system [^T[, the same wave functions that are involved in, e.g., the rare processes 
K*{892)^ J and B K [71, 72 1. Experimental observations oi B ip K decays 

also provide input to and tests of theoretical and phenomenological methodologies and 



hypotheses. Some of these issues are discussed in Section 1.2.3 
1.2.3 The Factorization Ansatz 

The factorization Ansatz model was first used by Bauer, Stech, and Wirbel (BSW) to de- 
scribe exclusive semileptonic decays of heavy mesons The techniques were subsequently 
extended to encompass nonleptonic heavy- flavoured meson decays [^^. A fundamental as- 
sumption in the BSW model is the distinction, noted in Section 1.2.2| , between two disparate 
time (or distance) scales in the decay dynamics. In the nonleptonic weak decay of a i? me- 
son, the h quark decays with a time scale of r ~ \/Mw ~ 10^^^ s. It is assumed that all 
other partons in the system (including the spectator quark, sea quarks, and soft gluons) 
suspend any interaction until the longer QCD time scale, r ~ 1/Aqcd ~ 10^^'^ s, indicating 
the time when confinement forces become important, takes effect. 

Using the mode B^ J/ipK~^ as an example, the amplitude of the decay {^A) can 
be expressed as a superposition of local operators, Oj, with scale-dependent short-distance 
Wilson coefficients, Ci{fj,), in the form 



{J/i^,K+\n,s\B+) 

-^^Q(^)(JMi^+|o.|s+), 



(1.12) 



where Gp '^s the Fermi coupling constant and TicS is the effective Hamiltonian embodying 



the 0{as) hard-gluon corrections to the weak decay of the b quark [75|: 



Weff = -^^cb {Clifl) [{Cbr (s'cU + C2ifl) [{s'br (ccu } 



(1.13) 



The primed fields represent weak eigenstates, which were defined in Equation 1.2. Similar 



to the form of Equation 1.3, Equation 1.12 uses the notation 



(929l)"^E«"2A7"(l 
A 



(1.14) 



for the current terms, where A is the colour index. Currents proportional to the Wilson 



coefficient Ci(/i) are charged, as exemplified in Figure 1.3, and currents proportional to 
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C2(/u) are effectively neutral. In principle, gluonic penguin^ interactions can also contribute 
to Equation 1.13| ; however, these are highly suppressed due to the fact that the production 
of a J/^p meson from gluons necessarily involves no less than three gluons [76|. 

The factorization Ansatz is an approximation whereby one of the currents in the (cur- 
rent) X (current) form expressed in Equation 1.13| is assumed to be asymptotic, thereby 
enabling the factorization of the amplitude in terms of hadronic (H) currents instead of 



quark currents ||7J, ^ : 

MiB+ ^ J/i;K+) = -9Za,V,,[{J/^p\icc)^\0) • (i^+|(s-'6)^|S+)} , (1.15) 
where the quantity 02 is a new scale-independent coefficient native to spectator-internal 



mechanisms (see Figure 1.3) and constructed via a linear combination of the Wilson coeffi- 
cients evaluated at a scale defined by /x ~ m^. 



02 ~ C2(mb) +^Ci{mb). 



(1.16) 



The new parameter ^ is an ad hoc colour factor, often taken naively to be ~l/3 due to 
the three degrees of colour freedom^. Since the colour structure in Equation |1.1(: can 
easily be destroyed by long-distance soft-gluon dynamics, 02 is typically treated as a free 
parameter 173, 17^ . Equation 1.15 suggests that a knowledge of the B ^ K decay rates 



can furnish estimates of the magnitude of the 02 parameter |75, 79 1. 

The asymptotic (J/'0|(cc)^|O) single-particle matrix element in Equation 1.15 is pro- 
portional to the J/V' decay constant, thus reducing the amplitude calculation to a deter- 
mination of the hadronic form factors constituting the {K~^\{s'b)fj\B~^) matrix element. 
Several different approaches to modeling the hadronic form factors, which typically employ 
measurements from experiments involving semileptonic decays, exist in the literature and 
will not be discussed here in detail. One elegant and successful method has been to exploit 
the heavy-quark symmetries that arise in the rub 00 limit |8C]. The application of heavy- 
quark symmetries to the determination of B ^ K form factors, however, is hampered by 



^Penguin decays have the characteristic that the W boson is reabsorbed by the quark hne from which 
it was emitted, yielding an eflective net flavour-changing neutral current process. A photon, gluon, or 
boson is emitted from the resulting loop. 

^As the choice of the notation a2 suggests, there is another scale-independent coefficient, ai, defined by 
ai ~ C'i{mb) + £,C2(rnb), which represents a separate 'spectator-external' class of _B-meson decays. 
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the relative lightness of the s quark. Moreover, the observed absence of tree-level flavour- 
changing neutral-current decays complicates the use of experimental inputs in estimates of 
B ^ K form factors. Highly model-dependent Ansdtze have been used to estimate B ^ K 
form factors from, e.g., D K^*^ Iv semileptonic decay measurements [81, p^]. Two prin- 



cipal sources of significant theoretical uncertainty are present in all these estimations: the 
uncertainties in the numerical values of the form factors at zero momentum transfer be- 
tween the B parent and the K daughter {q^ = 0) and the uncertainties in the assumed 
parameterizations of their cp' extrapolations fs^ . 

It is important to point out that, unlike for semileptonic decays where the amplitude 
can be decomposed into leptonic and hadronic currents, there is no theoretical justifica- 
tion for the application of the factorization Ansatz to colour-suppressed B ^ tp K decay 



modes |84, Tests of the validity of factorization for B ^ K decays are challenged 
by the difficulty in isolating factorization effects from form-factor effects. Part of the am- 
biguity can be removed by confronting the data with calculated observables that involve 
ratios of form factors, thereby purging some reliance on the absolute values of the form 
factors, but nevertheless retaining a dependence on assumptions about their q'^ extrapo- 
lations |8^]. Tests of factorization have consisted of requiring that the models reconcile 



ratio-of-branching-fraction measurements {e.g., B{B J /ip K)/B{B — > ip{2S)K)) with 
measurements of longitudinal polarization fractions {e.g., —{B^ J/ip K*{892)^)). Fac- 
torization models have been shown to be unable to account simultaneously for measurements 



of these two observables [70, 82 [. While inadequacies in the form- factor approximations are 



clearly possible, it has been suggested that the discrepancies may be due to nonfactorizable 



contributions to the decay amplitudes [g^j 87, 88, 8£, 90[. 
1.2.4 Theoretical Predictions 

The foregoing synopsis of factorization techniques points out some of the difficulties and 
uncertainties in constructing a reliable theoretical and phenomenological picture of exclusive 
S-meson decays to charmonium and strange-meson final states. Although the present 
tendency is for the experimental data to drive phenomenological investigations of form- 
factor assumptions and nonfactorizable contributions to the oi and a2 parameters, there 
nevertheless exist some recent branching-fraction predictions in the literature. 
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Channel 


DL dncniiig-r i dc l 
Deshpande 71] 


lUn JT I cCLlCLlUilo 

Deandrea Q 


P) l^^'J I 
Chene isl 


B+ J/ip K+ 

J/V'K*(892)° 
B+ 7p{2S) K+ 
B^ i){2S)K*{892f 


0.84 
1.63 
0.33 
1.27 


1.1 ±0.6 
1.6 ±0.5 
0.37 ±0.19 
0.74 ±0.23 


n/a 
n/a 
0.52 
0.76 


Year Published 


1990 


1993 


1997 



Table 1.2: Theoretical branching-fraction predictions, based on the factorization Ansatz, of 
the decay modes under study. The Deandrea et al. predictions assumed i?-meson lifetimes 
of tq+ = T^o = 1.4 X 10^^^ s and the Cheng et al. predictions used the world-average 
measured lifetimes: = 1.62 x 10^^-^ s and r^o = 1.56 x 10~^^ s. 



The three sets of factorization predictions considered in Table |l.2| were selected because 



they all included treatments of the tp{2S) final states. The Deshpande et al. results [71 
made use of form factors from BSW, who calculated solutions to a relativistic harmonic 
oscillator potential model and assumed a single-polef^g^ dependence jT^, 74 1. The Deandrea 
et al. calculations employed form factors that were derived from experimental exclusive 
semileptonic D-meson decay measurements and that were also taken to possess a monopole 
dependence [31|. The Cheng et al. predictions]^ for the 1^(23) final states were based both 
on an 02 parameter that was computed from experimental measurements of the analogous 
J/^/j modes and on form factors that were calculated explicitly over the entire physical 
range [p6[. 



1.3 Dissertation Overview 



This dissertation describes searches for the decay channels listed in Table 1.2 and relates 



details of their full reconstruction and branching-fraction (B) measurements, as observed 



in 1.8-TeV proton-antiproton collisions. Figures 1.4 and 1.5 schematically illustrate the 
topologies considered in the analysis of B^ and B^ mesons, respectively. Both species of 
charmonium meson are sought in their dimuon modes, and the ip{2S) meson is supplemen- 



^''The gen eric monopole form-factor formula is F(q^) = F{0)/{1 — /m'^), where m represents the pole 
mass |7|, 111. 

^^Note that the Cheng et al. predictions were compared in Reference with preliminary versions 
of the measurements described in this dissertation. 
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B- 



^(25) K+ 



Figure 1.4: Schematic diagrams of the decay modes reconstructed in this study. 



B^ 



J/il) (892)0 



'>p{2S) K*{892)^ 



-A'+ TT" 



Figure 1.5: Schematic diagrams of the B^ decay modes reconstructed in this study. 



tally reconstructed in its hadronic-cascade ip{2S) J/ipir^ vr^ mode (refer to Section 1.1). 

Just as the examination of relative, as opposed to absolute, branching fractions can 
simplify theoretical approaches to S-meson decay ( ||8^ and Section 1.2.3| ) , the use of candi- 
date event yields to compute ratios of branching fractions can have several benefits on the 
experimental side. Consider an expression describing the requisite ingredients to measure 
the experimental absolute branching fraction of, for example, the decay B^ — > ip{2S) , 
where the ip{2S) meson decays in its ip{2S) J/'ipTT~^Tr~ channel: 



B{B^ ip{2S)K^ 



(1.17) 



2eAfu ■(y{pp^ b{b) X) ■ JCdt- ^(-0(25) ^ J/ip tt+tt')- B{J/ij /x+/i-) ' 
where N{B^ —>■ il){2S) K^) denotes the yield of candidate decays reconstructed, e repre- 
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sents the product of several reconstruction efficiencies (Chapter |5|), A is the geometric and 
kinematic acceptance correction factor (Chapter is the fragmentation fraction, or the 

probabihty that a b quark will hadronize with a u quark into a meson, a (pp ^ b(b) X) 
is the 6-quark hadroproduction cross section {X denotes the underlying event), jCdt sig- 
nifies the time-integrated luminosity of the data sample, and B{'tp{2S) — > J /ij^ ti^ tt") and 
B{J/ip ;U+/i^) are charmonium branching fractions. The factor of two in the denomi- 



nator of Equation |1.17] accounts for the fact that both and B meson candidates are 
reconstructed. 



As the schematic decay chains in Figures 1.4 and 1.5 connote, several topological sim- 



ilarities exist amongst the decays under scrutiny. The construction of ratios of branching 



fractions like the one expressed in Equation 1.17 can therefore exploit these congruities 
by accommodating the cancelation of several common reconstruction efficiencies and sys- 
tematic uncertainties. Moreover, other quantities with sizeable measurement uncertainties 
concomitantly divide to unity in these ratios. Since the data sample (Section ^?^) is com- 
mon to modes in both the numerator and denominator of a ratio of branching fractions, 
the JCdt factors are no longer important. Depending on which modes make up a given 
ratio, some of the charmonium and i^*(892)'^ branching-fraction factors can divide out of 
the expression, thereby reducing the overall systematic uncertainty. The cancelation of the 
a (pp b{b) X) factors is expected to be especially beneficial, as theoretical models (Sec- 
tion |4.1| ) and experimental measurements of the 6-quark differential hadroproduction cross 
section versus pT, both of which rely on some assumptions, have been observed to have 
consistent shapes but only marginally consistent normalizations p2| . 

An example of a branching-fraction ratio to be measured in this study is given by 



extending Equation 1.17 to the expression 



13{B+ ^^P{2S) K+) _ N{B+ ^{28) K+) e'j^A' fd B{K*(892f K+tt") 



B{BO ^ J/^K* (892)0) N {BO ^ J/^K* (892)0) B(ij(2S) ^ J/i; 7r+ tt^)' 

(1.18) 

where the primed quantities refer to the B^ J/tp K*(892)0 reconstruction, the subscripts 
R indicate that the efficiency products have been reduced from their absolute values due 
to the cancelation of common factors in the ratio, and fd is the fragmentation fraction for 
the b B^ hadronization process. Direct measurements of branching-fraction ratios such 
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R) ^ B{i)/B{j) 


B+ J/i; K+ B^ J/i; K*{892)^ B+ ip{2S) K+ 


J/'ipK*{892f 
tjj{2S)K*{892f 


jjJ/^K* (892)0 
^ J/?/) K+ 

pV(25) K+ „^(2S) K+ 
^J/ipK+ ^J/^K* {892)0 

r,tP{2S) K* (892)0 pV(2S') K* (892)" pV'(2S) K* (892)0 
^J/ipK+ ^J/ipK* (892)0 ^ip(2S)K+ 



Table 1.3: The branching-fraction ratios measured for the various i?-meson final states. The 
ratio Rl is located in the i^^ row and the j^^ column, and the i and ?' indices refer to the 



numerators and denominators of the ratios, respectively. Note that the ratios containing 
ip{2S) mesons are composed of contributions from two separate i?-meson reconstructions. 



as that in Equation |1.18| may be made with the assumption that fu = fd- Measurements 
of fu/ fd that assumed isospin symmetry have confirmed this hypothesis in pp collisions, 
up to an uncertainty of 21% [^]. Reference |Q also assumes /„ = fd on the grounds 
that the B'^ and B^ meson masses are nearly equal and that the CLEO collaboration has 



measured fu/ fd = 1.13ib0.20 |94|. Finally, the world-average [30| B^ J/'ijj K'^ branching 
fraction will be used to extract measurements of those absolute branching fractions listed 



in Table 1.2 from the appropriate measured ratios of branching fractions. 

Table 1.3 lists the branching- fraction ratios investigated. Numerators and denomina- 
tors containing ip{2S) mesons are composed of contributions from two separate S-meson 
reconstructions, namely those involving the ip{2S) /i"*" /i~ and ip{2S) J/ij^Ti^ vr^ decay 
modes. 

Chapter ^ describes the acceleration and detection apparatus used to produce B mesons 
and record their decays, respectively. Techniques invoked to reconstruct i?-meson candi- 
dates and reject background processes are discussed in Chapter |^ whereas Monte Carlo 
methods employed to determine the geometric and kinematic acceptance corrections re- 
ceive treatment in Chapter ^. Chapter ^ traces the reckoning of several efficiencies and 
their associated systematic uncertainties, both of which are used in Chapter ^ to calculate 
ratios of branching fractions from the observed yields of the decays under study. The impli- 
cations of the measurements are discussed further in Chapter |6| and conclusions are offered 
in Chapter |^ 



Chapter 2 

Experimental Apparatus 



The apparatus employed in this study resides at the Fermi National Accelerator Laboratory 
(FNAL) in Batavia, Illinois, USA. The laboratory, which is commonly known as Fermilab, 
is owned by the United States Department of Energy and is operated under a contract with 



the Universities Research Association, Incorporated. Section |2.1| in this chapter briefly 
describes the sequence of accelerators, culminating in the Tevatron synchrotron, that ulti- 
mately accelerate and collide beams of protons against those of antiprotons at centre-of-mass 
energies of 1.8 TeV. 

In this study, observation of the ensuing collision products was achieved through the 
use of one of two general-purpose particle detectors situated at different interaction re- 
gions on the Tevatron collider ring. The Collider Detector at Fermilab (CDF detector) is 
an azimuthally and forward-backward symmetric device that consists of several tracking, 
calorimeter, and muon subsystems. Section provides an overview of the CDF detector 



and Sections 2.3 through 2.6 describe those subsystems apposite to the present analysis: 
the tracking, muon, trigger, and data acquisition systems. 



2.1 The Fermilab Tevatron pp Collider 

The acceleration of protons and antiprotons to energies of 900 GeV is accomplished at 
Fermilab by a synergism of six particle accelerators. The Cockroft- Walton pulsed ion 
source begins the sequence by converting gaseous H2 molecules to H~ ions, which are 
subsequently subjected to a 750-keV electric potential. The H~ ions then enter a 150-m 
linear accelerator, or Linac, where they are accelerated to energies of 400 MeV by a sequence 
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of drift-tube induced oscillating electric fields |96, 97 1. Refer to Figure 2.1 for a schematic 
diagram that depicts the Linac and the other accelerators. 

At the end of the Linac stage, the ions are guided into the Booster, an alternating 



gradient synchroton |95] with a circumference of ~470 m, in which they make ~ 16 000 
revolutions and acquire an energy of 8 GeV. During injection into the Booster, both electrons 
are stripped from the ions by passing the ions through a carbon foil. The 8-GeV 
protons are injected from the Booster into the ~6.28-km circumference Main Ring proton 
synchroton where they are accelerated to typical energies of ^^150 GeV under the guidance 
of 1014 conventional copper-coiled dipole and quadrupole magnets [p^]. Please refer to 



Figure |2J . 

The Tevatron, located directly below and in the same tunnel as the Main Ring, is another 
synchrotron, but one that is distinct from the others in that it has magnet coils consisting 
of superconducting niobium-titanium (Nb-Ti) alloy filaments embedded in copper instead 



of the conventional copper coils used in the Booster and Main Ring magnets [P9| , IOC ] . The 
increased magnetic fields produced by the Nb-Ti magnets enable the Tevatron to accelerate 
protons to energies of nearly 1 TeV. 

Antiprotons are produced using 120-GeV protons that are extracted from the Main Ring 
and are made to strike a 7-cm thick nickel or copper target. A liquid lithium lens focuses 
the antiprotons and directs them to the Debuncher, which is a ring 520-m in circumfer- 
ence where the antiproton beam aperture and energy distribution are reduced by means 



of stochastic cooling |101] and debunching [102| techniques, respectively. The antiprotons 
are then transferred to the Accumulator ring, which is concentric with the Debuncher, for 
storage and further cooling. Once enough antiprotons have been accumulated and the Teva- 
tron has already been filled with 150-GeV proton bunches, antiprotons in the Accumulator 
are 'shot' into the Main Ring, boosted to 150 GeV, injected into the Tevatron in counter- 
rotation to the proton bunches, and then accelerated along with the protons to 900 GeV. 
The counter-rotating beams are made to collide at interaction regions such as B0 (shown in 
Figure [2T|) where 1.8- TeV collisions occur near the geometric centre of the CDF detector. 
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DO detector 



Figure 2.1: Schematic diagram of the ~6.28-km circumference Tevatron pp colhder and 
its affihated accelerators. For simphcity, the Main Ring and Tevatron are diagrammed as 
coplanar. This figure is not drawn to scale. 
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2.2 An Overview of the CDF Detector 

The Collider Detector at Fermilab (CDF detector) is a general-purpose device designed to 
study the physics of pp collisions at centre-of-mass energies near 2.0 TeV. A comprehensive 



description of the CDF detector and its subsystems is given in References |103, 104, 105 



106| , 38, 107 1 and citations therein. 



The basic design goals of the CDF detector, pictured in Figure 2^, were to identify 
leptons and measure the momenta and energies of particles originating from the B0 in- 
teraction region. Since the phase space for high energy hadronic collisions is typically 



described by rapidity (refer to Section 1.1.3), transverse momentum, and azimuthal angle, 
it is natural that the CDF detector components have an approximately cylindrical symme- 
try and uniform segmentation in pseudorapidity and azimuthal angleQ. Tracking detectors, 
which detect charged particles and measure their momenta, reside nearest the interaction 
region and inside a ~1.4-T magnetic field. The field is generated by a large electromagnet 
that consists of 1 164 turns of Nb-Ti/Cu superconductor that constitute a solenoid 4.8-m in 
length, 1.5-m in radius, and 0.85 radiation lengths in radial thickness. The tracking systems 
surround an evacuated beryllium beam pipe that is 3.8-cm in diameter, has walls 0.5-mm 



thick, and forms part of the Tevatron. Section 2.3 describes the CDF tracking systems in 
some detail. 

The detector is divided into a central region (|?7| < 1.1), two end plug regions (1.1 < 
\ri\ < 2.4), and two forward-backward regions (2.2 < \r]\ < 4.2). Refer to Figure |2.3| for 
a schematic elevation-view of one quadrant of the CDF detector. The tracking volume 
and solenoid are surrounded by sampling calorimeters that measure electromagnetic and 
hadronic energy flow from the collision point for particles with \ri\ < 4.2. The calorimeter 
systems are segmented into projective rj-ip 'towers', each of which points back towards 
the nominal interaction region and has an electromagnetic shower counter in front of a 
corresponding hadronic calorimeter cell. In high energy pp collisions, such a projective 
tower geometry is motivated by the importance of reconstructing jets, which are defined as 



^The CDF coordinate system is right-handed with x pointing out of the Tevatron ring, y vertical, and 
2 in the proton beam direction. The polar angle, 9, is measured with respect to the proton direction; the 
pseudorapidity, rj, is defined by 77 = — ln(tan(6l/2)), with 9 measured assuming a z-vertex position of zero; 
the azimuthal angle is represented by i/p and defined with respect to the plane of the Tevatron; and the 
transverse displacement coordinate is denoted by r. 
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Figure 2.2: An isometric view of three quarters of the CDF detector. 

cohections of particles that have similar trajectories in rj-ip space and are typically assumed 
to originate from a single high energy quark or gluon. 

Each of the three main CDF detector regions has an electromagnetic calorimeter; these 
are the central electromagnetic (CEM), the plug electromagnetic (PEM), and the forward 
electromagnetic (FEM) calorimeters. Behind the CEM there are two hadronic calorimeters: 
the central (CHA) and wall hadronic (WHA) calorimeters. The PEM and FEM have a 
corresponding plug (PHA) and forward hadronic (FHA) calorimeter behind each of them, 
respectively. Table 12^1] summarizes some selected properties of the CDF calorimeter systems. 

Proportional chambers located between the solenoid and the CEM constitute the central 
preradiator detector, which samples the early r-ip development of electromagnetic showers 
induced in the material of the superconducting solenoid coils. Other proportional chambers 
with strip and wire readout are located inside the CEM calorimeter at a depth of ~6Xo, 
the approximate point where an electromagnetic shower is most fully developed. These 
central electromagnetic strip detectors measure the positions of showers in both the z and 
r-ip views. 
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Figure 2.3: A schematic side-elevation view of one quadrant of the CDF detector. The 
CDF detector, with the exception of the central muon upgrade and extension subsystems, 
is largely cylindrically symmetric about the interaction region, which is located on the 
beamline at the right side of the figure. 

2.3 The Tracking Systems 



The reconstruction of exclusive S-meson decays relies heavily on precise measurements 
of the daughter particle decay vertices, momenta, and charges. The CDF detector's main 
tracking capabilities consist of four distinct but complementary tracking subsystems. These 
systems, listed in order of increasing distance from the interaction region, are the silicon 
vertex detector, the vertex time projection chamber, the central tracking chamber, and the 



central drift tube array. Figure |2.3| illustrates the positions of the tracking subsystems, both 
relative to each other and to the rest of the CDF detector. 



2.3.1 The Silicon Microstrip Vertex Detector 



The silicon microstrip vertex detector (SVX) [|I04| , 105 ] enables the identification in the r-cp 
plane of secondary vertices displaced from the pp collision point resulting from the weak 
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Property 


CEM 


CHA 


WHA 


PEM 


PHA 


FEM 


FHA 


\v\ 


- 1.1 


0-0.9 


0.7-1.3 


1.1 - 2.4 


1.3 - 2.4 


2.2-4.2 


2.3-4.2 


Ari 


~ 0.1 


~ 0.1 


~ 0.1 


0.09 


0.09 


0.1 


0.1 


Aip 


15° 


15° 


15° 


5° 


5° 


5° 


5° 


Active 


plastic scintillator 


gas chambers with cathode pad readout 


Absorber 


Pb 


Fe 


Fe 


Pb 


Fe 


Pb/Sb 


Fe 


Thickness 


18Xo 


4.5Ao 


4.5Ao 


18 - 21Xo 


5.7Ao 


25Xo 


7.7Ao 



Table 2.1: Selected properties of the CDF calorimeter systems. Shown are the pseudorapid- 
ity coverage and segmentation, the azimuthal segmentation, the active medium, the type 
of absorber, and the thickness in radiation lengths (Xq) and interaction lengths (Aq) of the 
electromagnetic and hadronic calorimeters, respectively. 



decays of b quarks. Installed in the CDF detector in 1992, the SVX was the first detector 
of its kind to be operated successfully in a hadron collider environment. In 1993, a more 



radiation-hard and low- noise version of the SVX, the SVX' |10(:], was commissioned for the 
1994-1995 Tevatron collider run|^ (refer to Section 3.1). The SVX consists of two identical 
cylindrical modules, one of which is pictured in Figure each comprising four concentric 
cylindrical layers with radii of 3.0, 4.3, 5.7, and 7.9 cm. The SVX' has the same overall 
configuration as the SVX, except that the innermost layer has a slightly smaller radius of 
2.9 cm. Since the luminous pp interaction region is rather elongated in the z direction (with 
a Gaussian distribution having a standard deviation of o" ~ 30 cm), approximately 40% of 
pp collision vertices lie outside the acceptance of the SVX, which has an active length of 
51.1 cm. 

In both of the SVX barrels, the four layers are each segmented into twelve 'ladders' that 



subtend approximately 30° in azimuth and are oriented parallel to the beam axis. Figure 
depicts a typical ladder situated in the third layer of the SVX. The ladder substrates are 
fabricated from a light-weight foam (Rohacell) reinforced with strips of carbon fibre. Three 
single-sided 8.5-cm long silicon microstrip detectors are electrically bonded together with 
aluminum wire along the z direction to form a 25.5-cm active silicon region on each ladder 



^Unless noted otherwise, references to the SVX apply to the SVX' as weU. 
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Figure 2.4: Isometric view of one of the two silicon microstrip vertex detector (SVX) barrels. 
The dummy-ear sides of both barrels are conjoined (with an effective gap of 2.15 cm) at 
the z = position inside the CDF detector. 

module. 

Silicon microstrip detectors are a kind of solid ionization chamber in which incident 
charged particles dislodge electrons via ionization. In their most basic form, the detectors 
consist of an n-doped silicon wafer, typically 300-//m thick, with strips of p-type material 
on one side |108| ] . If a reverse-biased potential is applied to a strip, a p-n junction diode is 
set up, and electron-hole pairs created by the passage of a charged particle migrate to their 
respective electrodes, thereby manifesting an electronic signal that resolves the location of 
the particle's trajectory as a function of the strip separation, or 'pitch'. 

The silicon strip pitch of the inner three SVX layers is 60 ^m and that for the outermost 
layer is 55 fim. The average position resolutions for the SVX and the SVX' were measured 
to be 13 fim and 11.6 fim, respectively, and the high transverse momentum (asymptotic) 
impact parameter resolution was determined to be 17 /zm for the SVX and 13 ^um for the 
SVX'. Adjacent ladders in a given layer slightly overlap each other to provide full azimuthal 
coverage; this is achieved with a 3° rotation of the ladders about their major axes. The 
SVX' has all four of its layers overlapped; however, the innermost layer of the SVX suffers 
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Figure 2.5: Layout of a typical ladder module in the third layer of the SVX. Three single- 
sided silicon microstrip detectors are wire bonded together to constitute each ladder module. 



from a 1.26° gap in ip between adjacent ladder modules. 

As shown in Figure |2.5| , the outside end of each ladder has a small circuit board that 
contains the front end readout chips, which each serve 128 channels. Because the ladder 
widths increase with increasing r, the number of readout chips on a given ladder module 
depends upon the layer in question. The innermost layer, for example, has two readout 
chips per ladder module whereas the outermost layer has six chips per ladder module. The 
total number of instrumented strips in the SVX is 46 080. 

The silicon strips used in the SVX' are AC coupled whereas those for the SVX are DC 
coupled. The SVX' consequently benefits from a marked reduction in noise compared to the 
SVX, for which strip-to-strip leakage current variations have to be subtracted by cycling 
the front end readout chips through two successive charge integrations. An additional 
advantage of the AC-coupled SVX' design is that the readout chip preamplifiers, which 
have 40% more gain than their SVX counterparts, will not saturate, even after significant 
radiation damage has increased the silicon microstrip leakage currents appreciably. The 
SVX' readout chips were fabricated using radiation hard 1.2-^m CMOS technology. They 
therefore have a design absorbed-dose limit of 10 kCy compared with 200 Gy for the SVX, 
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which uses more conventional 3-/im CMOS electronics [109|. 

The readout electronics typically generate ~53.1 W of heat in each of the two SVX 
barrels. Cooling pipes transport chilled de-ionized water at a temperature of 13 °C and a 



flow rate of 10 g/s to the beryllium bulkhead (see Figure 2A) and the readout circuit boards 



(see Figure 2^) to maintain an operating temperature near 20 °C. The cooling circuit runs 
at a subatmospheric pressure to minimize the potential damage due to a breach in the 
cooling pipes. Controlling the temperature not only minimizes leakage currents in the silicon 
microstrips and prevents damage to the front end electronics, but it also discourages thermal 
gradients in the mechanical support structure that can distort the internal alignment of the 
SVX. 

2.3.2 The Vertex Time Projection Chamber 

A vertex time projection drift chamber (VTX) surrounds the SVX (refer to Figure 2.3). It 



was designed to measure the trajectories of charged particles in the r-z plane in a pseudora- 
pidity range \rj\ < 3.0. The VTX is important for the determination of the z position of the 
primary vertex and the identification of multiple interactions in the same beam crossing. 



The VTX resolution of a primary vertex location along the beamline, nominally 2 mm [ 110 |, 
depends on the number of detected tracks originating from that location and the number 
of primary pp interactions in the event. 

The VTX, which extends 132 cm in each z direction and has a radius of 22 cm, consists 
of 28 drift modules, each containing two drift regions separated by an aluminum central 
high voltage grid. Endcaps on each side of the drift modules are azimuthally segmented into 
octants and are rotated in (p by 15° with respect to adjacent modules (in z) to eliminate 
inefficiencies near module boundaries. Within each octant, 16 or 24 sense wires, depending 
on the z position of the module, are oriented tangentially, thereby providing tracking infor- 
mation in the r-z view. The z location of a track with respect to a given wire in a given 
module is determined by the drift time, and the r information is determined from the radial 
location of the wire. The electric field is maintained near 1.6 kV/cm and the gas used is 
a 50:50 admixture of argon and ethane. Figure p.6| is an event display diagram showing 
hits from charged particles in two VTX octants. The VTX vertex-finding algorithm has 
calculated a primary vertex z position in this event based on a fit to extrapolations of tracks 



2.3. THE TRACKING SYSTEMS 31 



reconstructed from sense wire hits. The z vertex position is represented by the seriffed cross 



in Figure 2.6. 



2.3.3 The Central Tracking Chamber 

The most prominent subsystem in the CDF detector is the central tracking chamber, or 
CTC. It is the only CDF tracking device that can perform three dimensional momentum 
and position measurements, both of which are essential to the reconstruction of exclusive 



S-meson decays. The CTC, as indicated in Figure p.3| , surrounds the VTX and SVX 
subsystems and has a coaxial bicylindrical geometry with a 3 201.3-mm length (including 
the endplates), a 2 760.0-mm outer diameter, and a 554.0-mm inner diameter. Aluminum 
is used in the construction of the outer cylinder; carbon fibre reinforced plastic constitutes 
the inner cylinder wall. 

The CTC is a drift chamber that contains 84 layers of 40-//m diameter gold-plated 
tungsten sense wires arranged into nine 'superlayers', five of which have their constituent 
sense wires oriented parallel to the beam axis (axial superlayers), and four of which have 
their wires canted at angles of either +3° or —3° with respect to the beamline (stereo 
superlayers). The innermost and outermost sense wires have radii of 309 mm and 1 320 mm, 
respectively. The axial and stereo superlayers alternate with increasing radius and each 
consists of twelve and six sense wire layers, respectively. The configuration is illustrated in 



Figure 2.7, which shows the wire slot locations in the aluminum endplates. The majority 
of the CTC pattern recognition is done using data from the axial layers, which provide 
tracking information in the r-ip view. The stereo layers furnish tracking information in the 
r-z view. 

The superlayers are functionally segmented into open drift cells. A drift cell contains 
a superlayer of (either 12 or 6) sense wires alternating with (either 13 or 7) stainless steel 
potential wires, which serve to control the gas gain on the sense wires. Two planes of 
stainless steel field wires on either side of the sense wire superlayers define the fiducial 
boundaries of each drift cell and control the strength of the electric field in the <40-mm 
drift regions. The number of cells in each superlayer increases with r such that the drift 
distance, which translates into a maximum drift time of ~800 ns, is approximately constant 
for all cells in the CTC. To keep the electric field uniform throughout the fiducial volume 
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Figure 2.6: Event display of two 45° octants in the vertex time projection chamber (VTX). 
The beamhne runs through the middle of the diagram, across the page. Shown are the 28 
drift modules containing hits that form tracks due to the traversal of charged particles from 
the interaction region. The hits represent the times that the drifting ionization electrons 
arrive near the sense wires. The crosses along the beamline represent z vertex candidates 
found by the VTX software; the single seriffed cross represents an especially high quality 
vertex that consequently defines the zero location of the event pseudorapidity coordinates, 
as indicated by the r\ scale denoted around the outside of the figure. The two rectangles 
drawn near the geometric centre of the VTX represent active regions of the SVX. This event 
is a dramatic example of a case where the primary vertex is so distant from z = that the 
SVX is of little consequence. Also interesting is the helical trajectory of a particle in the 
upper VTX octant at ~ —3.5, a result of the large axial magnetic field. 
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Figure 2.7: One aluminum endplate of the central tracking chamber (CTC), viewed from 
along the beam axis. The wire slot locations for the alternating axial and stereo superlayers 
are apparent. 
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of every drift cell, extra stainless steel shaper and guard wires are located near the cell 
perimeters, bringing the number of wires in the CTC to a total of 36 504. This translates 
to a total wire tension of 245 kN and a combined wire length that is in excess of 110 km. 



As is evident in Figure |2.7| , the CTC drift cells are tilted such that the angle between 
the radial direction and the electric field direction is approximately 45°. Such a large 
cant angle is necessary to offset the 45° Lorentz angle, which results from the combined 
effects of the ~1.4-T magnetic field, the argon-ethane-alcohol gas mixture used (in the 
proportions 49.6% : 49.6% : 0.8%), and the relatively low ~1.35-kV/cm electric field. The 
drift trajectories in the CTC are therefore approximately parallel to the azimuthal direction. 
Every high transverse-momentum track passes close to at least one sense wire in each 
super layer. 

Preamplifiers mounted on the endplates of the CTC are connected to the sense wires, 
whose analog signals are amplified further, shaped, and discriminated by circuitry mounted 



on the solenoid return yoke (see Figure O). The discriminator signals undergo time-to- 
digital (TDC) conversion in a counting room located at the end of 70 m of fiat cable. The 1- 
ns resolution TDCs can record > 7 hits per wire per event. The CTC double track resolution 
is < 5 mm due to the approximately 100-ns minimum separation between two resolved hits. 
The CTC has a single hit resolution of < 200 nm, and the overall momentum resolution of 
the combined SVX-CTC system is Spt/pt = (0.0009 px) + (0.0066) , where pT is the 
transverse momentum measured in units of GeV/c. Section |5.2| and Appendix ^ provide 
further discussion of the performance of the CTC, and Figure |3.2| contains an event-display 
diagram of a sample event showing reconstructed CTC track candidates. 

2.3.4 The Central Drift Tube Array 

The central drift tube array, or CDT, is situated at a radius of 1.4 m, between the outer 
cylinder of the CTC and the inner wall of the solenoid cryostat, as indicated in Figure |2.3| . 
The CDT system consists of stainless steel circular tubes; these are 1.27-cm in diameter, 
3-m in length, and 2 016 in number. Closely packed into three layers, the tubes are each 
strung with 50-/xm diameter stainless steel anode wires. By virtue of its charge division 
capability on the anode wires, the CDT can provide tracking information in both the r-ip 
and r-z views. For the analysis described in this thesis, CDT tracking information was 
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not used explicitly in the reconstruction of particle tracks; however, the CDT was used to 
identify cosmic ray muons as coincident hits with Aip ~ 180°. Cosmic ray muons were used 
to perform the initial relative alignments of the SVX, VTX, and CTC subsystems within 
the CDF detector. 

2.4 The Muon Chambers 

The ability to identify muons and their trajectories is essential to the reconstruction of J ftp 
and ip{2S) mesons in the dimuon channels. Muon identification can be achieved by exploit- 
ing the relatively high muon critical energy^, which is several hundred GeV in materials 
such as iron [^], significantly higher than the critical energy for other ionizing particles. 
This ability of the muon to penetrate matter thus motivates the location of the muon 
subsystems in the outer regions of the CDF detector that can only be reached by those 
charged particles that originate from the interaction region and that penetrate the inter- 
vening material. This material, consisting primarily of the calorimeters, serves to filter out 
the majority of hadrons and electrons before they reach the muon subsystems. Refer to 



Figure 2.3 for the locations of the three central muon subsystems used in this analysis: the 
central muon detector (CMU), the central muon upgrade detector (CMP), and the central 
muon extension (CMX). A map of the r]-(p muon detection coverage in the central region is 



shown in Figure |2.8| . The forward muon toroid subsystem, shown in Figures and |2.3| , is 
not used in this study due to its poor intrinsic momentum resolution and the lack of overlap 
in acceptance between it and the CTC and SVX tracking systems. 

2.4.1 The Central Muon Detector 

The CMU covers the region 55° < |^| < 90° and resides on the outer edge of the central 



hadronic calorimeter, 3 470 mm from the beam axis, as indicated in Figure 2.9. Each 12.6° 



azimuthal wedge comprises three modules, each subtending 4.2° in ip. A CMU module. 



shown in Figure 2.1C, consists of four towers, each with four layers of rectangular drift cells. 
The outermost and second innermost cells in each tower are oriented such that their sense 
wires lie on a radial that originates from the centre of the CDF detector. The innermost 
and second outermost drift cells lie on another radial that is offset from the first by 2 mm 



■^The muon critical energy is the energy at which losses due to radiation and ionization are equal |30 
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Figure 2.8: Coverage of the central muon subsystems in pseudorapidity (ry) and azimuth 
{ip) The lack of CMX coverage at 99 ~ 90° and (p ~ 270° results from the interference 

due to the Main Ring bypass beampipe and the concrete collision hall floor, respectively. 



2. 4. THE MUON CHAMBERS 



37 




Figure 2.9: The layout of a central muon (CMU) wedge with respect to a central calorimeter 
wedge in both the azimuthal (left figure) and polar (right figure) views. 



at the midpoint (in r) of the CMU. The offset cells in each tower resolve the side of the 



radial, in azimuth, on which the track passed. As indicated in Figure 2.1C, the absolute 
difference in drift electron arrival times for a pair of cells having sense wires on the same 
radial determines the angle between the candidate muon track and that radial. This angle 
can be related to the transverse momentum of a muon candidate and is therefore exploited 



by the trigger system (refer to Section 2.5.1 ). 



A drift cell in the CMU, shown in Figure 2.11 , is rectangular with dimensions 63.5 mm 
X 26.8 mm x 2 261 mm and has a single 50-^m stainless steel sense wire strung through 
its centre. The drift cells are operated in limited streamer mode using a 50:50 admixture 
of argon and ethane gas, and potentials of +3 150 V on the sense wires and —2 500 V on 
the I-beams, which are electrically isolated from the top and bottom aluminum plates by 
0.62 mm of insulation. The position of a muon candidate track along the sense wire (z) 
direction can be determined with a resolution of 1.2 mm using charge division electronics. 
The position resolution in the drift ((/?) direction is 250 /xm. 
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To pp interaction vertex 



Figure 2.10: Layout of a central muon detector (CMU) module showing the four towers, 
each with four layers of rectangular drift cells. The outermost and second innermost cells 
in each tower are oriented such that their sense wires lie on a radial that originates from 
the geometric centre of the CDF detector. The other two drift cells are offset to determine 
which side of the radial the track passed. The quantities t2 and t4 represent drift electron 
arrival times; their difference, |t4 — t2\, determines the angle between the candidate muon 
track and the radial, thus providing a crude but fast measurement of transverse momentum 
that can be used in a low level trigger. Analogous information from ti and yields a second 
independent measurement. 



29.0 mm 




Figure 2.11: Layout of a central muon detector (CMU) drift cell, showing the 0.79-mm 
aluminum top and bottom plates and the aluminum I-beams that separate adjacent towers. 
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2.4.2 The Central Muon Upgrade Detector 

An average of 5.4 pion interaction lengths lies between the CMU and the pp collision re- 
gion, resulting in approximately 1 in 220 high energy hadrons traversing the calorimeters 
unchecked. This 'noninteracting punch-through' results in an irreducible false muon back- 



ground rate. The central muon upgrade detector (CMP), shown in Figures p]3 and 2.3, was 



commissioned to contend with this punch-through hadron rate [111]. The CMP surrounds 
the central region of the CDF detector with 630 tons of additional steel. The geometry 
is box-like, with the return yoke of the solenoid providing the absorption steel on the top 
and bottom, and two retractable 60-cm thick slabs arranged as vertical walls on each side. 
The additional absorption material brings the number of pion interaction lengths to 7.8 at 



9 = 90°. Figure 2.8 illustrates the variation in pseudorapidity coverage with azimuth caused 
by the geometry of the CMP. 

The active planes of the CMP consist of four layers of half-cell staggered single- wire drift 
tubes operating in proportional mode. Each drift cell, of which there are 864 in the CMP, 
consists of a rectangular extruded aluminum tube 25.4-mm high, 152.4-mm wide, and with a 



length that depends upon where the tube is mounted. Figure 2.12 is a schematic drawing of 
the components of a CMP proportional drift cell. The anode, a 50-//m gold-plated tungsten 
wire, is biased to a potential of -1-5 600 V, the wide central field-shaping cathode pad is 
biased to -1-3 000 V, and the eight narrow field-shaping strips have decrementally decreasing 
voltages from the centre of the cell out to the edges. The maximum drift time is 1.4 fis. 

2.4.3 The Central Muon Extension 

The central muon extension, or CMX, provides additional pseudorapidity acceptance in 



the region 0.65 < \rj\ < 1.0. Shown in Figures |2.2| and the CMX modules possess 
geometries that correspond to sides of the frusta of two cones, each with a base at z = 
and an axis along either the proton or antiproton direction. The azimuthal coverage of the 
CMX is not continuous; due to the floor of the collision hall, there is a 90° gap in ip at the 
bottom of the CDF detector, and, due to the Main Ring bypass beampipe, there is a 30° 
gap at the top of the detector. The 1 536 proportional drift cells that constitute the CMX 



modules are shorter than, but otherwise identical to, those used in the CMP (Figure 2.12). 



No additional absorber was added between the CMX and the interaction region; however, 
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Figure 2.12: Schematic drawing of a proportional drift cell used in both the central muon 
upgrade detector (CMP) and the central muon extension (CMX) ]111| ]. 
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Figure 2.13: Proportional drift cell layout in a 15° module of the central muon extension 
(CMX) (top figure). Also shown is a close-up view of the staggered cell configuration about 
one of the threaded rods used to fix the CMX to its support arch (bottom figure) |111|]. 



the smaller polar angle through the hadronic calorimeter and magnet return yoke yields a 
shielding thickness of 6.2 pion interaction lengths at 6 = 55°. 

The CMX is organized into four stacks, two on the proton side and two on the antiproton 
side of the CDF detector. Each stack consists of eight modules, which each subtend 15° in 
azimuth. A module comprises 48 proportional drift cells that are grouped in eight half-cell 



staggered layers of six tubes each. Refer to Figure 2.13 for an illustration of the interleaved 
geometry necessary to arrange the rectangular cells along a conical surface. The invisibility 
of intermodule boundaries is an advantage of this interleaved configuration. 

The maximum CMX drift time, 1.4 fis, is such that the spread of arrival times due to 
background particles is short by comparison. Background rejection and a high speed trigger 
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are provided by an array of scintillation counters mounted on the inner and outer sides of 
each CMX module. Four such scintillators, each with its own photomultiplier tube, are 
located on both sides of every 15° CMX module, for a total of 256 scintillation counters. 
Background effects are vetoed in the trigger by requiring that both the inner and outer 
scintillators adjacent to a CMX hit produce pulses that are coincident with the pp beam 
crossing to within a few nanoseconds. 



2.5 The Trigger Systems 

In the course of data acquisition for this study, proton and antiproton bunch crossings in the 
Tevatron collider occurred every 3.5 /xs, corresponding to a crossing frequency of 286 kHz. 
With typical instantaneous luminosities of £ = 10'^^ cm~^s^^ and a measured pp total cross 



section of utot = 80.03ib2.24 mb | 112 | at y/s = 1.8 TeV, at least one pp interaction per beam 
crossing was expected. Moreover, the low signal-to-noise ratio of the bb hadroproduction 
processes (refer to Section |1.1.3| ) made the implementation of a trigger system essential. 

A CDF event, which amounted to the digitized information from a single beam crossing 
that could be read out from the CDF detector at a given time, had a data length of ~165 kB. 
Such an event size could only be reliably written out to several 8-mm magnetic tapes at 
a rate of approximately 10 Hz. This constituted the principal limitation to the CDF data 
acquisition rate and necessitated a trigger system that could both accommodate the pp 
interaction rate and select interesting physics events with a ~30 000:1 rejection factor. The 



CDF trigger |103| , 113 1 consisted of three successive levels, each of which imposed a logical 
"or" of a limited number of programmable selection criteria that collectively reduced the 
data rate exposed to the next higher trigger level. The reduction in the rate presented to 
the higher trigger levels provided time for more sophisticated analyses of potential events 
with the accrual of less dead timeQ. 

2.5.1 Level 1 

The Level 1 trigger required less time than the 3.5-/is beam crossing period to reach a 
decision on whether or not a given event was suitable for consideration by the higher trigger 



^In this context, 'dead time' refers to the amount of time that the CDF detector was unable to consider 
subsequent pp colhsions. 
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levels; it therefore incurred no dead time. Such speed of operation was achieved by the 
analog read-out and processing of data from selected detector components with dedicated 
FASTBUS-based electronics. At an instantaneous luminosity of £ = 5 x 10^'^ cm^^s^^, the 
Level 1 trigger acceptance rate was approximately 1 kHz Although it could have been 
configured to base its decisions on information from several different CDF subsystems, the 
Level 1 trigger primarily used signals from the calorimeters and the muon systems. 

The calorimetry component of the Level 1 trigger considered analog signals on dedicated 
cables from the scintillator photomultiplier tubes in the central calorimeter subsystems and 



from cathode pads in the plug and forward calorimeter subsystems (refer to Table 2T). For 
the purposes of the trigger, the calorimeters were logically segmented into 'trigger towers' 



with Aip = 15° and Arj = 0.2. For each subsystem listed in Table |2.1| , the individual tower 
minimum-energy thresholds could be specified to the trigger, which summed the deposited 
energies, weighted by the polar angle to determine the transverse energy E't = E sin 9, 
for all those trigger towers that were above these thresholds. If the total i?T) measured in 
this manner, exceeded a given global threshold, then the Level 1 trigger accepted the event. 
There was also a similar Level 1 calorimetry trigger that had significantly reduced tower 
energy thresholds, but was prescaled by a factor of 20 or 40, depending on the data-taking 
period (Run lA or Run IB, respectively; refer to Section pTip . 

The muon component of the Level 1 trigger exploited the relative drift electron arrival 
times (At) between pairs of drift cell layers in a given CMU module, as described in Sec- 
tion |2.4.1| . The two cells constituting each of these pairs were separated in r by one drift 



cell, as shown in Figure 2.1C. The trigger logic operated on objects, called 'muon stubs', 
that were defined by the existence of any wire pair in a 4-tower 4.2° muon detector module 
(see Figure |2.10| ) with a At less than a value corresponding to a given minimum transverse 
momentum (px) requirement. 

Out of a total of seven Level 1 triggers that involved muon candidates, two were directly 
relevant to this analysis because they specifically identified dimuon candidates in the central 
region of the CDF detector. One of these dimuon triggers (two_CMU_3pt3) required that 
two CMU muon stubs exist, whereas the other (two_cmu_CMX_3pt3) required that the 
two stubs each be in either CMU or CMX modules. In the case of the CMX stubs, the 



trigger required coincident hits in the CMX scintillators (see Section 2.4.3). The Level 1 
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system made no further requirements on the positions of the two muon stubs, except for a 
criterion that they be located in noncontiguous modules. That is, at least one muon module 
lacking a muon stub must have resided between the two modules where muon stubs were 
observed. If this was not the case, then the two adjacent stubs were merged into a single 
muon stub for the purposes of the Level 1 trigger. The minimum requirement on each 



stub in these two triggers was nominally 3.3 GeV/c. Section 4.4 discusses the efficiencies of 
the Level 1 low-px central muon triggers. 

2.5.2 Level 2 

In a pp beam crossing for which the Level 1 trigger did not fire, a timing signal from the 
Tevatron announcing the occurrence of the next beam crossing would cause the stored sig- 
nals in the CDF detector to be cleared in preparation for the next crossing. If the Level 1 
trigger did fire, then subsequent timing signals were inhibited from clearing information 
stored in the CDF detector for a period of up to 20 ^s, during which the Level 2 trigger 
made its decision and five disregarded beam crossings could occur. At an instantaneous lu- 
minosity of £ = 5 X 10'^'' cm~^s~^, the typical Level 2 trigger output rate was approximately 
12 Hz ||. 

With the increased processing time, the Level 2 trigger system could perform simple 
tracking calculations and determine basic topological features of the event by consider- 
ing, with greater sophistication, the same dedicated calorimetry and muon signals used in 
Level 1. Specifically, expeditious electromagnetic and hadronic transverse-energy clustering 
was performed at Level 2 by applying 'seed' and 'shoulder' thresholds to all the calorimeter 
trigger towers. If a given tower energy exceeded the seed threshold, which is higher than 
the shoulder threshold, then a cluster was formed. Adjacent trigger towers were iteratively 
included in this cluster if they both exceeded the shoulder threshold and were not part of 
another cluster. 

High speed track pattern recognition was achieved in Level 2 with the central fast tracker 
(CFT), a hardware track finder that detected high-px charged particles in the CTC (Sec- 
tion |2.3.3| ). The CFT measured transverse momentum and azimuth, since it only examined 
hits in the five axial CTC superlayers. For a given traversal of an axial superlayer by a 
charged particle, the CFT considered two types of timing information: prompt and delayed 
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hits. Prompt hits, gated <80 ns after the beam crossing time, were due to the short drift 
times caused by charged particles traversing the plane of sense wires in a superlayer. Pairs 
of delayed hits, one on each side (in ip) of a given superlayer, were recorded by a gate that 
occurred 500-650 ns after the beam crossing. The absolute prompt and delayed drift times 
provided information on a track's trajectory, whereas the relative drift times furnished mea- 
surements of curvature, and hence p^. After all the drift hits were recorded, the CFT sought 
to construct tracks by first examining hits in the outermost superlayer. For each sense wire 
in the outer layer with a prompt hit, the CFT looked to the inner layers for 'roads', or 
hit patterns, that matched patterns in a look-up table that had eight pT bins and two ip 
bins, one for each sign of curvature. The pT bins ranged^ from ~3 to ~30 GeV/c, and the 
transverse momentum resolution was Spi/pT ~ 0.035 x p^, where px is in units of GeV/c. 

The Level 2 trigger system organized the energy clusters, CFT tracks, and muon stubs 
into clusters called 'physics objects'. These included jets, SE't (total transverse energy), 
electrons, photons, taus, muons, and neutrinos (whose signature is missing transverse en- 
ergy, pr)- For the majority of the data-taking period, a custom-built 'Jupiter' module 
accessed the clusters and made a Level 2 decision to accept or reject events. The Jupiter 
module had two separate processor boards, one that loaded clusters into memory, and an- 
other that checked the triggers by imposing several requirements on the physics objects. In 
the last year of data acquisition, these processors were replaced with commercially available 
AXP (Alpha) processors manufactured by Digital Equipment Corporation. 

All of the dimuon selection triggers at Level 2 imposed a matching criterion between at 
least one of the two Level 1 muon stubs and a CFT track. Early in the data-taking period, 
the requirement was that the stub and the extrapolated track have an azimuthal separation 
in the transverse plane that was <15°. This criterion was later tightened to A(p < 5° to 
reduce further the trigger rate due to accidental coincidences. The various dimuon Level 
2 triggers used in this analysis are listed, along with their prerequisite Level 1 triggers, 
in Appendix Section iA discusses the efficiencies of the Level 2 low-px central muon 
triggers. 



^The CFT px-bin thresholds were changed between Runs lA and IB. 
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2.5.3 Level 3 



The Level 3 trigger |114] was a flexible, high-level, software-based computer processor 'farm' 
that could reconstruct several events in parallel. When the Level 2 trigger accepted an event, 
the channels in the CDF detector with valid data were digitized and read out by the data 
acquisition (DAQ) system. The DAQ electronics subsequently transported the event data 
to the Level 3 processor farm. Over the course of the data-taking period, both the Level 3 
trigger system and the DAQ system (see Section [2^ ) underwent several significant changes. 
Although most of these changes were effected in the interval between the Run lA and Run 
IB collider running periods, not all of them were implemented for physics data taking from 
the beginning of Run IB. Unless specifically noted otherwise, the following descriptions of 
the Level 3 trigger and the DAQ system pertain to the upgraded configurations used for 
data taking later in Run IB. 

The Level 3 computing farm consisted of 64 commercial processors that were manufac- 
tured by Silicon Graphics, Inc. and that ran under the IRIX operating system, a flavour of 
UNIX. Half of these processors were R3000 Power Servers and half were R4400 Challenge 
machines. The farm processors received data fragments read out by the DAQ system for a 
given beam crossing and 'built' these fragments into a contiguous event. Prior to the logging 
of the events to disk or 8-mm magnetic tape, the processors reconstructed and characterized 
these events for later selection using a configurable trigger table and optimized executables 
of compiled FORTRAN computer codes. A given event data buffer did not necessarily re- 
side on the same farm computer as the reconstruction executable and processor that were 
operating on it. 

In Run lA, only 48 R3000 Power Servers were used, and the event fragments were 
already built into events by the Run lA DAQ system prior to their reception by the Level 3 
farm (refer to Section ^.6| ). Every processor had two dedicated, but separate, buffers, which 
each had the capacity to contain an entire event. A separate 'farm steward' computer 
communicated with the rest of the DAQ system, controlled the initiation and cessation of 
event processing on the farm CPU^, maintained performance statistics of farm activities, 
and provided Level 3 status information. In the Run IB trigger system, however, the duties 



^CPU is an acronym for 'central processing unit' 
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of the farm steward were absorbed into other programmes executing on the farm computers. 

For the purposes of analyses involving central muons, the Level 3 executables recon- 
structed muon stubs and CTC tracks using algorithms that were largely identical to those 
employed in the off-line event reconstruction (refer to Section 3.2); however, because three- 



dimensional track reconstruction constituted most of the Level 3 execution time, only the 
faster of two tracking algorithms used in the off-line code was engaged in the trigger. Two 
Run lA dimuon triggers were used to form the data samples for the present analysis. 
One of these, that which contained dimuon candidates from the decay J/ip ^ fi^ , re- 
quired two oppositely charged muon candidates with a combined invariant mass in the 
range 2.8 — 3.4 GeV/c^. The other trigger, which accepted dimuon candidates from the de- 
cay ip{2S) — > fj,^ , had no opposite-charge requirement and selected dimuon candidates 
with an invariant mass in the range 2.8 — 4.0 GeV/c^. In Run IB, both the J ftp and ip{2S) 
dimuon modes were accepted by the same trigger, one that imposed no opposite-sign charge 
requirement and had an invariant mass window of 2.7 — 4.1 GeV/c^. 

In addition, the Level 3 dimuon triggers used in this study placed position matching 
requirements between the muon stubs and their associated CTC tracks. The algorithm ex- 
trapolated the CTC track to the appropriate muon subsystem and determined the difference 
in position between the projected track and the muon stub in both the r-ip plane and the 
z direction, correcting for energy loss and multiple scattering as a function of px- The Run 
lA J/ip tracks were required to match the muon stubs within 4 standard deviations (a) of 
the combined multiple scattering and measurement uncertainties, whereas the ip{2S) tracks 
had a matching requirement of ~6(7. The Run IB Level 3 muon matching requirement was 

~3(T. 



2.6 The Data Acquisition System 

The data acquisition (DAQ) pipeline, at its lowest, or 'front end', level, began with the 
readout of analog signals from channels in the various subsystems of the CDF detector. In 
general, these analog pulses were preamplified and transported to an electronics crate on 
the detector where they underwent further amplification and, depending on the subsystem, 
pulse shaping and discrimination. Digitization of these analog signals was accomplished 
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with either analog-to-digital or time-to-digital conversion electronics. 

Dedicated signal cables communicated synoptical event information from the front end 
electronics to the Level 1 (LI) and Level 2 (L2) triggers, as illustrated schematically in 



Figure 2.14. The Level 1 and Level 2 event acceptance decisions were coordinated by a 
programmable fastbus device known as FRED|^ which acted as the interface between the 
trigger and the trigger supervisor board. The trigger supervisor, also a fastbus module, 
initiated and monitored the readout of data from the front end electronics using fastbus 
readout controllers, or FRCs, which were single-width modules that contained MIPS R3000 
processors. The FRCs sent their data over a 16-bit parallel scanner bus to six scanner 
CPUs, or SCPUs. The SCPUs, which were VME^-based Motorola 68030 processors runnmg 
the VxWorks operating system, transported the event fragments to the Level 3 processor 
nodes via a commercial 256 Mbit/s serial Ultranet hub. The scanner manager, also a 
VME-based 68030 CPU, used a fibre optic reflective memory network to control the flow 
of data between the FRCs and the Level 3 trigger system, ensuring that all fragments of 
a given event were destined for the same Level 3 node. The trigger supervisor interface. 



also shown in Figure 2.14 , facilitated communication between the scanner manager and the 
trigger supervisor. It consisted of two FRCs that resided in the same crate as the trigger 
supervisor and were connected to the scanner manager via a scanner bus. 

Events accepted by the Level 3 trigger were passed via Ultranet to the consumer server 



(see Figure 2.14 ), a process running on a dedicated two-processor Silicon Graphics, Inc. 
Challenge L machine. The consumers, which received events from the consumer server over 
ethernet, were diagnostic applications that monitored trigger rates, luminosity conditions, 
detector performance, and a selected number of known high-rate physics processes. Another 
consumer, an on-line interactive event display system, provided CDF control-room personnel 



with graphical physics and detector performance diagnostics (refer to Figures |2.6| and 3.2 



for examples of two of the system's displays). Data logger programmes, one for each output 
data stream, ran on the consumer server machine and wrote accepted events to local disk. 
The disk-resident data logger events were subsequently written to 8-mm tape by a tape- 
staging programme that also ran on the consumer server computer. 



^The acronymic or abbreviational origins of FRED are unknown |113|. 

®VME stands for 'Versa Module Eurocard' and is a crate-based electronics package scheme. 
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Figure 2.14: A schematic drawing of the principal elements of the CDF data acquisition 
pipeline |115]. The individual components are described in the text. Scramnet is the fibre 
optic reflective memory network. 
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In contrast, the Run lA DAQ system lacked FRCs, using custom-built scanner modules 
instead. A hardware event builder |116|, rather than SCPUs, constructed complete events 
and communicated them to the Level 3 trigger, fastbus, in lieu of the fibre optic reflective 
memories, delivered all system control messages, and the role of the Run IB scanner manager 
was served in Run lA by a buffer manager running on a MicroVAX II computer. 



Chapter 3 

Selection of Candidate Events 



This chapter discusses the means by which charmonia and i?-meson candidates are extracted 
from the backgrounds that are characteristic of hadron-hadron cohision environments. First, 
the data set is defined. This is fohowed by a general description of the software machinery 
used to reduce the data sample. Finally, specific details about the criteria used to identify 
candidate signal events and reject background sources are presented. The charmonia and 
S-meson candidate invariant mass distributions will be presented in Chapters |5| and |^, 
respectively. 

3.1 The Data Sample 

The experimental data used in this investigation were recorded, using the CDF detector, in 
the years 1992 through 1995. During this period, the Fermilab Tevatron collided protons 
with antiprotons at a centre-of-mass energy of ^/s ~1.8 TeV. 

The 1992-1995 Tevatron running period, which is generally known as Run 1, comprised 
two separate data-taking intervals: Run lA and Run IB. Run lA consisted of approxi- 
mately nine months of physics running, which commenced on 26 August 1992 and ended on 
30 May 1993. The data collected during Run lA correspond to a time-integrated luminosity 
of jCdt = (19.5 ± 1.0) pb-^ 

In the intervening period between Run lA and Run IB, the CDF detector underwent 
several improvements. As described in Chapter |2|, the SVX was replaced with the SVX', 
new trigger requirements and hardware were implemented, and a significantly improved 
data acquisition system was brought on-line. Most components of the upgraded trigger and 
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data acquisition hardware were piggybacked onto the existing system to facihtate, if the 
need arose, the immediate reversion to their analogous Run lA counterparts. In this way, 
the AXP Level 2 trigger processors and the upgraded data acquisition system underwent 
testing early in Run IB, but were only used for physics data collection several months into 
the Run IB Tevatron running period. 

Physics-quality data collection in Run IB began on 18 January 1994 and ended on 
24 July 1995. The recorded time-integrated luminosity was JCdt = (89 it 7) pb^^, yielding 
a consolidated Run 1 total of JCdt = (109ib7) pb^^. It should be noted that some operating 
conditions were not constant over the course of the Run 1 data-taking period. The tempera- 
ture and humidity in the collision hall, for example, exhibited a strong seasonal dependence. 
The CDF solenoid, described in Chapter |2|, was operated at slightly different magnetic field 
strengths (refer to Section ^^ ). The most significant operating variations during Run 1, 
however, were those associated with the performance of the Tevatron accelerator. While the 
centre-of-mass energy remained constant at 1.8 TeV, the instantaneous luminosity delivered 
to the B0 pp interaction region varied from a mean value of £ = 3.5 x lO^'' cm~^s~^ in Run 
lA to a mean value of £ = 8.0 x 10'^'^ cm^^s^^ in Run IB. Peak instantaneous luminosities 
in Run IB reached figures of £ ~ 2.6 x 10^^ cm~^s^^. Throughout the Run 1 period, 
the Tevatron experienced several service interruptions due to causes that included power 
outages, quenching magnets, beam optics adjustments, a liquid N2 procurement delay, and 
scheduled maintenance and cost-saving shutdowns. 

The present study treated the entire Run 1 data sample. Where necessary, known run- 
dependent effects were accounted for in the analysis. With these corrections, the data set 
could be regarded as a single uniform sample. 

3.2 The Data Reduction Procedure 

From the Level 3 trigger system to the final selection at the physics analysis stage, sev- 
eral iterations of data processing and reduction occurred. The same event-driven software 
paradigm was applied at each of these stages. This uniform software control framework 
provided a flexible means of combining a virtually unlimited number of independent sub- 
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programmes, or analysis 'modules'^, into a single executable computer programme. For a 
given programme, the user could invoke, in any order, any subset of the constituent mod- 
ules. Events were handled sequentially by this subset of modules, with the output of one 
subprogramme used as input for the next one in the analysis path. A standard interface 
was used to modify the adjustable parameters for a given module. Modules also could be 
configured to function as filters, which would abort the analysis of an event in the middle 
of processing if that event failed any specified filter criteria. The possible input sources for 
this analysis framework included disk and magnetic-tape files, on-line event buffers, Monte 
Carlo generators, and user-written input modules. A single programme could support mul- 
tiple output streams, and any module could function in several analysis paths with different 
parameters used in each implementation. 

A typical event was first processed by this standard software framework on-line in the 
Level 3 trigger system (refer to Section 2.5.31 ). A single optimized computer programme 
running on a trigger farm node could implement all of the Level 3 triggers and pipe selected 
events to any of three output streams. 



Off-line, events written to 8-mm tape by the data acquisition system (refer to Section p. 61) 
entered the 'production and splitting' farms, which consisted of 64 Silicon Graphics, Inc. and 
37 IBM RS6000 computers running the UNix operating system. The production programmes 
contained modules that reconstructed higher-level physics objects from low-level detector 
output data structures, much in the same way as the Level 3 trigger did. Since run- 
dependent, as opposed to event-dependent, information was known at the off-line production 
stage, and since execution speed was less critical, the production processes could make use of 
better calibration and alignment constants, more comprehensive run-condition information, 
and more sophisticated tracking algorithms than were possible in the Level 3 trigger system. 
The production farm nodes also 'split' events into several data sets defined using physics 
analysis criteria and typically based on one or more Level 3 triggers. The rate capacity of 
the production and splitting farm computer system was ~ 1.3 x 10^ events/week. 

The third principal stage in the data reduction procedure was the application of analysis- 
specific algorithms to the reconstruction of the candidate physics processes under investiga- 
tion. In the present study, the existence of similarities in the reconstruction techniques of 



^ Analysis modules were typically coded in FORTRAN. 
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different meson candidates lent itself to a software design philosophy that employed a small 
number of highly-configurable and reusable analysis modules [117]. For example, it is clear 
that the reconstruction of the decays J/'ijj ^ A*^ /^^ and ip{2S) ^~ should use iden- 

tical coding algorithms governed by different kinematic parameters. To this end, a generic 
'dimuon finder' module was used to reconstruct J/ip and dimuon ip{2S) meson candidates. 
A separate 'two-track finder' module was used to seek decays of the form i^*(892)'^ vr^ 
by iterating through charged-track non-muon pair combinations in a given event. Finally, a 
higher-level 'charmonium-parent finder' module could be configured to seek combinations of 
charmonium candidates, often reconstructed using the 'dimuon finder', and either charged 
tracks in the event or meson objects identified using the 'two-track finder'. 

Particle types, kinematic selection criteria, and constraints used for geometric and kine- 
matic fitting could all be specified to the analysis modules at run time using the standard 
user interface. Communication between a given module and another module downstream in 
the analysis path was achieved through data structures identical to those used to represent 
the event information itself. This enabled the same module to be used more than once, but 
in different capacities, in a single analysis path. 

This modular philosophy was demonstrated by, for example, the reconstruction of the 
decay %1){2S)K^ via the '4>{2S) J / tp ir'^ it~ channel. First, the 'dimuon finder' 

module was used to find J/ip candidates in the event. Next, the 'two-track finder' module, 
configured to seek candidate dipions with the expected kinematic criteria, was invoked. 
The 'charmonium-parent finder' module then combined the J ftp and dipion candidates into 
■0(25) candidates, applying the appropriate selection criteria and fit constraints. Further 
down the analysis chain, the 'charmonium-parent finder' module was subsequently reused, 
with very different operating parameters, to form meson candidates by constructing 
combinations of the ip{2S) candidates with hitherto-unused tracks in the event. This recur- 
sive recycling of modular coding logic enhanced the internal consistency of the analysis. 



3.3 Magnetic Field Considerations 



Since the present branching-fraction study relied heavily on charged-particle trajectory re- 
construction, an understanding of the magnetic field conditions for a given event in the 
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CDF detector was essential, although somewhat less so than for mass measurement anal- 
yses. In Run lA investigations of the meson |118|] and W'^ boson |107| masses, the 



absolute momentum scale was calibrated by using the decays J/tp fj,^ , T fj.^ ^~ , 
and — > /i^ fi~ . The ensuing nominal magnetic field value used for Run lA was 1.4127 T. 
For the Run IB portion of the data sample, a nominal magnetic field value of 1.4116 T 
was used; however, this value was corrected on a run-dependent^ basis using a database of 
~1 200 magnetic field measurements performed over the course of Run IB. These corrections 



were typically <0.17% [119], and measurements of the magnetic field in the central detector 
had uncertainties of 2 x 10^^ T |107] ]. Local residual magnetic field nonuniformities, both 
in magnitude and direction, were surveyed in the CDF central tracking volume before the 
solenoid magnet was installed [ [120| ; these were corrected for in the reconstruction of track 
helices in both Runs lA and IB. 

3.4 Primary Vertex Considerations 

A knowledge of the primary vertex, or point of origin, of the decay process under examina- 
tion was important to the analysis. Primary vertex information was used in the calculation 
of primary-secondary decay lengths, momentum-pointing fit constraints, and isolation se- 
lection criteria. Two distinct aspects of the primary vertex determination will be discussed 
here: how the positions of vertices in a given event were measured, and how a single vertex 
was selected in events where more than one vertex was observed. 

The transverse and longitudinal components of the primary vertices were measured in 
two very different ways. The longitudinal (z) coordinate was established on an event-by- 



event basis using data from the VTX detector (refer to Section 2.3.2), where vertex quality 
was determined on the basis of the number of VTX hits used to identify the vertex. As 
the measurement uncertainty calculated by the vertex-finding software was deemed to be 
unreliable due to the resolution of the VTX subsystem, a fixed uncertainty of az = 0.3 cm 
was assumed for all events. The transverse (x-y) coordinates of the primary vertices in 
an event were typically calculated using the measured run-averaged beam position. The 
rationale for this was that the transverse beam position typically varied less than 10 nm 



'run', in this context, refers to a period of uninterrupted data collection, typically lasting several 

hours. 
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in either the x or y directions over the course of a single data-taking run and that any 
event-by-event transverse coordinate measurements would be biased by fluctuating track 
multiplicities and event topologies in individual events. The slopes and intercepts of the run- 
averaged beam position were therefore combined with the event-by-event z locations of the 
vertices in an event to determine the transverse positions of those vertices. The transverse 
coordinate measurement uncertainties were fixed to {crx,cry) = (25,25) ^m, corresponding 
to the observed circular beam spot size in the transverse plane. It should be noted that, in 
the unusual cases where no run-dependent beam position information or no VTX z-vertex 
information was available for a given event, the lacking vertex coordinates were computed 
using the available SVX and CTC track information for that event. 

During the latter stages of Run IB, when instantaneous luminosities often exceeded 
those in Run lA by an order of magnitude, the primary vertex multiplicities also increased 
dramatically. Whereas in Run lA the average number of high-quality vertices in a given 
event was ~1.6 with ~3% of events having at least four such vertices. Run IB events 
averaged ~2.9 high-quality vertices with ~5% of events having at least eight such vertices. 
Once a pair of muon candidates had been identified using the 'dimuon finder' machinery 
described in Section |3.2| (which was a procedure that did not directly employ primary 
vertex information), the z coordinates of the two candidates were used to select a single 
vertex. Specifically, of those vertices possessing the highest quality classification for the 
given event, the vertex that had the shortest longitudinal displacement from either of the 
two muon candidates was chosen as the primary vertex corresponding to the muon objects. 
The position coordinates of muon candidates, as opposed to those of other tracks used in 
this analysis, were used to select primary vertices because the two muon candidate tracks 
constituted an unambiguous part of the final state under study. 



3.5 Track Quality Criteria 

The imposition of quality requirements on the tracks used in the analysis was intended 
to reduce those backgrounds arising from poor track measurements in the CDF detector. 
Track candidate fits reconstructed for this study were required to have used at least four 
hits in each of at least two axial CTC superlayers (refer to Section p. 3. 3 for a description 
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of the CTC). These track fits also had to use at least two hits in each of at least two stereo 
CTC superlayers. No requirement was made on which two of the five axial and four stereo 
superlayers were to be used in the fit. 

In this study, information from the VTX and CTC subsystems was employed in the 
reconstruction of track paths. For all such tracks, the helical trajectories were extrapolated 
back into the SVX where associated hits were sought using a road algorithm. If a sufficient 
number of good SVX hits was found, then the track was refit using all of the relevant VTX, 
CTC, and SVX information and the resulting track helix was used. Performance disparities 
in the two silicon microstrip vertex detectors used in Runs lA and IB (refer to Section [2.3.1 ) 
motivated two different associated SVX hit requirements: >3 hits in the SVX, and >2 hits 
in the SVX'. Moreover, track fits that made use of SVX information were only considered 
by this analysis if the SVX x^/hit, defined as the increase in the track fit per SVX hit 
due to the inclusion of SVX hits in the CTC track fit, satisfied the condition x^/hit < 6.0. 
This requirement was similar to one used in CDF studies of the top quark ||3^ . 

Tracks possessing transverse momenta pT ~ 250 MeV/c were not reconstructed in the 
production stage of the data reduction process (refer to Section . Useful measurements 
of tracks with px < 250 MeV/c would have been difficult due to the number of track helices 
that subtended >360° in azimuth while inside the CTC ('loopers') and due to the dearth of 
available hits populating the outer superlayers of the CTC. A study of the low-px pattern 
recognition efficiency (see Appendix j^) indicated that a requirement of > 400 MeV/c 
would ensure that candidate tracks would be in a pT region with a relatively constant and 
measurable tracking efficiency. This pT requirement was imposed for all tracks prior to any 
corrections for multiple-scattering and energy- loss effects. 

Similarly, tracks with trajectories at high absolute pseudorapidities tended to deposit 
fewer hits in the CTC, and therefore were reconstructed less efficiently. Requirements on 
the radii at which extrapolations of track helices intersected one of the endplate planes 
of the CTC, "rcTC) were used to remedy this problem by diminishing the pseudorapidity 
acceptance. For example, a cut of t^-^q > 132.0 cm on this CTC exit radius would have 
required that the track in question traverse all nine superlayers of the CTC. The tracking 
efficiency study outlined in Appendix j^, however, showed that a requirement of r^j^c ^ 
110.0 cm, corresponding to the radial position of the outer edge of the second-outermost 
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axial superlayer, defined a set of tracks with high efficiency while minimizing the reduction 
in geometric acceptance. This requirement was imposed for all non-muon track candidates 
after multiple-scattering and energy corrections had been implemented. Muon candidates 
were not subjected to an r^^*^ cut in order to accept candidates identified by the CMX 
subsystems, which represented ~22% of the total central muon fiducial acceptance. Since 
the muon candidates in all the channels had comparable r] distributions, any muon tracking 
inefficiencies resulting from this lack of a CTC exit radius criterion divided to unity in the 
final calculations of the ratios of branching fractions. 

3.6 Optimizing the Kinematic Selection Criteria 

The kinematic selection criteria were chosen to achieve accurate and unbiased measurements 
of the i?-candidate event yields while maximizing both the rejection of the background and 
the statistical significance of the signal. Wherever possible, uniform criteria were used to 
reconstruct the different decay modes. This served both to maximize internal consistencies 
in the study and to reduce the magnitude and number of systematic uncertainties. 

In order to minimize any bias caused by fluctuations in the background levels in either 
the signal or sideband regions of the invariant mass distributions, all of the i?-meson decay 

modes were used in the determination of the kinematic selection criteria. The figure of merit 

Ns 

that was maximized was the signal significance, S = — The N^ symbol denotes 

^/Ns + Nb 

the number of signal events calculated with a Monte Carlo procedure (refer to Chapter |^ 
for a given decay mode and a given time-integrated luminosity. The expected number 
of background events, N^, was estimated from the data by extrapolating the background 
rate in the observed i?-meson invariant mass sidebands to the background yields under the 
signal regions, defined to lie within three standard deviations of the observed means of the 
S-meson resonances. The resultant optimal kinematic selection requirements, which are 
detailed in the remaining sections of this chapter, varied little within each of the charged 
and neutral -B-decay mode categories. 

The optimum kinematic cuts determined by maximizing S for each of the decay channels 
were confirmed by maximizing a different measure of statistical significance derived entirely 




was taken to be the event-yield uncertainty returned 



3.7. MUON CANDIDATE SELECTION 



59 



by the fit of the given invariant mass distribution to a single fixed-width Gaussian signal 
and a linear background parameterization. Although this latter method was more sensitive 
to statistical fluctuations, its conclusions were consistent with those of the S'-maximization 
technique. 



3.7 Muon Candidate Selection 

The minimum transverse momentum required of a muon to traverse the central calorimeters 
and the solenoid magnet at r/ ~ and reach the CMU subsystem was ~1.4 GeV/c. A raw 
kinematic requirement of pT > 1-4 GeV/c was therefore placed on all muon candidates in 
the analysis. Because the numerators and denominators in the ratio-of-branching- fraction 
calculations involved muons with identical selection criteria, there was no need to demand 
that all candidate muons populate the plateau region of the trigger efficiency parameteriza- 
tions (refer to Section 4.5.2| ). As a result, the fact that the lower edge of this requirement 



fell in inefficient regions of the measured Level 1 and Level 2 px-dependent trigger efficiency 
distributions was of little consequence. A higher requirement would have unnecessarily 
weakened the statistical significance of those i?-meson decays with a ip{2S) candidate in 
the final state. 

A major source of muon background was that due to charged kaons and pions decaying 
to yield muons within the CDF tracking volume. In some cases, the charged kaon or pion 
tracks were reconstructed in the CTC, and the daughter muons registered stubs in the muon 
systems. A second major source of muon background was that due to hadronic 'punch- 
through' particles (see Section 2.4.2), namely hadrons that passed through the calorimeters 



and entered the muon systems. 

In order to reduce these backgrounds, a track candidate in the CTC, when extrapolated 
out to the muon chambers, was required to match the position of a muon stub. This condi- 
tion was only met if the muon stub and the extrapolated CTC track in question matched 
within three standard deviations of the multiple-scattering and measurement uncertainties 
in both the transverse (r-99) and longitudinal (2) planes. In cases where a CTC track could 
be extrapolated to stubs in more than one muon subsystem. Boolean "or" operations were 
implemented to combine the appropriate matching requirements. In Run lA, the muon 



60 CHAPTER 3. SELECTION OF CANDIDATE EVENTS 



matching criteria for J/i/; dimuons had an efficiency of (98.66 it 0.04)% and increased the 
J/ip signal-to-background ratio from 3.08 to 3.61 in the 3.0 — 3.2 GeV/c^ invariant mass 
interval |121|; comparable effects were expected for Run IB. 



3.8 Charmonium Reconstruction 

The three charmonium decay modes used in this study were J/V' /U^, ^(25) fi^ ^~ , 

and '4>{2S) J/tpn^ it~, where in the latter case the J ftp meson was reconstructed in its 
dimuon mode. These channels were used over other charmonium decay modes because the 
fi^ final state could easily be identified using the CDF trigger system while also making 
possible the rejection of several background processes. This section describes the selection 
of charmonium candidates; the resulting inclusive charmonium invariant mass distributions 



are presented in Section 5.1, where they are used in the determination of certain efficiency 
corrections. 

3.8.1 Dimuon Charmonium Decays 

Dimuon charmonium candidates were formed using the 'dimuon finder' module (see Sec- 



tion 3.2) by considering all the muon candidates in a given event that met the criteria 



outlined in Section 3.7. The two candidates constituting a muon pair were required to 
possess charges of opposite sign. 

Background contributions were reduced significantly by performing a least-squares fit 
of the two muon track candidates and applying a 'vertex constraint', which forced the two 
tracks to originate from a common point in space |122| ]. In this fit, an initial approximation 
to the track parameters was found, corrected for multiple scattering and dE/dx (energy- 
loss) effects, and then adjusted under the vertex constraint so as to minimize the y^. The 
confidence level, CL{x^)i of the fit was required to exceed 0.01. 

When J/tp and ip{2S) candidates were used in the reconstruction of an exclusive B 
final state, an additional fit was performed on the dimuon system to improve the invariant- 
mass resolution of the B candidate. In this case, the fit was done with the simultaneous 
application of a vertex constraint and a 'mass constraint' [|122(| , which required that the 
dimuon mass equal the appropriate world average J/tp or ip{2S) mass of 3.09688 GeV/c^ 
or 3.68600 GeV/c^, respectively [pO|. The confidence level of each of these vertex-plus-mass 
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constrained fits was also required to satisfy the condition CL{x^) > 0.01. 
3.8.2 Hadronic Cascade ^^(25) Decays 

Candidates for the decay ip{2S) — > J/tp'jr~^TT~ were chosen by combining J/tp dimuons, 



which were selected as described in Section 3.8.1, with dipions identified using the 'two- 



track finder' module described in Section 3.2. The two pion candidates were required to 



meet the track quality criteria detailed in Section 3.5, to have charges of opposite sign, 
and to have an invariant mass, prior to any vertex or vertex-plus-mass constrained fits, 
in the range 0.35 < M{tt^tt^) < 0.61 GeV/c^. The lower limit of this mass range was 
motivated by the expected dipion invariant mass distribution for ip{2S) decays, discussed 



in Section |4.2.1| . The upper limit corresponded to the maximum kinematically-allowed 
dipion invariant mass, defined by the ip{2S) — J mass difference with an allowance for 
measurement uncertainty. 

In order to reduce background effects, a vertex-constrained fit was performed on the 
four-track dimuon-dipion system, with the two candidate muon tracks simultaneously con- 
strained to form an invariant mass equal to the world average J/V' mass [^0| . The confidence 
level of this fit was required to exceed 0.01. As in the case of the dimuon charmonia de- 



scribed in Section 3.8. Ij , when an hadronic cascade decay of the ip{2S) meson was used as 



a daughter process in the reconstruction of another exclusive decay, the four-track system 
was simultaneously subjected to an additional mass constraint, namely the world average 
tp{2S) mass [^0|. Again, the probability of this fit was required to satisfy the criterion 
CLix^) > 0.01. 

3.9 5-Meson Candidate Reconstruction 

The reconstruction of the decay modes B+ iij{2S) K+ , V(25') ivT* (892)°, S+ 

J/tpK'^, and B^ — > J /ipK* (892)^ was achieved by combining one of the three sought- 
for charmonium meson candidates in a given event with either K'^ or K*{892)^ candidate 
mesons, as appropriate. In spite of the facts that six different decay chains were recon- 
structed and that the final-state charged-particle track multiplicities forming these chains 
ranged from three to six tracks, the methods used to identify candidate S-meson decays 
were kept as uniform as possible. This was due in large part to the modularity provided 
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by the data reduction procedure described in Section 3.2. The fohowing sections describe 
the selection of kaon candidates, the selection of i?-meson candidates, and the treatment of 
multiple invariant mass hypotheses. 

3.9.1 Kaon Candidate Selection Criteria 

The CDF detector lacked the ability to differentiate kaon candidates from pion candi- 
dates, a fact that necessitated the consideration of all eligible tracks and the assignment 
of an appropriate kaon or pion mass hypothesis to these tracks. A consequence of these 
'blind' hypothetical mass assignments was the amassment of significant combinatorial back- 
grounds. Fortunately, the charged particle inclusive cross section in pp collisions is a rapidly 



falling function of transverse momentum |123], thus enabling the use of reasonably efficient 



minimum-pT criteria to remove a significant fraction of the combinatorial background. 

The pt > 400 MeV/c cut (see Section |3.5D , imposed to avoid low-px tracking ineffi- 
ciencies (see Appendix also served to reduce combinatorial backgrounds caused by the 



lack of particle identification. The optimization scheme described in Section |3.6| , however, 
supported the application of higher cuts to combat these combinatorial backgrounds. 
The criteria applied were therefore pt{K~^) > 1.5 GeV/c for kaon candidates used in the 
reconstruction of candidates and pT^{K*{892)^) > 2.0 GeV/c for the two-track K-n 
candidates used in the reconstruction of candidates. Note that the two tracks consti- 
tuting each K*{892)^ candidate did not have individual pT cuts imposed on them beyond 
the universal 400 MeV/c track quality criterion. 

In a given event, K*{892f candidates were reconstructed by considering all track pairs 
with charges of opposite sign. For each pair of tracks considered, both mass assignment 
hypotheses were initially retained. To reduce combinatorial backgrounds further, the invari- 
ant mass of each K-tt candidate, M{K^Tr^), was required to fall within a window centred 
at the world average K*{892)^ mass, 896.10 MeV/c^ The size of this invariant mass 
window was 160 MeV/c^, a requirement that was estimated to be ~80.5% efficient under 



the assumptions that the i^T* (892)" resonance could be described by a Breit-Wigner [124| 



line shape with an intrinsic width of F = 50.5 MeV |3C] and that the experimental resolution 
was significantly less than F. Finally, no constrained fits were used to select and reconstruct 
-R'*(892)*^ candidates explicitly, although the two constituent tracks did undergo some ad- 
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justments in subsequent global constrained fits of i?-meson candidates (see Section 3.9.2). 
The invariant mass window criterion described above was imposed on the K-n candidates 
after they had undergone global constrained fits. 

3.9.2 5-Meson Candidate Selection Criteria 

B meson candidates were reconstructed by forming combinations of charmonium candidates 
(J/V', i^C^S) dimuon, or ^{2S) hadronic cascade modes) with either or K*(892)^ candi- 
dates. In a manner similar to that used in the reconstruction of charmonia, a least-squares 
kinematic fit was performed on the B candidate daughter tracks under the constraints that 
all the tracks originate from a single common secondary vertex and that the charmonium 
candidates possess an invariant mass equal to the applicable world average mass |Q . It was 
possible to constrain all the B decay daughters to a single decay point because the distances 
traveled by charmonium and K*{892)^ mesons before decaying were negligible compared to 
the decay vertex resolution of the CDF detector |122| ] . 

In the case of i?-meson reconstruction, a 'momentum-pointing' constraint was used in 
the kinematic fits. This additional constraint required the flight path direction of the B 
candidate to be parallel to its momentum in the transverse (r-ip) plane. The pointing 
constraint was not performed in three dimensions due to the large uncertainty on the z 



component of the primary vertex position (refer to Section 3^). The confidence level of the 
global least-squares fit of each B meson candidate, with vertex-plus-mass and momentum- 
pointing constraints applied, was required to exceed 0.01. 



The kinematic selection optimization, discussed in Section 3^, resulted in a > 
6.0 GeV/c transverse momentum requirement for B^ candidates and a pT > 9.0 GeV/c 
requirement for B^ candidates. 

Motivated by the understanding that b quarks were expected to 'fragment' in a way that 
imparted most of their momentum to the ensuing B meson (References ||6^, ^] and Sec- 
tion p..2.1| ), an isolation criterion was imposed on the B meson candidates. This requirement 
was expressed using the variable 

R 

J2 Pi-PB 

Ib = (3.1) 



64 CHAPTER 3. SELECTION OF CANDIDATE EVENTS 



where pB was the 3- momentum of the i?-meson candidate and the pi were the momenta 
of additional particles, other than those constituting the B candidate, contained within a 



cone of radius R = y (A^?)^ + (Ar/)^ < 1.0 and with its axis collinear with the S-candidate 
3-momentum direction. In addition, non-B tracks were only included in the sum in Equa- 



tion 3.1 if their longitudinal displacement parameter, zq, lay within 5 cm of the primary 
vertex location corresponding to the B candidate in question. This improved the efficiency 
of the isolation criterion by not discarding B candidates that were unisolated due to tracks 



from other primary vertices in the event. Figure |3.1| depicts the distribution of Ib for 
candidate B^ — > J/ipK^ decays, after the combinatorial background under the B^ de- 
cay signal was statistically subtracted using events in the J /ip invariant mass sideband 
regions. The criterion imposed was Ib < 7/13, which resulted from the optimization proce- 



dure described in Section 3.6. Figure |3.1| illustrates that the Ib requirement was relatively 
efficient. 

The long lifetimes of and 5° mesons, (1.64 ± 0.05) ps and (1.55 ± 0.05) ps |o|, 
respectively, made possible the rejection of background with a cut on the proper decay 
length, ctb, defined by 

CTB = 2 — "^B, (3.2) 

Pt 

where xt was the distance between the primary and secondary (S-candidate) vertex, pro- 
jected onto the transverse plane. The quantity ttib represented the invariant mass of the can- 
didate B meson. In the present analysis, a proper decay length requirement of ctb > 100 fim 
was imposed on S-meson candidates. 

Figure ^]2| depicts a sample event display of a candidate i?-meson decay that was identi- 
fied as such by the data reduction algorithms and selection criteria described in this chapter. 

3.9.3 Multiple Invariant Mass Combinations 

Multiple S-meson candidates that passed all the selection criteria for a given event were 
a significant source of background. Such additional combinations usually resulted from 
other tracks that formed a vertex with the two muons in the event and were therefore 
more prevalent in reconstructions with higher daughter-track multiplicities, such as those 
involving the ip{2S) J/ip vr^ vr^ mode. Another obvious source of multiple combinations 
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Figure 3.1: The distribution of the isolation variable Ib for candidate — > J/tp decays. 
A background subtraction was performed using the sidebands in the J/ ^ invariant mass 
distribution. The arrow marks the point below which candidates were accepted. 
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Figure 3.2: A sample event-display diagram of an observed candidate decay. The 
view in this figure is similar to that in Figure 2^. CTC sense wire hits are visible for 
all nine superlayers, and the transverse projections of the reconstructed candidate track 
trajectories are indicated by smooth curves. The inner circle denotes the boundary between 
the CTC and VTX, and the bars outside the CTC represent energy deposition in the 15° 
central calorimeter wedges. Near the top of the figure are listed the missing transverse 
energy {^) and its azimuthal direction, the total transverse energy for the event, and the 
maximum energy recorded in a single calorimeter wedge. The "+" symbols represent stubs 
in the central muon chambers. The data reduction algorithms and event selection criteria 
identified this event as a candidate for the decay — > ip{2S) , where the ip{2S) meson 
was reconstructed in its ip{2S) — > J/tpTT~^7r~ mode and the J/V' meson was reconstructed 
in its J /^l) ^ n'^ ii~ mode. The identified track candidates are indicated in the figure. The 
azimuth and pseudorapidity coordinates in the bottom right corner are those of the track 
with the highest pT, the candidate, a close-up view of which is illustrated in the left 
panel of the figure. 
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was in the case of i^*(892)'^ candidates, where the K-n mass assignment ambiguity and the 
size of the K*{892)^ mass window cut required by the natural width of the i^'*(892)'' meson 
increased the probabihty of additional invariant mass combinations in the event. 

In the case of the i^*(892)*' candidates, three treatments of the K-tt mass hypothesis 
ambiguity were considered: (a) random selection of one of the mass hypotheses, (b) selec- 
tion of the hypothesis with M{K^-k~) nearest the accepted world average K*{892)^ pole 



value |3C], and (c) parameterization of the relative contributions and resolutions of the right 
and wrong mass assignment hypotheses in the signal region. All of these techniques would 
have required Monte Carlo calculations. 

In this analysis, the last of the above three treatments was adopted in conjunction with a 
general procedure that made use of the probability of the global constrained least-squares 
fit, CL{x^), to remove multiple combinations. The prescription used was the following: 

1. All selection criteria were applied, including the M(K^Tr^) requirement. 

2. For a given event, only the candidate possessing the highest global CL{x^) value 
returned by the constrained fit was retained. 

3. For events where both i^*(892)° K -TT mass hypotheses survived all of the selection 
criteria, the two CL{x^) values were often similar. If one of the M{K^tt~) combina- 
tions satisfied the highest-CL(x^) criterion, and a second combination that differed 
only in the K-ir mass assignment existed, then the second entry was also retained. 

4. The signal observed in the B candidate invariant mass distribution of those modes 
containing a i^*(892)'' candidate was fit using two Gaussian distributions with their 
relative widths and amplitudes constrained by ratios determined using the Monte 
Carlo calculations described in Chapter ^. The constraints are listed in Table 3.1 



Refer to Section S.l for a synopsis of the 5-meson candidate event yields. 
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B Modo 


cc Modo 




(TUT 


S° ^ J/il;K*{892f 




0.0683 


3.623 


^iP{2S) (892)0 


V'(25) ^ /x+ /X- 


0.0456 


5.750 


5° ^-0(25) ii:*(892)° 




0.0429 


6.327 



Table 3.1: The relative amplitudes (^) and widths {a) of the right {RC) and wrong (WC) 
K-TT combinations, as determined by fits to Monte Carlo invariant mass distributions using 
double-Gaussian signal parameterizations. 



Chapter 4 

Geometric and Kinematic 
Acceptance 



In branching fraction analyses, a knowledge of the proportion of candidates that go unob- 
served due to both the fiducial detector geometry and the kinematic selection criteria is 
essential. Since the data themselves could not be used to determine these acceptances and 



efficiencies, it was necessary to make use of Monte Carlo |125| calculations. 

This chapter describes the Monte Carlo machinery used to generate bottom quarks, 
hadronize these quarks into B mesons, decay the B mesons into final-state particles, simulate 
the signatures left by the particles in the CDF detector, and simulate the first two levels of 



the triggers relevant to the study. Neither the underlying event (refer to Equation 1.4) nor 
multiple pp interactions were modeled in these Monte Carlo studies. 

4.1 The Monte Carlo Generation of B Mesons 

Single b quarks were generated according to an inclusive transverse momentum {kT{b)) 



spectrum based on a next-to- leading order QCD calculation |54] that used the Martin- 



Roberts-Stirling MRSD0 parton distribution functions |l2q], a renormalization scale of 



H = fiQ = J ml + A;^(6), and a 6-quark mass of m;, = 4.75 GeV/c^. Refer to Section 1.1.3 



for an account of 6-quark hadroproduction. The b quarks were produced in the rapidity 
range —1.1 < yt < 1-1 with > 5.0 GeV/c and fragmented into B mesons according 



to a model that used the Peterson fragmentation function [S5| with the Peterson param- 
eter defined to be 0.006 [^] (refer to Section 1.2.1| ). Flavour was conserved in both the 
production and fragmentation of b quarks. 
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4.2 The Monte Carlo Decay of B Mesons 

Decays of Monte Carlo generated B mesons into charmonium and kaon final states were per- 



formed using a modified version of the CLEO QQ Monte Carlo programme [127]. Properties 



of the relevant particles, including mass, lifetime, and intrinsic width, were updated in the 



programme to reflect the current world- average values [30|. For the pseudoscalar — > vec- 
tor + vector decays, the decay helicities were nominally set to the central value of the world- 
average longitudinal polarization fraction measured for the decay J/tp K*{892)^ , 



Tl/T = 0.78 ± 0.07 [0, H, |T28[. Finally, a customized QQ Monte Carlo matrix element 
was constructed to model correctly the observed kinematics of the pions in the charmonium 
decay ip{2S) J / tp ir^ . The following section describes this customization. 

4.2.1 The ^(25) J/V'tt+tt^ Monte Carlo Matrix Element 

The dipion kinematics in the decay ip{2S) J/ipii^T:~ were studied to search for tech- 
niques to remove background efficiently and to ensure that the pion tracks were simulated 
correctly. Whereas the generated QQ transverse momentum distribution of the dipion sys- 
tem was consistent with that in the data, the simulated QQ dipion invariant mass was not 
in agreement with the data. 

The default QQ Monte Carlo programme determined the kinematics of the J/'^vr+vr" 
decay products using pure three-body phase space; however, at the time of the discovery 
of the ip{2S) J/ipir'^ ir" decay channel, it was observed that the angular distributions of 
the pions were isotropic |32|. In subsequent phenomenological investigations of the decay 
amplitude, the absence of any observed angular correlations made apparent the fact that the 



amplitude had a strong dependence on the invariant mass of the vr^ vr^ pair ||129| , |130| , 131 

A set of customized routines was therefore written to model the ip{2S) J/ip'K^'K~ 
matrix element in the QQ Monte Carlo programme using the phenomenological prescrip- 
tion. For the purposes of the algorithm, the decay was considered to have the logical form 
ip{2S) — > J/tpV, where V was taken to be a dipion pseudostate {V vr+vr^). In the first 
step, the invariant mass of V was calculated using the V mass distribution |131|], 



which is in a form similar to that used in Reference [132]. The mass of the V pseudostate 
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Figure 4.1: The dipion invariant mass distribution calculated using a phenomenologically- 
prescribed QQ Monte Carlo matrix element (solid histogram). The dashed histogram illus- 
trates the results of a pure phase-space calculation. 



was extracted using the Von Neumann acceptance-rejection Monte Carlo method |3C]. The 
second and third steps in the algorithm consisted of decaying, using phase space alone, 
the modes '0(25') — > J/^jV and V vr+vr", respectively. The J/il) — > fJ-'^ fJ-^ decay was 
performed using the default QQ matrix element for that decay. 

Figure illustrates the phenomenologically-inspired and phase-space QQ dipion invari- 
ant mass distributions, as determined from Monte Carlo calculations. The correspondence 
between a parameterization of the phenomenological prediction and the data is pictured 
in Figure where a sideband subtraction was used to remove the background to the 
ip{2S) candidate decays. A recent precision measurement of the ip{2S) J/ipn'^'K" 
branching fraction reports a similar M(7r"^7r~) distribution |133(| . A conservative cut of 
M(7r+7r-) > 0.35 GeV/c^ (refer to Section ^^2]) was applied to the dipion candidates prior 



to the application of any vertex, vertex-plus-mass, or momentum-pointing constrained fits. 
The M(7r'^7r~) requirement was found to be efficient and capable of rejecting a significant 
number of background candidates. 
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Figure 4.2: The observed background-subtracted distribution of the dipion invariant mass 
(points) in decays of the form ip{2S) — > J / ip it~^ tt~ . The arrow indicates where a cut at 
0.35 GeV/c^ was subsequently imposed. The sohd curve represents a phenomenological 



prediction due to Pham et al. |13l|] ; the broken curve describes a pure phase-space param- 
eterization. 
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4.3 Simulation of the Detector Response 

Once Monte Carlo B mesons were generated and decayed into their daughter states, a 
fuh Monte Carlo simulation of the CDF detector response was invoked. The simulation 
produced raw data structures that were in most cases identical to those read out from 
the DAQ system using collider data (refer to Section |2.6| ); this enabled the subsequent 
processing of Monte Carlo events in a manner almost identical to that used for the data. 
The corrections for local nonuniformities in the magnetic field, described in Section 



however, were switched off at the detector simulation and production (see Section 3.2) 
processing stages. 



As outlined in Section 3.1, the CDF detector geometry changed between Runs lA and 
IB. For the purposes of the detector simulation in the present study, the Run IB geometry 
information was used to represent the CDF detector for the entire Run 1 period. Within 
the Monte Carlo statistics used for this analysis, this assumption appeared to be reasonable 
for the calculation of relative geometric acceptances and kinematic selection efficiencies. 

Greater care was necessary in the characterization of the beam profiles near the primary 
vertex. The simulated transverse beam profile was tuned to match approximately that 
observed in the data. The longitudinal beam profile, however, was more important to 
this study due to its implications for the fiducial SVX acceptance. The z location of the 
simulated primary vertex was forced to be distributed according to the sum of two Gaussian 
probability density functions with means and standard deviations of 2.3 cm and 2.1 cm, 
and 18.4 cm and 35.6 cm, respectively. The relative normalization of these two Gaussian 
distributions was approximately 2:3. 

The Monte Carlo events were processed by the same production machinery (Section |3.2D 
used on collider data. This ensured that the data and the Monte Carlo events were both 
subjected to any biases, should they exist, inherent to the reconstruction algorithms. 

4.4 Simulation of the Level 1 and 2 Triggers 



After the CDF detector response to the Monte Carlo decay fragments had been simulated, 
the events were passed through a Level 1 and Level 2 dimuon trigger simulation module. 
(Refer to Sections 2.5.1 and 2.5.2 for descriptions of the Level 1 and Level 2 triggers, respec- 
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tively.) This module used parameterizations of the measured muon trigger efficiencies to 
determine the probabihty that a given candidate dimuon event satisfied the trigger require- 
ment. The Monte Carlo events were then selected using a Monte Carlo procedure based on 
this probability. 

Figure |4.3| illustrates the Level 1 CMU and CMX trigger efficiency parameterizations as 
functions of muon transverse momentum. These efficiencies were assumed to be functions 
of pt alone and to remain unchanged during the course of Run 1. In the case of the 
CMX efficiency, the veto scintillation counters (refer to Section 2.4.3| ) degraded the plateau 
efficiency uniformly. 

The Level 2 dimuon trigger efficiencies were found to depend not only on pi, but also 
on charge, pseudorapidity (??), azimuth ((^), and time- integrated luminosity (/ £ dt). Due to 
the geometry of the CTC, positive muon tracks were accepted by the trigger more efficiently 
than negative tracks. The rj dependence of the efficiency was observed to be parabolic due 
to the increased hit efficiency of higher-r/ tracks]^. The ip dependence of the efficiency was 
observed to be sinusoidal due to the offset of the beam from the geometrical centre of 
the CTCQ Finally, the decline in the CTC hit efficiencies as a function of JC dt caused a 
corresponding degradation in the Level 2 CFT efficiency. This effect also caused an increase 
in the curvature of the r/-dependent efficiency parabola. 



Figure 4.4 depicts the Level 2 CMU and CMX trigger efficiency parameterizations as 
functions of muon transverse momentum. The curves shown are based on 8.89 pb~^ of 
positively charged muon candidates that fell in a particular IB run range (65 001 — 66 000) 
and the lowest px bin (~2.2 GeV/c) of the CFT. 

Note that the Level 3 dimuon trigger efficiency, which was measured to be 0.97ib0.02 and 
independent of pT, was not explicitly simulated in the geometric and kinematic acceptance 
calculations. 



^Muon tracks with large |r7| values deposited more charge on the CTC wires; this increased the pulse 
height, which in turn e nhanc ed the hit efficiency. 

^The CFT (Section 2.5.2) did not measure the impact parameter, thereby introducing a false curvature 
term, which varied with ip. 
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Figure 4.3: The Level 1 1ow-]3t CMU and CMX trigger efficiency parameterizations. The 
dashed-dotted curves are the central values and the solid curves represent shifts of the 
plateau efficiencies and pt values by one standard deviation. 
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Figure 4.4: The Level 2 low-px CMU and CMX trigger efficiency parameterizations for 
positively charged muons in the IB run range 65 001 — 66 000 and the lowest-pT CFT bin. 
The dashed-dotted curves are the central vahies and the solid curves represent shifts of the 
plateau efficiencies and pT values by one standard deviation. 
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B Mode 


cc Mode 


Acceptance 


B+ J/i) K+ 
B+ 4^{2S) K+ 
B+ i){2S) 


Tp{2S) J/V'TT+vr- 


0.01920 ±0.00025 
0.02183 ± 0.00027 
0.00631 ±0.00015 


B^ J/V^ivT* (892)0 
B^ ^i){2S) K*{S92f 
B^ ^i){2S) K*{S92f 


ip{2S) 11+ fi- 
ip{2S) J/V'TT+vr- 


0.00777 ± 0.00016 
0.00828 ± 0.00017 
0.00385 ±0.00011 



Table 4.1: A summary of the absolute products of the geometric and kinematic acceptances, 
calculated for each decay mode for B mesons produced with /ct(^) > 5.0 GeV/c and —1.1 < 
Uh < 1.1. The grouping is on the basis of common selection criteria, and the uncertainties 
given are due to Monte Carlo statistics alone. 



4.5 The Acceptance Calculations 



In the final step of the default Monte Carlo procedure, the events were processed by the 
same analysis-specific user algorithms used on the collider data. The reconstructed Monte 
Carlo events were used to estimate the geometric and kinematic acceptances for the decays 



under study. The products of these quantities are listed in Table 4.1. The Monte Carlo 



statistical uncertainties on the absolute acceptance values listed in Table 4.1 ranged from 
1.2% to 2.9%. 



Table 4.2 lists the appropriate relative geometric and kinematic acceptances necessary 
for the branching-fraction ratio calculations. The syntax in Table |4^ follows that introduced 
in Table |1.3| , but with particular entries for the two different iIj{2S) modes. The trends in 



Table 4.2 are consistent with differences in daughter-track multiplicity and cc meson mass. 
When the relative acceptance values were used in the branching-fraction calculations, a 
systematic uncertainty of 3.0% was applied to account for differences between two differ- 



ent CDF detector simulations and between the Run lA and IB detector geometries [134|. 
This systematic uncertainty in the acceptance-measurement method precluded a need for 
increased Monte Carlo statistics. 
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Relative 
Acceptance 


B+^J/il)K+ B^ ^ J/i)K*{d,mf B+^il){2S)K+ 
J/il) — > fJ-+ ^J-^ J/V' /^^ A*^ V'(2'S') — > tJ''^ fJ-'^ 


rO T/ih K*(RQ'?')^ 
J/V- 11^ M" 


0.405 ±0.010 

1.137 ±0.020 2.810 ±0.068 
0.329 ± 0.009 0.813 ± 0.025 0.289 ± 0.008 
0.431 ±0.010 1.066 ±0.031 0.379 ± 0.009 
0.201 ±0.006 0.496 ±0.018 0.177 ±0.006 


i^{2S) ^ 11+ fi- 


B+ V(25) K+ 


^{2S) K*{892y> 
ij{2S) fi+ 


B^ iIj{2S) K*{m2f 



Table 4.2: A summary of the relative products of the geometric and kinematic acceptances, 
calculated for each ratio of decay modes for B mesons produced with /ct(^) > 5.0 GeV/c 
and —1.1 < yh < 1.1. The indicated uncertainties represent the Monte Carlo statistics 
alone. 

4.5.1 Differential Production Cross Section of Generated b Quarks 

A significant source of systematic uncertainty on the relative acceptance calculations re- 



ported in Section |4.5| was that due to the assumed generated 5-quark momentum spectrum. 
Although the measured S-meson differential cross section has been observed to have a 
shape in agreement with theoretical predictions, the absolute rate was near the limits of 
that predicted by typical variations in the theoretical parameters [p^]. An estimate of the 
acceptance uncertainty due to the calculated 6-quark momentum spectrum was taken by 



varying the default Monte Carlo generation parameters detailed in Section 4.1 from their 
nominal values. The &-quark mass was varied from its central value of mb = 4.75 GeV/c^ 
to 4.50 GeV/c^ and 5.00 GeV/c^, and the renormalization scale was varied from a nominal 
value of /i = /iQ to /io/4 and 2fiQ. The effects of these variations on the acceptance ratios were 
examined, and one half of each maximum deviation was taken as the appropriate systematic 



uncertainty on the relative acceptance. Table 4.3 summarizes the systematic uncertainties 



attributed to the generated differential production cross section of the b quarks. 
4.5.2 Trigger Effects 



Even though the numerators and denominators of the branching-fraction ratios in Table 1.3 
involve muon pairs that originate from common data samples and the pt spectra of the B+ 
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6-Quark Production 
Spectrum Systematics (%) 


B+ J/'4)K+ 
J/ip 11+ 11- 


B" J /'ip K* (892)^ 
J/i) 11+11- 


B+ tp{2S) K+ 
^Ij{2S) 11+ ^i- 


1 111, K*(fi.QO\^ 


4.0 














2.4 


5.9 




tp{2S) J/ijTT+TT- 


2.8 


5.0 


5.5 


5^ ^ ip{2S) K*{892f 
ij{2S) n+ 11- 


5.5 


1.9 


5.8 


'il){2S) K*{892f 
^(25) ^ J/^7r+ TT- 


6.7 


3.4 


9.6 



Table 4.3: A summary of the systematic uncertainties (expressed in %) on the relative 
geometric and kinematic acceptance due to variations in the differential production cross 
section of the generated b quarks. 

and B° mesons are expected to be comparable, some care is required due to the possibility 
of topological effects when considering relative trigger acceptances in branching-fraction 
ratio studies. In approximately 75% of the selected events, the two muon candidates that 
were identified as charmonium daughters were also the muon candidates identified by the 
dimuon trigger system. In most of the remaining events, an additional muon candidate in 
the event satisfied the dimuon trigger requirements. These 'volunteers' were included in the 
analysis in order to maximize the sensitivity of the data sample. 

In order to determine systematic uncertainties on the relative acceptance values that 
accounted for candidates that were triggered in ways not treated by the dimuon trigger sim- 
ulation, the relative topology dependence of the unsimulated triggers was investigated by 
examining some representative yield ratios. Three different topological scenarios with trig- 
ger implications were identified as (a) mass-difference efi'ects in ratios involving the dimuon 
decays of both J/ip and ip{2S) daughters, (b) polarization effects in ratios involving both 
i^*(892)'' and K+ daughters, and (c) effects in ratios involving both of the aforementioned 
phenomena. Of the six decay topologies in the analysis, the three with the most populous 
statistics {B+ J/^pK+, B° J/^K*(892)0, and B+ i){2S)K+ (dimuonic)) were 
used to form three yield ratios representing the three topological scenarios described above. 
The relative yields for these ratios were examined for the full data sample, the subsample 
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Trigger 
Systematics (%) 


B+^J/il)K+ B^ ^ J/i)K*{d,mf B+^il){2S)K+ 
J/il) — > /i^ /^^ J/V' A*^ A*^ V'(2'S') — > /i"*"/^^ 


rO T/ih K*(RQ'?')^ 
J/V- M"*" M" 


2.0 

6.0 7.0 
2.0 2.0 6.0 
7.0 6.0 2.0 
2.0 2.0 7.0 


i^{2S) ^ 11+ fi- 


B+ V(25) K+ 


^{2S) K*{892y> 
ij{2S) fi+ 


B^ iIj{2S) K*{m2f 



Table 4.4: A summary of the systematic uncertainties (expressed in %) on the relative 
geometric and kinematic acceptance due to the effects of the dimuon triggers on the relative 
acceptance ratios. 



that passed the dimuon trigger simulation, and the subsample that failed the trigger simula- 
tion. Within the available statistics, the yield ratios remained consistent for all three trigger 
samples. Comparison of the yield ratios between the full sample and the subsample that 
passed the dimuon trigger simulation suggested a 6% systematic uncertainty for ratios with 
a dimuon ip{2S)/J/'ip mass difference, a 2% uncertainty for ratios with a i^*(892)'^/i^^ po- 
larization difference, and a 7% uncertainty for ratios with both differences. The application 



of these results to the acceptance ratios is summarized in Table 4.4. 

It is useful to investigate the sensitivity of the relative acceptance measurements to 
uncertainties in the Level 1 and Level 2 CMU and CMX trigger efficiency parameterizations. 



Table 4.5 details the effects of the one-standard-deviation parameterization shifts, pictured 



in Figures |4.3| and 4.4, on the relative geometric and kinematic acceptance values. A 



comparison of Tables 4.2 and 4.5 indicates that, for each ratio of decay-mode acceptances, 
the effect of varying the trigger efficiency parameterizations is insignificant when compared 
with the accuracy to which each relative geometric and kinematic acceptance ratio was 
determined. These variations were therefore not used in the calculation of the systematic 



uncertainties listed in Table 4.4. 
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+1(7 Trigger Shift 


B+ J/il)K+ 


B^ J/i)K*{m2f 


B+ ij{2S) K+ 


—1(7 Trigger Shift 






ij{2S) /i+ fi- 


^0 _^ J/ip K*{892)'^ 


0.403 






J / ih — > 11 


0.409 






B+ 7p{2S) 


1.139 


2.814 




Tp{2S) /u+ fi- 


1.138 


2.782 




B+ ij{2S) K+ 


0.326 


0.808 


0.287 


i){2S) J/ll^TT+TT- 


0.330 


0.807 


0.290 


^ tp{2S) K*{892)'' 


0.428 


1.061 


0.376 




0.435 


1.064 


0.382 


i;{2S) K*{892f 


0.199 


0.493 


0.175 


V'(25) ^ J/^TT+TT- 


0.202 


0.495 


0.178 



Table 4.5: A summary of the effects of +la (upper number) and —la (lower number) 
trigger efficiency parameterization shifts on the relative geometric and kinematic acceptance 
quotients. 



4.5.3 Helicity Distributions 

The longitudinal polarization fractions inherent to the S-meson decays with vector-vector 
final states can have implications for the relative acceptance measurements. The system- 
atic uncertainties associated with the longitudinal polarization fractions were estimated by 
varying the nominal measured world-average Tl/T value for the B° J/ipK*{892f mode 
(refer to Reference ||7^ and Section |4.2|) by its standard deviation, ±0.07. 

The value oiTi/T has not been measured for the decay B^ '^{2S) i^*(892)^; however, 
it is expected to be similar to that for the B^ — > J ftp K* (892)^ mode. Two standard 
deviations of the world-average Tl/T value for the B^ — > J /ipK* (892)^ measurements, 
±0.14, were therefore used as the Tl/T variation range for the B^ 'ip{2S) K*{892)^ 
modes. This variation range was applied to each of the ■0(25') i^* (892)*^ reconstructions, 
i.e., those involving '0(25') fji-^ t^-^ and ip{2S) J/tpn'^ir' decays. 



The consequent systematic uncertainties are summarized in Table 4.6. The 'absolute' 
systematic uncertainties are used in cases where the vector-vector acceptance is considered 
relative to that for one of the vector-pseudoscalar modes. 



4.5. THE ACCEPTANCE CALCULATIONS 



81 



Tl/T 

Systematics (%) 


fi^ ^ J/'iPK*{892f 


Absolute 
Acceptance 


J/ipK*{892f 

J/ip 




2.0 




4.2 


3.3 


^ 4^{2S) K*{892y 
i){2S) J/V'7r+ vr" 


2.2 


2.6 



Table 4.6: A summary of the systematic uncertainties (expressed in %) on the relative 
and absolute geometric and kinematic acceptance due to variations in the longitudinal 
polarization fractions (F^/r) of those decays with vector- vector final states. The column on 
the right lists the systematic uncertainties on the relevant absolute geometric and kinematic 
acceptance values. 



Chapter 5 

Efficiency Corrections 



The topological and candidate-selection similarities amongst the six decay channels in this 
study afforded the cancelation of several common efficiencies and systematic uncertainties 
in the branching-fraction ratios. Only those efficiencies that were unique to one of the 
decay modes in a given ratio or that differed between the channels required careful study. 
This chapter details the determination of event-yield corrections and their systematic un- 
certainties due to inefficiencies relevant to the detection and reconstruction of the candidate 
i3-meson decay modes. 

5.1 Constrained Fit Confidence Level Criteria 

The vertex and vertex-plus-mass constrained fit confidence level criteria imposed on the 
charmonium and i?-meson event selection were described in Section In the present 
section, several different techniques are used to study the effects of the CL{\^) criteria on 
the yield efficiencies. 

5.1.1 Dimuon CL{x^) Criteria 



In a branching-fraction ratio expression such as the one given in Equation 1.18| , most of the 
dimuon selection criteria, including the CL{x^) requirements, are expected to divide out 
of the calculation. It is nevertheless useful to investigate the effects of the dimuon CL(x^) 
criteria to test this assumption, examine the background rejection, and test the efficiency 
measurement techniques. 

Figures and illustrate the distribution of the CL{x^) variable for vertex con- 
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strained fits of J ftp and ip{2S) dimuon candidates, respectively. The shaded histograms 
consist of like-sign (//^ /u^) dimuon candidate combinations, normalized to the measured 
background yield. Dimuon candidates that clearly did not originate from a common vertex 
are expected to result in constrained fits with low CL{x^) values, as exemplified by the 
populous lower CL{\^) bins in the left-hand plots of Figures 5.1 and The relative 
'flatness' of the CL(x^) distributions over most of the abscissa suggests that the covariance 
matrices of the individual candidate track fits are realistic for most vertexing situations. 



Figures 5.3 and 5.4 are analogous to Figures 5.1 and |5.2| , with the difference that the 
former are due to fits with simultaneous vertex and mass constraints. Similar conclusions 
to those for the vertex-constrained fits may be drawn for the vertex-plus-mass constraint 
CL(x^) distributions. 

The CL{x^) > 0.01 criterion efficiency for vertex-constrained J/ if} candidate dimuon 
fits, may be calculated directly by examining the effect of the cut on the J/ip 

candidate yield. Figure depicts the J /ip invariant mass distribution before and after the 



application of a CL{x^) > 0.01 requirement. Similarly, Figure 5^ illustrates the analogous 
V'(2S') invariant mass distributions used in the calculation of (^^{2S)—>^i^i- 

An alternate method to measure the charmonium event yields is available by construct- 
ing normalized invariant mass distributions. Figures |5.7| and |5.8| depict the normalized 
invariant mass distributions, or 'pulls', before and after the CL{\^) criterion was imposed, 
for the J /tp and ip{2S) dimuon candidate distributions, respectively. Note that the widths 
of the Gaussian signals exceed unity in both cases, indicating some underestimation of the 



track helix uncertainties. The distributions in Figures 5.7 and |5.8| may also be used to 
estimate e}/^_^^ and e^^^s)^^^- 

The extra inefficiency arising from the additional constraint to the appropriate world 
average charmonium mass, e™, may be estimated from the normalized invariant mass dis- 
tributions. By relating the peak area obtained with a double-Gaussian fit in a ±2.58 
standard deviation window about the mean of the normalized invariant mass signal (with 
the C-L(x^) > 0.01 criterion on the vertex-only constrained fit already applied) to the entire 
signal fit area, the efficiency of an additional mass-constraint CL(x^) > 0.01 requirement 



can be computed. The right-hand plots in Figures 5.7 and |5.8| were used in this way to 
measure ej^^^^^ and £^(25)^/./. > respectively. 
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Figure 5.1: The CL{x'^) distribution of the vertex-constrained J/ijj dimuon fit. The un- 
shaded histogram in the left plot contains the combined signal and background events; the 



shaded histogram in the left plot represents the like-sign dimuon {fi ^ ) candidates, nor- 
malized to the background yield. The right plot is the difference of the two histograms in 
the left plot. 




Figure 5.2: The CL(x^) distribution of the vertex-constrained ip{2S) dimuon fit. The 
unshaded histogram in the left plot contains the combined signal and background events; 
the shaded histogram in the left plot represents the like-sign dimuon (/x^ /x*) candidates, 
normalized to the background yield. The right plot is the difference of the two histograms 
in the left plot. 
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Figure 5.3: The CL{x^) distribution of the vertex-plus-mass-constrained J/V' dimuon fit. 
The unshaded histogram in the left plot contains the combined signal and background 

events; the shaded histogram in the left plot represents the like-sign dimuon (/v,^ /v,^) can- 
didates normalized to the background yield. The right plot is the difference of the two 
histograms in the left plot. 




Figure 5.4: The CL(x^) distribution of the vertex-plus-mass-constrained ip{2S) dimuon 
fit. The unshaded histogram in the left plot contains the combined signal and background 
events; the shaded histogram in the left plot represents the like-sign dimuon {fi^ n^) can- 
didates normalized to the background yield. The right plot is the difference of the two 
histograms in the left plot. 
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Figure 5.5: The J ftp dimuon candidate invariant mass distribution before (left) and after 
(right) a vertex-constraint criterion of CLi^) > 0.01 was apphed. The fit is to a double 
Gaussian signal with linear amplitudes and a linear background term. 




Figure 5.6: The tp{2S) dimuon candidate invariant mass distribution before (left) and after 
(right) a vertex-constraint criterion of CL{x^) > 0.01 was applied. The fit is to a Gaussian 
signal and a linear background term. A requirement of pT > 2.5 GeV/c was imposed on 
each muon candidate to control the background. 



5.1. CONSTRAINED FIT CONFIDENCE LEVEL CRITERIA 



87 




Figure 5.7: The J/V' dimuon normalized candidate invariant mass distribution before (left) 
and after (right) a vertex-constraint criterion of CL{x^) > 0.01 was applied. The fit is to a 
double Gaussian signal with linear amplitudes and a Gaussian background term. 




Figure 5.8: The ip{2S) dimuon normalized candidate invariant mass distribution before 
(left) and after (right) a vertex-constraint criterion of CL{x^) > 0.01 was applied. The fit 
is to a Gaussian signal and a Gaussian background term. A requirement of pi > 2.5 GeV/c 
was imposed on each muon candidate to control the background. 
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Efficiency 


Technique 


Measurement 


gV 
gV 


mass 
normalized mass 

mass 
normalized, mass 


0.966 ± 0.003 
0.969 ± 0.004 
0.956 ± 0.036 
922 ± 039 




weighted mean 


0.967 ± 0.003 


m 

m 


normalized mass 
normalized mass 


0.963 ± 0.002 
0.958 ±0.035 


m 


weighted mean 


0.963 ± 0.002 



Table 5.1: Dimuon vertex and vertex-plus-mass CL(x^) > 0.01 criterion efficiency summary. 
The quantities e^^ and e™^ were computed by performing a weighted mean of the measured 
values. 



Table 5.1 presents a summary of the various CL(x^) criterion efficiencies measured in 
this section. As expected, the various measurements, using either J/^lJ or ip{2S) dimuon 
reconstructions, are largely mutually consistent, thereby supporting the assumption that 
dimuon CL{x^) selection effects divide out of the branching-fraction ratio calculations. The 
average vertex and mass constraint requirements were measured to be e^^ = 0.967 it 0.003 



and e™, = 0.963 ± 0.002, respectively. 



MM 



5.1.2 Dipion CL{x^) Criteria 

Unlike in the dimuon cases, charmonium reconstructions performed via the decay ip{2S) 
J/tpTT^TT^ did not necessarily have CL{x^) requirement efficiencies that divided out of 
the branching-fraction ratios. The vertex and vertex-plus-mass CL{\^) cut efficiencies for 
this hadronic-cascade charmonium mode were taken to be factorizable according to the 
expressions e;(2s)-.j/^^^ = ^It^ ' ^In and e;^(25)-^j/v^^ = ' respectively. The fact 
that the J/ip dimuon candidate was first vertex-plus-mass constrained separately, before 
the vr"^ vr" candidates were included in the fit, permitted this factorization of CL(x^) cut 
efficiencies. The principal advantage of factorizing the efficiency in this way is that the 
dimuon efficiency factors, e^^ and ej^, can still divide to unity in the branching- fraction 
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Figure 5.9: The CL{x^) distribution of the vertex-constrained vr'^ vr^ fit. The unshaded 
histogram in the left plot contains the combined signal and background events; the shaded 
histogram in the left plot represents the like-sign dipion (vr^vr^) candidates normalized to 
the background yield. The right plot is the difference of the two histograms in the left plot. 



ratio expressions. 



Figure 5.9 illustrates the distribution of the CL{x^) variable for vertex-constrained fits 
of J/ipTT^TT~ candidates performed after vertex-plus-mass constrained fits (and the corre- 
sponding CL{x^) criteria) were applied to the J/ip dimuon candidates alone. The shaded 
histogram was created from like-sign-pion {J/iJj vr^ vr^) candidate combinations, normalized 
to the measured background yield. The features in Figure are similar to those observed 
in Figures 5.1 and 5.2 in Section 



Figures 5.10 and 5.11 respectively depict the J/iljTT~^Tr~ invariant mass and normalized 
invariant mass distributions before and after the additional CL{x'^) > 0.01 criterion was 
applied to the 7r+ vr^ vertex-constrained fit quality. The techniques used to measure the ej^^ 
and e™jr efficiencies were similar to those described in Section |5.1.1| . Table ^.2| summarizes the 
e^^ and e™^ results. The uncertainties on the efficiencies were estimated by conservatively 
combining in quadrature the yield uncertainties before and after the CL{x^) > 0.01 criterion 
was imposed. 
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Figure 5.10: The J / ip Tr~^ Tr~ candidate invariant mass distribution before (left) and after 
(right) a vertex-constraint criterion of CL{x^) > 0.01 was apphed. The fit is to a Gaussian 
signal and a cubic polynomial background term. 




Figure 5.11: The J / tp n'^ Tr~ normalized candidate invariant mass distribution before (left) 
and after (right) a vertex-constraint criterion of CL{\^) > 0.01 was applied. The fit is to a 
Gaussian signal and a cubic polynomial background term. 
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Efficiency 


Technique 


Measurement 


'-■KTT 


mass 
normalized mass 


0.816 ±0.039 
0.858 ± 0.045 


TVTT 


weighted mean 


0.834 ± 0.039 




normahzed mass 


0.945 ± 0.031 



Table 5.2: Dipion vertex and vertex-plus-mass CL{x^) > 0.01 criterion efficiency summary. 
5.1.3 5-Meson CL{x^) Criteria 

The efficiency of the vertex and momentum-pointing CL{x^) > 0-01 criterion on the candi- 
date B decays was expected to be similar to those calculated for the inclusive charmonium 
reconstructions. In practice, limited statistics made it difficult to examine the effects of 
CL{x^) requirements on 'global' vertex-constrained fits of S-meson candidates. To accom- 
modate any possible noncancelation of the global B CL{x^) cut efficiencies for a given ratio 
of branching fractions, a systematic uncertainty of 2.0% was assigned on each ratio. This 
systematic uncertainty was estimated by taking the standard deviation of the measured 



and e™ values in Table 5.1 



5.2 Tracking Efficiencies 

An advantage of measuring ratios of branching fractions was that the tracking efficiencies 
largely divided out of the calculations. Such was the case for ratios involving the dimuonic 
charmonium modes and like species of strange meson in both numerator and denomina- 
tor. Ratios of branching fractions involving the decay ip{2S) J/ipir^n^ and unlike 
kaon species of the form K*{^Q2)^ /K^ , however, possessed residual dipion and single-pion 
tracking efficiency factors, respectively. 

Over the course of the entire Run 1 period, several known effects with potential implica- 
tions for tracking performance were noted: the instantaneous luminosity varied by amounts 
in excess of an order of magnitude (refer to Section |3.1[) , the propensity for multiple pri- 
mary vertices in a given event increased, the Tevatron beam optics were modified, the wire 
geometry in the CTC was altered by gravity-induced sagging, and the electrostatic and 
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gas-quality characteristics of the CTC changed. 

Motivated not only by the requirements of the present ratio-of-branching-fraction anal- 
ysis, but also by the needs of other CDF analyses that made use of low momentum tracking 



information, a detailed study [135] of single- and double-track CTC pattern recognition 
efficiencies was conducted. Appendix |^ summarizes the methodology and results of the 
pattern recognition efficiency measurements. The single-pion pattern recognition efficiency 
was measured to be 0.928 it 0.020 for > 0.4 GeV/c and r^y^c > ^0 cm. The double-pion 
efficiency was measured to be 0.881 it 0.043, also for > 0.4 GeV/c and v'q^^q > 110 cm. 

5.3 Proper Decay Length Criteria 

In this analysis, both the numerator and the denominator in each ratio of branching fractions 
contained or meson candidates. The measured world-average proper decay lengths 
(ctb) of and B^ mesons are (492 it 18) ^m and (465 it 15) ^m |3^, respectively. The 
fact that the accepted lifetimes for the two S-meson flavours were mutually consistent lead 
to the expectation that the ctb > 100 fim criterion (refer to Section [3.9.2) efficiencies would 



be similar for all the modes. The world-average B lifetimes, however, were short enough 
that the efficiency of a 100-/im proper decay length cut was potentially sensitive to the 
resolution of the ctb measurement itself, which, in turn, was sensitive to the momentum 
and vertex resolution of the final-state tracks. Here the different decay topologies could 
play a measurable role in the efficiency, owing to the ctb resolution dependence on the 
multiplicity of candidate tracks that made use of SVX (Section |2.3.1 ) tracking information. 



The sideband regions of the 5-meson candidate invariant mass distributions were ex- 
pected to be populated primarily by prompt, i.e., low proper-decay-length, processes. The 
regions 5.145 — 5.220 GeV/c^ and 5.340 — 5.415 GeV/c^, with the kinematic selection cri- 
teria relaxed to enhance the statistics, were therefore used to examine the ctb resolution 
as a function of the number of daughter-candidate tracks with SVX hit information used 
in their fits. The results for the B'^ and B^ decay modes are given in Tables and 5.4, 
respectively. 

The relative frequencies of occurrence of the different SVX multiplicity possibilities were 
estimated for each decay mode with a Monte Carlo calculation (described in Chapter ^), for 
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Number of 
tracks with 
SVX hits 


J/i) n+ 11- 


i}{2S) n+ II- 


B+ V(25) K+ 

j/i) 11+11- 





575 ±6 


783 ± 49 


520 ± 26 


1 


171 ±4 


181 ± 1 


189 ± 35 


2 


94 ± 1 


99 ±3 


243 ± 122 


3 


57 ± 1 


52 ± 1 


55 ±9 


4 






78± 10 


5 






49 ±2 



Table 5.3: Sideband ctb resolutions, measured in units of /xm, for each B+ decay mode and 
SVX track multiplicity. 



Number of 

tracks with 
SVX hits 


^ J/^K*(892)" 

J/ll) 11+ IX- 

ii:* (892)0 _^ K+TT- 


B"^^{2S) K*{892)" 
%1){2S) 11+11- 
iC* (892)0 ^ K+TT- 


B^^^P{2S) if* (892)0 

i;{2S) J/i^TT+TT- 
J/ip n+ IX- 
(892)0 ^ i;s:+7r- 





591 ±3 


762 ± 22 


584 ± 13 


1 


202 ±3 


210 ±7 


291 ± 52 


2 


117 ± 1 


102 ±3 


168 ± 21 


3 


68 ± 1 


68 ± 1 


74 ±8 


4 


49 it 1 


56 ± 1 


61 ±4 


5 






59 ±2 


6 






48 ± 1 



Tabic 5.4: Sideband ctb resolutions, measured in units of /xm, for each B decay mode and 
SVX track multiplicity. 



94 



CHAPTER 5. EFFICIENCY CORRECTIONS 



B Mode 


cc Mode 


Central 


-la 


+l(J 


Efficiency 


B+ - 


J/i) K+ 




.753 


.741 


.764 


.753 ± .012 


B+ - 


il>{2S)K+ 




.738 


.726 


.748 


.738 ± .013 


B+ - 


il>{2S)K+ 




.762 


.750 


.774 


.762 ± .012 


B^- 


J /%bK* (892)° 




.755 


.742 


.766 


.755 ± .013 


B^- 


■> V(25) i^*(892)0 


ip{2S) 11+ fi- 


.749 


.736 


.760 


.749 ± .013 


B^- 


■> V(25) i^* (892)0 




.761 


.748 


.773 


.761 ± .013 



Table 5.5: Measured ctb cut efficiencies in the B+ and B decay modes, showing the effects 
of variations of the world-average B lifetimes by their standard deviations. 



which the detector simulation had been tuned on data to determine the SVX hit efficiencies 



and hit association characteristics. The measured ctb resolutions in Tables 5.3 and 5.4 
were then convolved with ctb distributions that used the world average B'^ and B^ proper 
lifetimes in their decay constants, respectively. A toy Monte Carlo programme was 
employed to perform 10^ convolutions for each decay, and the resolution contributions from 
the various SVX track multiplicities were administered according to weighting probabilities 
determined by the relative frequencies calculated from the Monte Carlo samples. The 
procedure was then repeated for the world-average ctb values shifted high and low by their 
standard deviations (cj) |^^. The results of these shifts, along with the central values. 



are listed in Table 5.5. The systematic uncertainties on the efficiency values correspond 
to the maximum efficiency deviations observed in the —la and +la studies. As expected, 
the measured ctb efficiencies were similar, with all values near 0.75. The small efficiency 
differences amongst the various channels were consistent with differences in the daughter- 
track multiplicities. 



5.4 5-Candidate Isolation Criterion 



The 5-meson candidate isolation criterion, Ib 



R 

^ , ^ , < — , was defined and de- 

\pb\ 13' 

scribed in Section 3.9.2 . This requirement rejected significant contributions of background 
to the -B-candidate signals, namely backgrounds arising from invariant mass combinations 
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B Mode 


cc Mode 


Efficiency 


7~>-l- 1 / n {~^\ ry4- 

B^ 1^(25 ) 
B+ Tp{2S) K+ 


■0(26) fi^ ^ 

-0(25) ^ J/V^TT+TT- 


0.900 ± 0.063 

1 nnn +0.000 

1-000-0.300 
1 000+°°°° 


B^ J/V'K* (892)0 
5° ^ ^(25") (892)0 
50 ^ ^(2S) (892)0 


-0(25) ^ ^+^- 

-0(25) ^ J/V^TT+TT- 


n qQ7 +0.063 

n/a 
9 07 +0 093 

u.»U( _o.268 


Combined-Channel Measurement 


0.928 ±0.054 



Table 5.6: Measured isolation criterion efficiencies in the i?+ and i?° decay modes. The 
domination of background in the i?° tp{2S) i^* (892)0, ij){2S) //+ /i~ decay-chain recon- 
struction prevented an individual measurement; however, the invariant mass distributions 
for all six modes were used to compute the combined efficiency in the bottom row. 



of the dimuon candidates with non-i? tracks from the underlying event (which is repre- 



sented by X in Equation lA). The large rj-ip cone radius and the stipulation that the 
non-i? tracks that contribute to the summation in the isolation expression be consistent 
with having originated from the same primary vertex as the B candidate suggested that the 



isolation criterion would be relatively efficient, as exemplified in Figure |3.1| . Investigations 
of S-meson isolation requirements indicate that the efficiencies of these criteria are largely 
independent of the type of B decay mode and the pT of the B candidate |136]. 



As a consistency test of the observation [ 136 | that the i?-candidate isolation criterion 
was independent of decay topology (therefore obviating a need for individual efficiencies in 
the branching- fraction ratio calculations), the effective efficiency of the isolation require- 
ment was examined for each of the six S-meson reconstructions investigated in this study. 



Table 5.6 lists the efficiencies, which were measured by comparing the appropriate candidate 
yields after and before the isolation criterion was imposed. The combined efficiency of the 
isolation requirement was found to be 0.928 it 0.054 and the efficiencies for the individual 
topologies were observed to be mutually consistent. 
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Decay Mode 


Branching Fraction 


ip{2S) /i+ fi- 
ip{2S) J/tPtt+tt- 


(6.01 ±0.19) X 10-2 

(8.5 ±0.7) X 10-3 
(3.07 ±0.19) X 10-1 


(892)0 ^ i^+vr- 


2/3 



Table 5.7: Branching fractions of the daughter meson decay modes that were reconstructed 
in the present analysis. For the charmonia, the world- average branching fractions were 



used, with the assumption of lepton universality in the B{'ip{2S) fj. ) case [3C]. The 

K*{892)^ branching fraction was based on isospin symmetry. 



5.5 Daughter Branching Fractions 

Depending on the ratio of i?-meson branching fractions being calculated in this study, dif- 
ferent combinations of daughter branching fractions did not divide out of the calculation. 
These quantities are listed in Table |5.7| . The charmonium branching fractions were taken 
to be the most recent accepted world-average values |Q; however, the B{J/ip A*^ /t^ ) 
average did not take into account the measurement reported in Reference [|137[| . For the 
value of B{ip{2S) — > fi^ fi~), the world average for B{Tp{2S) — > e^e~), which has a sig- 
nificantly smaller uncertainty, was used under the assumption of lepton universality. The 
B{K* (892)^ 7r~) value is based on isospin symmetry. The uncertainties on the world- 

average charmonium branching fractions constitute the largest contributions to the total 
systematic uncertainties in the present analysis. 



Chapter 6 

Branching Fraction Measurements 



6.1 B-Meson Candidate Event Yields 

The S-meson candidate event yields in this study were computed by fitting binned invariant 
mass distributions to signal and background lineshapes using a maximum log-likelihood 
technique. In order to avoid any contributions of partially reconstructed B mesons to the 
lower-sideband fit region, invariant masses <5.15 GeV/c^ were excluded from all the fits. 



Figure 6A. features the observed J /if) {J /'>P fJ-^ M ) invariant mass distribution. 

Figure |6^ illustrates the analogous B'^ V'(25') invariant mass distributions, with both 
the il^{2S) Ai"*" M~ ip{2S) — > J/'ijj'K^ tt^ charmonium daughter modes. 

The fits for the B^ candidate event yields embodied the relative double-Gaussian pa- 



rameterizations described in Table 3T in Section 3.9.3 . Figure O shows the J/tp K*{892)^ 
{J/ip — > /U"*" invariant mass distribution, whereas Figure |6.4| presents both the ip{2S) 
fi~ and ip{2S) J/tpir'^ tt~ reconstructed iIj{2S) i^*(892)'^ invariant mass distributions. 



The event yields and their statistical uncertainties are summarized in Table 6.1. The 
relative yields of those final states with ip{2S) candidates are consistent with expectations 
from the acceptance and efficiency measurements. An uncertainty of 4.0% was assigned on 
the branching-fraction ratios to account for differences in the event yields due to different 
fitting intervals and different right-mass-combination counting techniques in the i^*(892)^ 
modes. The fit interval was varied from 5.1 — 5.6 GeV/c^ to 5.2 — 5.6 GeV/c^ and an algo- 
rithm that selected vr^ combinations nearest the world-average K*{892)^ pole mass pO| 
were used in lieu of the default techniques (fitting over the interval 5.15 — 5.6 GeV/c^ and 
using the double-Gaussian method) to estimate this systematic uncertainty. 
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Figure 6.1: The J/ipK~^ invariant mass distribution. The fit is to a Gaussian signal Hne- 
shape and a linear background parameterization. 
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Figure 6.2: The ip{2S) invariant mass distributions using the ip{2S) ji'^ (left plot) 
and ip{2S) — > J/ipTT~^Tr~ (right plot) modes. Each fit is to a Gaussian signal lineshape and 
a linear background parameterization. 
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Figure 6.3: The J/ip K*{892)'^ invariant mass distribution. The fit is to a double Gaussian 
signal lineshape and a linear background parameterization. The indicated fit parameters 
describe the Gaussian distribution of those combinations with correct K-n mass hypotheses. 
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Figure 6.4: The i>i2S) K*{892f invariant mass distributions using the '4^{2S) (left 
plot) and ■0(25') J/ijj 7r~^ Tr~ (right plot) modes. Each fit is to a double Gaussian signal 
lineshape and a linear background parameterization. The indicated fit parameters describe 
the Gaussian distribution of those combinations with correct K-ir mass hypotheses. 
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B Mode 


cc Mode 


Event Yield 


B+ J/'4)K+ 
B+ tl;{2S) K+ 
B+ il){2S) K+ 


V'(25) n+ fi- 
ip{2S) J/tPtt+tt- 


856.7 ±38.3 
71.9 ± 13.4 
37.4 ± 7.4 


B° J/tpK*{892)^ 
B^ ^^{28) if* (892)0 
B^ 7p{2S) ET* (892)0 


ip{2S) 11+ fi- 
ip{2S) J/iPn+TT- 


378.8 ± 24.8 
20.9 ± 7.3 
29.1 ± 7.5 



Table 6.1: Summary of S-meson candidate event yields. The indicated uncertainties are 
statistical only. 

6.2 Sample Branching-Fraction Ratio Calculation 

The event yields, relative geometric and kinematic acceptances, efficiencies, daughter-meson 
branching fractions, and systematic uncertainties were combined into the calculations of the 
ratios of branching fractions. An example of the calculation of the ratio -Rj/^'^*^892)o ™ 



Table l.S can be illustrated by expanding the form of the expression in Equation 1.18 to 
treat both of the tp{2S) daughter reconstructions: 

B{B+ i;{2S) K+) 



(6.1) 



B{BO J/^K*{892)0) 

N{B+ ^i;{2S)K+,'^{2S) ^ fi+ fj-) + N{B+ ^^{2S)K+,'^{2S) ^ J/'^ir+TT~) ^ 

N{BO J/^P K*(892)0, J/iP n+ fi-) ' 

£tot(5+ ^ ^(25)E:+,?/'(25) ^ /i+^-) + e^°^{B+ ^^P{2S)K+,%b{2S) ^ J/ipTT+TT- 
where the symbols represent the event yields and the quantities represent products of 
absolute geometric and kinematic acceptances, efficiencies, and daughter-meson branching 



fractions^ . Resolving the expression in Equation 6T into the relevant measured factors gives 
the form 

B{B+ il>{2S)K+) _ 
B{B^ J/V' if* (892)0) - 

N{B+ V(25) K+,iP{2S) 11+ I1-) + N{B+ tP{2S) K+,iP{2S) J/^ tt+ vr") 

A^(^o ^ j/^ iC*(892)0, J/^ ^ fi+ ^-) ^ 



^The method of combining the results from the two ^(2S) channels was motivated by the low efBciencies 
involved in identifying candidate events, thereby requiring the use of Poisson probability distributions. 
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gtrk . ^cr^^o ^ j/^K*(892)°) • e(i^*(892)0 ^ i^+TT-) X 



A{B^ J/4^K*{892f 



+ 



A{BO J/V'K*(892)C 



.trk 



-1 



where e^'"'^ represents the single-pion pattern recognition efficiency (refer to Section |5.2| and 
Appendix |C|), e'^'^ represents the proper-decay-length requirement efficiencies (refer to Sec- 
tion |5.3| ), the quantities A represent the geometric and kinematic acceptances, e^.^ and e™^ 
represent the CL{x^) requirement efficiencies (refer to Section 5.1.2| ), and e^^^ represents the 



double-pion pattern recognition efficiency (refer to Section 5.2). The numerical ingredients 



entering into the calculation of Equation |6.2| are listed in Table 6.2. 



6.3 Branching- Fraction Ratio Results 



The results for the ratios of the four S-meson final states are presented in Table 6.3. Where 
appropriate, the ^{2S) fJ-^ and ip{2S) — > J / tp ir'^ 7r~ contributions were consolidated 

to produce single measurements using approaches similar to that of the sample calculation 



in Section 6.2. In ratios that contained unlike species of B meson in the numerator and 



denominator, the equality fu = fd was assumed (refer to Section |1.3|). It is interesting to 



note from Table 16. 31 that R 



M2S)K+ 
J/V' K+ 



^jf^K^{m2'f ' which suggests that the flavour of the 



light spectator quark (refer to Figure 1.3) has little effect on the decay width. 

Figure compares the measured branching-fraction ratios with predictions derived 



using two different factorization approaches. As discussed in Section 1.2.3 , several factors 
divide out of the theoretical expressions for the decay-rate ratios, including the CKM ele- 



ments and the parameter a2 in Equation 1.15 . The Neubert et al. [ 138 | predictions used 
the same form factors as BSW which were calculated using a relativistic harmonic 

oscillator potential model; however, whereas the BSW form factors were given a single-pole 



dependence on q in References [73, 74], the Neubert et al. |13^ ] predictions were based 



on most of the form factors having a dipole behaviour. As mentioned in Section 1.2.4, the 
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Quantity 


Value 


N [B^ ^ ip^zb) ,^j[lb) ^ 11 ) 
N{B+ ^j{2S) K+,%bi2S) J/%bTT+TT-) 
N{B^ J/V' i^*(892)0, J/ib fi+ 11-) 


71.9 ± lo.4 
37.4 ± 7.4 
378.8 ± 24.8 


A{B+ ^(25) K+,ij{2S) fi+ fi-) 
A{BO J/%b K*(892)0) 

A(R+ -^ih('?'^\ 7/;('9S"l 7r+ TT^^I 

A{BO J/V'K*(892)0) 


2.810 ±0.068 
0.813 ±0.025 


ecr/^+ ^ ^(25) K+,^h(2S) ii+ u,-) 
e^^{B+ i){2S)K+,i>{2S) J/V'vr+vr-) 
gCr^^O ^ (892)0) 


0.738 ±0.013 
0.762 ±0.012 
0.755 ± .013 


trk 

TTTT 


0.928 ± 0.020 
0.881 ± 0.043 




0.834 ± 0.039 
0.945 ± 0.031 


B{^Pi2S) ^ fi+ fi-) 
B{il;{2S) J/V^TT+vr-) 
B{J/^P ^ fi+ 

Oyi\ yov-Lj — > J\ 71 j 


(8.5 ±0.7) X 10-3 
(3.07 ±0.19) X 10-1 
(6.01 ±0.19) X 10-2 


6-Quark Production Spectrum Systematics 

Trigger Systematics 

Helicity Systematic 

Simulation Systematic 

Fitting Procedure Systematic 

Global CL^x^) Systematic 


5.9%, 5.0% 
7.0%, 2.0% 

2.0% 

3.0% 

4.0% 

2.0% 



Table 6.2: Summary of event yields, geometric and kinematic acceptances, efficiencies, 
daughter-meson branching fractions, and systematic uncertainties used in the sample cal- 



culation of R 



J /ipK* (892)0 



by way of Equation 6.2. 
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Branching-Fraction Ratio 



Measurement 



R 



R 



J/^K* {892)0 
J/ipK+ 

J/ipK+ 

tp{2S)K*{892Y' 
J/i)K+ 



R 



R 



R 



R 



.i'(2S) K+ 
J/V'ii'* (892)0 

■il){2S)K*{892f 
J/ipK* {892)0 

ij{2S)K* {892)0 
i!{2S) K+ 



1.76±0.14[stat] ±0.15[syst] 
0.558 ± 0.082[stat] ± 0.056[syst] 
0.908 ± 0.194[stat] ± 0.100[syst] 
0.317 ± 0.049[stat] ± 0.036[syst] 
0.515 ± 0.113[stat] ± 0.052[syst] 
1.62 ±0.41 [stat] ±0.19[syst] 



Table 6.3: The branching-fraction ratios measured for the various S-meson final states, 



corresponding to the elements of Table 1.3. Note that the ratios containing il^{2S) mesons 
are composed of contributions from two separate i?-meson reconstructions. The first un- 
certainties are statistical and the second are systematic. 



Deandrea et al. |8^ predictions employed monopole form factors that were derived from 
D-meson decay measurements. 



Inspection of Figure |6.5| indicates that neither of the two phenomenological approaches is 
discounted by the data. The agreement between prediction and data, however, is favoured 



somewhat in the Neubert et al. |138] case, for which all but one of the predictions are 
consistent with the measured ratios to within one standard deviation. The Deandrea et 
al. [^l| predictions are only consistent within one standard deviation of the measured values 
for the last two ratios in Figure |6.5| ; in the first four entries, the Deandrea et al. predictions 
agree with the data within two standard deviations. Note that Deandrea et al. predict a 
somewhat lower — > ip{2S) rate, which accounts for much of the observed difference. 

Although the consistency between data and prediction in Figure |6.5| does not constitute 
proof for the validity of the factorization Ansatz in colour-suppressed i?-meson decays, the 
data do appear to favour the approach where the majority of the form factors have a dipole 
dependence in q"^ and where semileptonic D-meson measurements were not used to estimate 
the heavy-to-light form factors. 
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1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 

• B{B^ ^ J K* (892)^) 

B{B+ ^J/tPK+) 


O A 


B(B+^V(2S) K+) 
B{B+^J/iP K+) 


• 

O A 


B(b0 _^ ■4,{2S) ii:*(892)0) 
B(B+ — >■ JM K+\ 


O A 


B(B+->V(2S) i4r+) 
B{B^-^Jlt\) if*(892)0) 


— • — 

O A 

1 1 1 1 1 1 


S(B0_^^(2S) A:*(892)0) 
B(BO^J/-0 if* (892)0) 

• 6(50^ -0(25) ii:*(892)0) 

A O B(B+^^(2S)iir+) 
1 1 1 1 1 1 1 1 1 1 1 



0.0 1.0 2.0 3.0 

Ratio of Branching Fractions 



Figure 6.5: Comparison of the measured branching-fraction ratios (filled circles) with theo- 
retical predictions by Neubert et al. [|138(| (triangles) and Deandrea et al. [^l| (diamonds). 
The CDF measurements were taken from Table |6.3| , with the statistical and systematic 
uncertainties combined in quadrature. In ratios involving both and mesons, = fd 
was assumed. 
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6.4 Derived Absolute Branching-Fraction Results 

The measured ratios of branching fractions listed in Table were used to extract absolute 
branching fractions for three of the four final states by employing the world average value 
for B{B^ — > J/t/j K^), which is the best measured quantity^. The technique of normalizing 
to a single world-average branching fraction was adopted in the present study for several 
reasons: a weighted-mean approach would have required a more complicated method of 
combining correlated systematic uncertainties, the world-average values for those modes 
other than B~^ — > J/^ had comparatively large uncertainties, no previous measurements 
using the CDF data sample were used in the determination of the world average value of 
B{B^ J/tl^K^), and the extracted absolute branching-fraction results derived using a 
single world-average normalization could be updated easily with subsequent improvements 
to the world-average value of B{B^ — > J/ip K~^). 
The world-average value [^] 

B{B+ J/iP K+) = (1.01 ± 0.14) X 10^3 (6.3) 



was therefore used with the first three entries in Table |6.3| to compute the derived abso- 
lute branching fractions listed in Table |6.4| . The calculations of the 'branching fraction' 
uncertainties for the two modes involving ip{2S) daughters assumed that the world-average 
charmonium branching-fraction uncertainties were independent of the uncertainty in Equa- 



tion 6.3; this was the most conservative assumption. 

The derived absolute branching fractions for the decays B^ — > ipC^S) and B^ 
ip{2S) K* (892)^ are compared with measurements and limits from colliders in Fig- 

ure |6.6| . These results are the world's most precise measurements of these branching frac- 
tions. In the case of the mode B~^ iJj{2S) , the result from the present study is 
consistent within one standard deviation of the measurement reported by the CLEO II 
collaboration[|139|] , but is only marginally consistent with the measurement reported by the 



ARGUS collaboration |T4|. For the B° ^{28) K*{m2f channel, the CDF result is seen 



to be consistent with the previous measurement and limitsQ. 



The structure of the matrix in Table 



1.5 , however, suggested that an alternate approach could have been 
effected by performing weighted averages of multiple ratio measurements and normalizing to world-average 
values for the modes other than the one for which an absolute value was to be derived. This approach was 
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T 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I 



ARGUS 1990 



CLEO 1992 (90% CL) 



CLEO 1994 



CDF 



i/{2S) K*(892y 



ARGUS 1990 (90% CL) 

CLEO 1992 

CLEO 1994 (90% CL) 



CDF 



J I I 



J I I 



J I I 



10"' 



1 0'' 1 0'' 

Branching Fraction 



10'^ 



Figure 6.6: A comparison of the derived CDF B{B+ ^{23) K+) and 5(5° 
il)(2S) K* {S^2)^) absolute branching fractions with measurements and limits from the AR- 



GUS CLEO lUH, and CLEO II |3|] experiments. The hatched bars denote 90% CL 

upper hmits and the error bars represent the statistical, systematic, and branching-fraction 
uncertainties added in quadrature. 
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Branching Fraction (B) 


CDF Measurement [xlO^^] 


PDG [xlO-^] 


B° J/?/; (892)0 

B+ ip{2S) K+ 

B° i){2S) (892)0 


1.78±0.14[stat] ±0.29[syst] 
0.56 ± 0.08[stat] ± 0.10[syst] 
0.92 ± 0.20[stat] ± 0.16[syst] 


1.49 ±0.22 
0.69 ±0.31 
1.4 ±0.9 



Table 6.4: The derived CDF absolute branching fractions calculated by relating the entries 



in Table |6^ to the world-average value in Equation |6^. The first uncertainties are statistical 
and the second are systematic (including the uncertainty in Equation 3.S and, if applicable, 
the uncertainties in the world-average values of the charmonium branching fractions [^0| ). 
The right column lists the previous world-average values compiled by the Particle Data 
Group (PDG) 1^]. 



The derived absolute B ipi^S) branching fractions may also be compared against 
the predictions detailed in Table |1.2| , as is illustrated in Figure |6.7| . Whereas the less 
recent Deshpande and Trampetic |71] predictions do not lie within one standard deviation 
of the measured CDF values, the Deandrea et al. [^] and Cheng et al. |86| predictions are 
consistent with the measurements. The general agreement amongst the measurements and 
predictions in Figure |6.7| is relatively good in spite of the large model dependencies and 
uncertainties in the predicted values. 



taken in References 121]. 

^The reversion of the CLEO result from a measurement to a limit (refer to Figure 
to statistical effects |141|. 



6.6) has been attributed 
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T 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I 

f (2S) K'' 

▲ Deshpande et al. 1990 



Deandrea et al. 1 993 



Cheng et al. 1997 



CDF 



B° ij{2S) K*(892)° 



Deshpande et al. 1990 



Deandrea et al. 1993 



Cheng et al. 1997 



CDF 



J I I 



J I I 



J I I 



10 



-4 



1 0"' 1 0"' 

Branching Fraction 



10 



-1 



Figure 6.7: A comparison of the derived CDF B{B+ i;{2S)K+) and BiB'^ 
i/j{2S) K* {892)^) absolute branching fractions with theoretical predictions due to Desh- 



pande and Trampetic |71|, Deandrea et al. |81|, and Cheng et al. ||8^. The error bars in 
the Deandrea et al. results represent only the uncertainties on the experimental inputs to 
their form factor estimations. 



Chapter 7 

Conclusions 



Exclusive nonleptonic decays of B mesons into strange mesons and J/ip or tp{2S) mesons 
were studied using the CDF detector to observe a sample corresponding to JCdt = (109 it 
7) pb~^ of 1.8-TeV proton-antiproton collisions produced by the Fermilab Tevatron collider. 
The investigated decays were 

B+ J/ip K+ 

B+ ^{2S) K+ 

5° ^ ^{2S) (892)0, 

which were reconstructed via the daughter meson decay modes 

J/tJj n+ li- 

K*{892f ^K+TT-. 

With statistical significances of ~5 and ~3 standard deviations, respectively, the obser- 
vations of the decay modes B+ tp{2S) K+ and 5° iIj{2S) K*{892f constituted the 
first measurements of these processes in a hadron collider environment and amounted to 
the world's largest recorded sample of these decays. 
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7.1 Branching Fraction Measurements 

Several similarities between the decay modes were exploited to make possible the precise 
measurement of the relative branching fractions between the i?-meson decay modes studied. 
The calculated relative geometric and kinematic acceptances, measured efficiency correc- 
tions, observed yields, world-average daughter-meson branching fractions, and systematic 
uncertainties were used to compute the ratios of branching fractions. For ratios that in- 
volved both and B° mesons, equal production rates were assumed for both meson types. 
The measured ratios are 

1.76 ± 0.14[stat] ± 0.15[syst] (7.1) 
0.558 ± 0.082[stat] ± 0.056[syst] (7.2) 
0.908 ± 0.194[stat] ± 0.100[syst] (7.3) 
0.317 ± 0.049[stat] ± 0.036[syst] (7.4) 
0.515 lb 0.113[stat] ± 0.052[syst] (7.5) 

1.62 ± 0.41[stat] ± 0.19[syst], (7.6) 

where the first uncertainties are statistical and the second are systematic. 

The branching-fraction ratio measurements were compared with two different sets of 
theoretical predictions that each employed the factorization Ansatz, but that used form fac- 
tors derived in different ways with different paranietcrizations in q'^. All of the predictions 
were consistent with the measured ratios to within two standard deviations, thereby sup- 
porting, but not proving, the applicability of the factorization Ansatz to colour-suppressed 
B ^ tp K decays. 

Of the two sets of predictions, one was consistent with the measurements to within 



B{B° J/V^i^* (892)°) 
B{B+ J/i)K+) 

B{B+ V;(25) K+) 
B{B+ J/il)K+) 

B{B^ i){2S) K*{892f) 



B{B- 


^ J/^K^) 


B{B- 


^ ^ V(25) K+) 


B{B^ - 


J /tl;K* {892)^) 


B{B^^ 


V'(25) if*(892)0) 


B{B^ J/V' if* (892)0) 


B{B^ -> 


(/.'(25) if* (892)0) 


B{B+ 


^ip{2S) K+) 



7.2. FUTURE PROSPECTS 
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one standard deviation for all but one of the six branching-fraction ratios. This set of 
predictions, due to Neubert et al. [138|, assumed a dipole behaviour for the majority of 



the form factors, which were calculated from solutions to a relativistic harmonic oscillator 
potential model |73, 74]. The somewhat less- favoured set of predictions assumed single-pole 



form factors that were estimated from measurements of semileptonic D-meson (heavy-to- 



light) decays |81|. 



Relative branching-fraction measurements were combined with the world-average branch- 
ing fraction for the decay i?+ J/ip Q, (1.01 ±0.14) x 10^'^, to extract the following 
absolute branching fractions: 

BiB° J/V'i^*(892)0) = (1.78±0.14[stat] ±0.29[syst]) x lO^^ 

B{B+ 1^(25) K+) = (0.56 ± 0.08[stat] ± 0.10[syst]) x 10"^ 

B{B° ipilS) K*{892f) = (0.92 ± 0.20[stat] ± 0.16[syst]) x 10^^ 

where the first uncertainties are statistical and the second are systematic. The derived 
absolute branching fractions for the two S-meson decays to ip{2S) final states were shown 
to be consistent with both recent measurements and limits at e"*" e~ colliders and recent 
phenomenological predictions. 

The branching-fraction measurements reported here are consistent with the world- 
average values compiled by the Particle Data Group (refer to Table |67^ ) . The measured 
branching fractions that involve ip{2S) final states constitute the world's most precise mea- 
surements of these quantities. 

7.2 Future Prospects 

Using the same data sample as that described in the present study, measurements of sev- 
eral other similar exclusive i?-meson decay processes are feasible. These processes in- 
clude the analogous decays S° J/ipK^, B° ^{28) K^, B+ J/^jj K*{892)+ , and 
B^ ip{2S) i^*(892)"'', which, when combined with the results of the present work, could 
yield a measurement of the ratio of fragmentation fractions, fu/ fd- Decays of the B^ meson 
into cc final states, B^ .J/ip 0(1020) and B^ tp{2S) (j){1020), also have similar topologies 
and reconstruction techniques. Relative branching fractions of some of these modes were 
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reported for Run lA data |143, |9^. Other recent work [|128(| also used the Run lA data set 
to measure the longitudinal polarization fractions in the decays J / il) K* {S%2)^ and 

S° J/%1) (/)(1020), a study that can be extended to use the full Run 1 sample and to mea- 
sure Tl/T for the other vector- vector decays mentioned above. In particular, improved mea- 
surements of the quantity —{B^ J / ifj K* {S92)^) are needed for comparison with mea- 
surements from e'^ colliders pH, |144|] . Finally, the large samples of ip{2S) J/ipir^ ir^ 



and ij{2S) A*^ decays in Run 1 (refer to Chapter |5|) could be used to measure the 

relative branching fraction between the two charmonium modes, thereby reducing the total 
systematic uncertainty of the i?-meson results reported in this dissertation. 

The Fermilab Tevatron is currently being upgraded with a 'Main Injector' to oper- 
ate at a centre-of-mass energy of ^/s = 2.0 TeV with instantaneous luminosities of up to 
2 X 10'^^ cm~^s^^. Scheduled to begin in 1999, the next run of the Tevatron (Run 2) is 
expected to accumulate a data sample with time-integrated luminosity JCdt = 2 fb^^, a 
twenty-fold increase over that for the sample used in this study. Commensurate with the 
improvements to the Tevatron, the CDF detector is being rebuilt into the upgraded 'CDF 
ir detector [145|. The tracking systems will have redesigned and integrated silicon and 
drift-chamber detectors with improvements in the combined pattern recognition, momen- 
tum resolution, and pseudorapidity coverage. The muon systems will be extended to cover 
nearly all angles in azimuth and pseudorapidities in the range |?7| < 2.0. The data acquisi- 
tion system is expected to handle 132-ns bunch crossings, and the trigger will have no dead 
time with the added improvements of tracking information at Level 1, impact parameter 
discrimination at Level 2, and 300-IIz operation at Level 3. The CDF II and Tevatron 
upgrades are expected to augment significantly both the scope and reach of inquiry into the 
properties and interactions of the bottom quark well into the next millennium. 



Appendix A 

Level 2 Dimuon Triggers 



This appendix describes the Level 2 dimuon triggers that were used to construct the data 



samples studied for this analysis. Refer to Section 2.5.2 for a description of the Level 2 



system. For each Level 2 trigger, the associated Level 1 prerequisite dimuon trigger is 



given. A description of the Run lA and Run IB running periods is given in Section 3J- 

A.l Run lA Level 2 Dimuon Triggers 

1. TWO_CMU_ONE_CFT: Two CMU Level 2 muon clusters are necessary, with one of 
the clusters required to match a CFT track. The Level 1 prerequisite trigger is 
two_cmu_3pt3. This trigger is dynamically prescaled. 

2. cmx_one_cmu_CFT: One CMX Level 2 muon cluster and one CMU Level 2 muon 
cluster are necessary, with the CMU cluster required to match a CFT track. The Level 
1 prerequisite trigger is two_CMU_CMX_3pt3. This trigger is dynamically prescaled. 

3. cmu_OR_CMX_one_cmx_cft: Either two CMX Level 2 muon clusters or both a CMU 
cluster and a CMX cluster are necessary, with a CMX cluster required to match a 
CFT track. The Level 1 prerequisite trigger is two_cmu_CMX_3pt3. This trigger is 
dynamically prescaled. 

A. 2 Run IB Level 2 Dimuon Triggers 

1. TWO_CMU_TWO_CFT: Two CMU Level 2 muon clusters are necessary, with each cluster 
required to match a CFT track. The Level 1 prerequisite trigger is two_CMU_3pt3. 
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This trigger is not prescaled, and the two muon clusters are required to be in non- 
contiguous modules. If one cluster is in the +z region of the CDF detector and the 
other cluster is in the —z region, then the two clusters are required to have different 
(f values. 

2. cmx_cmu_two_CFT: One CMU Level 2 muon cluster and one CMX Level 2 muon 
cluster are necessary, with each cluster required to match a CFT track. The Level 1 
prerequisite trigger is TWO_cmu_CMX_3pt3. This trigger is dynamically prescaled. 

3. TWO_CMX_TWO_CFT: Two CMX Lcvcl 2 muon clusters arc necessary, with each cluster 
required to match a CFT track. The Level 1 prerequisite trigger is TWO_CMU_CMX_3pt3 
This trigger is dynamically prescaled, and, if both CMX muon clusters are in the same 
hemisphere (in z) of the CDF detector, then they are required to be separated by at 
least one wedge. 

4. two_CMU_ONE_CFT_6tow: Two adjacent CMU Level 2 muon stubs, which together 
form a single cluster that spans six or more calorimeter trigger towers, are necessary. 
The single cluster is required to match a CFT track, and the Level 1 prerequisite 
trigger is TWO_cmu_3pt3. This trigger is dynamically prescaled and is intended to 
offset losses due to the discontiguities imposed in the TWO_CMU_TWO_CFT trigger. 
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CDF Track Characterization 



This appendix introduces the mathematical formulation [122, |l4q ] , adapted for use in the 
present study, that describes helical trajectories of charged particles in the CDF detector. 
Five parameters are used to define a track's helix: the curvature, c; the azimuthal angle, ipo, 
at the point of closest approach to (x,y) = (0,0); the cotangent of the polar angle, cot 6; 
the impact parameter with respect to {x,y) = (0,0), do; and the local z coordinate with 
respect to {x,y) = (0,0), zq. In Section |B]|, the global coordinates of a particle's trajectory- 



are derived from the track helical parameters; in Section B.2, a method to derive the track 



helix using only point-of-origin and momentum information is outlined. 

B.l Global Coordinates from the Track Helix 

The curvature of a track is defined as 

c = ^, (B.l) 
2ro 

where tq is the radius of curvature and Q denotes the sign of the electric charge of the 
particle. It is assumed here that \Q\ = 1. 

In order to determine the global coordinates of a track in the CDF detector, the az- 
imuthal deflection of the particle's momentum, Aip, is defined as 

EE i^, (B.2) 

where s is the arc length of the helix, projected onto the transverse {x, y) plane. By 



substitution of Equation B.l into Equation B.2, the global azimuthal angle, 93, is defined in 
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radians as 



ip = (Aip + ifo) mod 27r = (2cs + ipo) mod 27r. 



(B.3) 



If the coordinates (xo,yo) define the transverse location of the centre of the track helix, 
then the global transverse coordinates of the helix are 



X = xo + roQ sin (p 
y = yo-roQ cosip. 



(B.4) 



The same form as that in Equation |B.4| may be used to calculate the coordinates of the 
point in the track trajectory that is nearest, in the transverse plane, to the geometric centre 
of the detector, {x',y'): 



y 



xo + roQ sin ipQ 
yo-roQ cos 930- 



(B.5) 



Since the magnitude of the impact parameter, \dQ\, is defined as |do| 
its sign is defined as 



dp 

\do\ 



^o + yo- ^0 



the coordinates {x',y') may also be written as 



and 



(B.6) 



x' = —do sin(^o 
y' = do cosv?o- 



(B.7) 



Solving for (xo,yo) in Equations B^ and B^ yields 



Xo = -{roQ + do) sin(^o 
yo = (roQ + do) cos 990- 



(B.8) 



Substitution of Equation B.8 into Equation B^ , and translation of the zq helix parameter 
into its analogous global coordinate, furnishes the following expressions for a track's global 



B.2. HELIX FROM MOMENTUM AND POINT OF ORIGIN 117 



coordinates: 

X = roQ sin (f - {roQ + do) sin 930 

y = -ro Q cos (f + {ro Q + do) cos ipo (B.9) 

z = zo + s cot 9. 

B.2 Helix from Momentum and Point of Origin 

In studies that involved the embedding of Monte Carlo generated particle tracks into actual 
data events, a transformation from momentum, electric charge, and spatial-origin coordi- 
nates to local helical parameter coordinates was necessary. The transformations used, as 
well as any assumptions made, are detailed in the following. 
The curvature parameter was calculated with the formula 

c = g-ISSaSxlO-^X^^. (B,io) 

where B is the magnetic field in units of Tesla and pT is the transverse momentum, = 
\p\ sinO, in units of GeV/c. 

The azimuth parameter was determined with the expression 

{arccos ( — ) if p„ > 

yP^J , (B.ll) 

27r - arccos iipy<^ 

where px and py are the x and y components of the momentum vector, respectively, and an 
approximation]^ that 939 = 9^ is assumed. 

The cotangent of the polar angle was simply computed using 

cote = ^, (B.12) 

where Pz is the z component of the momentum vector. 

The calculation of the impact parameter, do, required care to avoid a numerical insta- 
bility. The impact parameter may be defined as 

do = Q(e-^o), (B.13) 



^The stated approximation is justified in track embedding studies where only a subset of the Monte 
Carlo tracks generated for a given event is considered and there is therefore no need to preserve the precise 
kinematics of the parent particle. 
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where ^ = \J Xq + y^. In practice, Equation B.13 is difficult to calculate since |do| i- This 
was remedied by rewriting Equation p.l3| in the form jl46[ 



Q{i- rp) (g + ro)/ro 



(B.14) 



The coordinates {x",y") were defined to represent the displacement of the track's point of 



origin. A substitution of an analogous form of Equation B.4 was made to convert Equa- 



tion B.14 to the form 
do = 

Further substitution using Equation |B.1| yielded the practical form 
dn = 



Q [{x'T + jy'T] + 2ro [y" cosy.p - x" simpo] 
ro + J {x" -roQ sin ipof + {y" + roQ cos v^o)^ 



(B.15) 



2c [(x")^ + (y")^] + 2 [y" COSV30 - x" sint^o] 

-roQ sin ipof + {y" + roQ cos ipo)'^ 



(B.16) 



l + 2cQ 

Following the same justification as that used in the ipo = ip assumption in Equation B.ll| , 

— s cot 9, was approximated as 



the value of zq, which according to Equation B.£ is zq 



(B.17) 



Appendix C 

Pattern Recognition Efficiencies 



This appendix summarizes a study |135| ] of the single- and double-track pion pattern recog- 
nition efficiencies in the CDF central tracking chamber (CTC) over the course of the Run 
1 data-taking period. 

C.l Data Sample 

The sample used in this efficiency study consisted of ~36 000 raw dimuon events from 40 
data-taking runs, 10 from Run lA and 30 from Run IB. The runs represented in this sample 



are detailed in Table C.l. 



The four loose run-range bins in Table C.l are delimited by the changeover period be- 



tween Runs lA and IB, a major reorientation of the Tevatron beam optics, and a prolonged 



shutdown of the collider, respectively. The left side of Figure |C.l| illustrates the instanta- 
neous luminosity profile of ~530 000 — > J/ip candidates, including background, and 
clearly features the natural divisions between the four run bins. Portrayed on the right side 



of Figure C.l is the profile of the event sample used in the present efficiency study. 



C.2 Track Embedding Procedure 

The crux of this study was the embedding of hits representing two Monte-Carlo-generated 
pions into raw data events. mesons were generated in a manner similar to that described 



in Section |4.1| . The CLEO QQ programme [|127| was used, as in Section [4.2| , to decay 
the mesons in the chain B^ ip[2S) K'^ , V'(25') — > J/ipiT'^TT~, and J/ip ^ fi^ ^~ . 
The dipion invariant mass distribution was generated according to the customized matrix 
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Experiment 


Run lA 




Run IB 




Run 




Low 


Medium 


High 


Raiigc- 


iO 100-17 835 


55 212-G0 900 


60 901-67 390 


G7 391-71 023 


R 


40 512 


55 571 


61024 


68040 


u 


41085 


56 002 


61548 


68 206 


n 


42 030 


57006 


63 171 


68 517 




43 016 


57513 


63 502 


68 758 


N 


44366 


58093 


64041 


69036 


u 


45 047 


58 564 


65 004 


69 550 


m 


46 031 


59042 


65 560 


69 809 


b 


46 519 


59 517 


66 020 


70019 


e 


47008 


60004 


66 539 


70 559 


r 


47 552 


60 597 


66 615 


71013 


J£dt [pb-^] 


20.90 


9.46 


47.08 


35.41 



Table C.l: A list of runs representing four logical run-range divisions of the Run 1 data- 
taking period. Each of these four divisions contributed ^^9 000 events to the sample. The 
bottom row of the table lists the time-integrated luminosities represented by the various 
run-range bins. 
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Figure C.l: Instantaneous luminosity vs. run number profile for (left) ^^530 000 5+ — 
J/'0K+ candidates and that for (right) ~36 000 raw dimuon events used in the present 
efficiency study. 
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element described in Section 4.2.1. For events where the vr^ and particles each passed 
a pt > 250 MeV/c transverse momentum requirement, the 3- momentum and charge of the 
two pions in the decay chain were extracted from the CLEO QQ output and stored for use 
in the embedding stage. 

The curvature (c), azimuth {<po), and polar angle (cot 0) helical parameters of the tracks 
to be embedded were derived from the 3-momentum and charge information furnished in 
the Monte Carlo track generation stage. The impact (do) and longitudinal displacement 
(zq) helical parameters were calculated with additional information on the primary vertex 
location of the particular event into which the tracks were to be embedded. Appendix |^ 
describes the method by which the track helices were determined. 

Once all five helical parameters for each of the two pion tracks were known, a set 



of routines [147| was used to embed hits into the CTC and VTX data structures. The 
values of the embedded parameters were recorded for later comparisons, and the full default 



production executable package (described in Section 3^) was used to reconstruct the tracks 
in the events. 

In order to determine the wire hit efficiencies to use in the embedding software, the 
mean numbers of hits used by the reconstructed track fits in each CTC superlayer were 
examined, as functions of instantaneous luminosity, for both embedded and data tracks. 
Figures C.2 and |C.3| show the mean numbers of hits per superlayer used in fits of tracks 



embedded with 100% wire hit efficiencies, as functions of instantaneous luminosity, for the 
axial and stereo superlayers, respectively. 

The mean numbers of hits per superlayer used in the embedded track fits were tuned to 
match those observed in data by iteratively adjusting the superlayer-by-superlayer wire hit 
efficiencies in the embedding routines. Table p.2| summarizes the reduced wire efficiencies 
calculated using this technique. 



It is evident from Figures |C.2| and C.3 that there was negligible sensitivity to the charge 



of the embedded track. Consequently, the wire efficiencies featured in Table C.2 were applied 
to both positively and negatively charged tracks. Within the axial and stereo superlayer 
categories in Table p.2| , the wire efficiencies were observed to diminish monotonically from 
the outer to the inner superlayers. Due to the geometry of the CTC and the pT distribution 
of charged tracks |12S], increased charged-track occupancies near the inner-superlayer radii 
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Figure C.2: The mean number of hits per axial superlayer used in fits of tt~^ (left) and vr^ 
(right) tracks embedded with uniform 100% wire hit efficiencies, as a function of instanta- 
neous luminosity. 
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Figure C.3: The mean number of hits per stereo superlayer used in fits of vr^ (left) and tt~ 
(right) tracks embedded with uniform 100% wire hit efficiencies, as a function of instanta- 
neous luminosity. 
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3 


30-35 
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18 - 29 


0.80 




1 


12 - 17 


0.82 







0-11 


0.74 



Table C.2: The wire efficiencies used to embed hits in the 9 CTC superlayers. These were 
calculated by tuning the mean number of fitted hits per superlayer to match in embedded 
and data samples. 

accounted for the bulk of the observed relative trends in the superlayer-dependent wire 
efficiencies. 

C.3 Definition of a "Found" Track or Track Pair 



The technique used to identify a "found" track attempted to account for interparameter 
correlations in the reconstructed helices. Specifically, for a given Monte Carlo track with 
embedded helical parameter vector a, where Oi £ {cot6,c,ZQ,dQ,LpQ), each reconstructed 
track in the event, S, that shared the same sign of curvature with ct, was used to calculate 



a matching x quantity via the expression [30 



X 



{a — S)'^ V (a 



a 



(C.l) 



where the matrix V was the 5x5 covariance matrix of the reconstructed track helix. For 
each embedded track charge, the reconstructed track in the event that formed the lowest 
matching value with the embedded track, Xmin' retained for consideration as a 
"found" track. 

There was some speculation about the relative effects of the stereo and axial contribu- 
tions to inefficiencies in the CTC. Tracks were initially identified in the CTC by seeking 
segments in the axial superlayers and forming 2-dimensional track objects in the r-ip view. 
Stereo reconstruction then took place to form 3-dimensional tracks, using the information 
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from the stereo superlayers and from the VTX, which provided the initial information on 
the zq hehx parameter of the tracks. It was therefore predicted that, under increased in- 
stantaneous luminosity conditions that precipitated higher primary vertex muhiphcities in 
the VTX, the rehability of the zq determination, and hence the stereo reconstruction effi- 
ciency, would suffer. Occupancy-related luminosity-dependent efficiency losses in the inner 
CTC superlayers were also expected to afflict the stereo more than the axial reconstruction 
because of the lesser number of superlayers and wires per superlayer in the stereo view (refer 



to Section |2.3.3| for a description of the CTC). 

The foregoing considerations about axial and stereo contributions to tracking perfor- 
mance motivated the construction of an axial-only matching quantity. Specifically, the 



vector a in Equation C.l was modified such that a, G (c, do, ^o) and the matrix V became 
the corresponding 3x3 track parameter covariance matrix. 

In the case of 3-dimensional track matching, a reconstructed track was designated as 
"found" if its matching value satisfied the criterion 

Xmin < 500. 

The analogous criterion in the case of 2-dimensional axial-only matching was 

xLn < 300. 

The efficacy of these cuts in discriminating between the embedded track and other tracks in 
the event is exemplified in Figure C.4, which plots the minimum matching x^, Xmirn versus 



the next-to-minimum matching x^, Xnext) the 3-dimensional case. Figure p.4| furnishes 
two useful observations: first, the matching selection criteria, indicated by the horizontal 
lines, were relatively efficient at determining the distinction between the found track and 
any other track in the event; and second, the majority of events where an embedded track 
was not found at all had Xmm ~ Xnext ' expected. 

The pt dependence of tracks that passed and failed the "found" -track identification 
criteria is pictured in Figure |C.5| , this time for 2-dimensional matching. These plots indicate 
that any inefficiencies in the Xmin matching criteria were largely confined to embedded tracks 
with extremely low px- These figures also suggest that the majority of embedded tracks 
that failed to be "found" had low p^. 
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Figure C.4: The minimum 3-dimensional matching x^i Xniim plotted against the 3- 
dimensional next-to- minimum matching x^, Xnexti 3 000 embedded vr"*" mesons (top) and 
7r~ mesons (bottom). The horizontal line designates the selection criterion, Xmin < ^^O- 




Figure C.5: The minimum 2-dimensional matching x^, Xmin' plotted against the embedded 
track transverse momentum, for 3 000 embedded tt'^ mesons (top) and 7r~ mesons 
(bottom). The horizontal line designates the selection criterion, Xmm < "^^O. 
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The embedding of both ir^ and vr" mesons into the same event enabled the determination 
of a "found" track pair. A track pair was considered to be "found" if, in the case of 
3-dimensional track matching, 

max(xmmk^].Xmink"]) < 500 

and, in the case of 2-dimensional track matching, 

max(xmink+],Xmm[7r"]) < 300. 

C.4 Pattern Recognition Efficiency 
C.4.1 Definition 

The single- and double-track CTC pattern recognition efficiencies may be expressed as a 



function of one of several observables. Equation C.2 defines the bin- by-bin pattern recog- 



nition efficiency for a given bin, A^, in an observable ^: 

(C.2) 



/ iVfound(e)rf? 

em- 



/ [A^found(e)+A^foWd(0] 
A? 



CTC \ CTCJ 



The A'^foundCO symbol represents the number of tracks (or track pairs) at a given value of 



^ that are identified as "found" according to the criteria discussed in Section C.3. The 
- ^found (^) symbol represents the corresponding "found"-track (or track pair) failures, and 
the quantity Tq^c CTC exit radius, or the radius at which tracks intersect one of the 

two endplate planes of the CTC. Note that the px or rfj^Q requirements in Equation C^ 



are not applied when ^ = p-^ or = r^^Q, respectively. Also, when paired track efficiencies 
are determined as a function of pT, ^ = min (pT[7r'''],PT[7r~]) in Equation |C.2| . 

An aggregate pattern recognition efficiency may be calculated in terms of ^ = for 
Pt > Px", as 



elpT > Pt 



E iVfound(ApT) 



E [A^found (Apt) + Apt) 



(C.3) 



^CXlt /^CXlt ^ 

CTC ^ V CTC^ 



with binomial statistics used to compute the statistical uncertainty on e. 
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Figure C.6: The efficiency dependence on the CTC exit radius for tt~^ (top) and vr (bottom) 
tracks. 



C.4.2 Efficiency Dependence on Kinematic Observables 



Employing the definition described in Equation C.2, the pattern recognition efficiency was 
studied for tracks in several different kinematic observables. Unless noted otherwise, p™™ = 
0.4 GeV/c in Equation 

Figure p.6| depicts the efficiency dependence on r^'rc both vr"*" and vr" tracks with 
no (?'cTc)™™ criterion. The shape of the turn-on motivated the use of a (j^cTc)™™ criterion 
less than the often-used 132-cm one, which is conservative and corresponds to the outer edge 



of superlayer 8. Figure C.t suggests that a conservative criterion of (?'cTc) ~ 

was appropriate for the provision of reliable tracking efficiency information for muon tracks, 

without any significant compromise in the CMX acceptance. The remaining figures in this 



appendix (except the left side of Figure C.13 ) have been prepared using (?'cTC 



110 cm, 



which corresponds to a radius just outside the outer edge of superlayer 6. 

Figure |C.7 presents an example of the px-dependent pattern recognition efficiencies for 
single vr"*" and vr" tracks and vr^ vr^ track pairs. The data points were calculated using 
Equation |C]| with pJf'"" = and (rgj,*c) = ^0 

cm. The indicated aggregate efficiencies 



128 APPENDIX C. PATTERN RECOGNITION EFFICIENCIES 




0.0 



0.20 



0.45 0.70 0.95 

Pt [GeV/c] 



1 .20 " 



1.2 

o 0.8 
c 

h OA 



£(Pt > 0.4 GeV/c) = 0.947 ± 0.003 



0.0 I— L 
0.20 



J L_l l_l l_J , I , l_J L 

0.45 0.70 0.95 

Pt [GeV/c] 



1.20 



1.2 



1.0 



O.f 



0.6 



0.4 



0.2 



0.0 
0. 



-f-f 



£(P,""" > 0.4 GeV/c) = 0.915 ± 0.005 



20 0.45 0.70 0.95 1.20 

minimum P, [GeV/c] 



Figure C.7: An example of the efficiency dependence on transverse momentum for single 
tracks (left) and track pairs (right), in this case for the Run IB "Low" run-range sample. 
The indicated measured aggregate efficiencies represent this run range alone; the uncertain- 
ties are statistical. 



were calculated using Equation p.3| , but with p™™ = 0.4 GeV/c; the uncertainties are 
statistical only. Note that the product of the two single-track efficiencies, 0.894 it 0.004, 
differs from the two-track efficiency by ~3.3 standard deviations, indicating that the pattern 
recognition efficiencies for two tracks in a single event are correlated. This observation is 



discussed further in Section C.4.3 



C.4.3 Efficiency Dependence on Occupancy-Related Observables 

The effects of CTC occupancy on the pattern recognition efficiency were traditionally con- 



templated in terms of instantaneous luminosity dependence. Figure |C.8| shows plots of the 
single- and two-track efficiencies as functions of instantaneous luminosity. 

For a given event, the instantaneous luminosity was only indirectly related to the true 
track occupancy in the CTC fiducial volume. A more direct measure of occupancy is 
provided by the multiplicity of high-quality primary vertices in an event. In Figure p.9| , 
the single- and two-track efficiencies as functions of the number of Class-12 vertices^ are 



^The highest quality z vertex identifiable using the VTX was deemed a 'Class-12' vertex. Class-12 vertices 
were defined as such when the number of associated hits in the VTX exceeded 180. 
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illustrated. Note that the efficiency dependence on primary vertex multiplicity is greater 
than that on instantaneous luminosity by more than a factor of two. 

It is useful to examine the efficiency dependence on variables that are considered to 
be even more correlated than primary vertex multiplicity with the true occupancy of the 
CTC. An example of one of these is the total number of fitted CTC tracks in each event. 
In spite of the fact that the CTC track fits themselves used primary z-vertex information, 
the CTC track multiplicity was expected to be more correlated with the CTC occupancy 
because the track multiplicity was measured in the CTC itself and, for example, would 
be largely independent of those primary vertices that represented tracks at high absolute 
pseudorapidities, outside the fiducial acceptance of the CTC. Figure p.lO| portrays the 
single- and double-track efficiencies as functions of the number of CTC tracks. Although 
the distributions are limited by poor statistics at low and high multiplicities, the efficiencies 
exhibit a significant dependence on CTC track multiplicity. 

The correlation between the pattern recognition efficiencies of two tracks in a given 



event was discussed briefiy in Section C.4.2. Figure C.ll compares the efficiency of a 



track pair, ^(vr^ ^ )> with the efficiency product of two individual tracks, e{TT^) e{iT^), as 
a function of Class-12 primary vertex multiplicity. For the purposes of this figure, a more 
restrictive Pt{tt^) > 0.5 GeV/c criterion was imposed to reduce the tracking differences 
between positively and negatively charged tracks; the exit radius cut, however, was kept at 
''cTC > 110 cm. The central values of the e(7r+ vr^) measurements were consistently greater 
than those for the e('7r+) ^(vr^) measurements, but the two quantities were only statistically 
inconsistent (by at least one standard deviation) in events with only one or two Class-12 
primary vertices. In spite of the fact that the events with only one or two Class-12 vertices 
constituted a major fraction of the total number of events in the sample, the degree of 
correlation between two tracks in a single event, while interesting, was unobservable when 
systematic uncertainties were considered. 

C.4.4 Time-Dependent Effects 



Section |5.2| mentioned several time-dependent factors that may have been detrimental to the 
CTC tracking efficiencies in Run 1. Time-dependent effects can be organized into two broad 
categories: effects due to occupancy and those due to aging and other sources. Although 
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Figure C.8: The efficiency dependence on instantaneous luminosity for single tracks (left) 
and track pairs (right). 
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Figure C.9: The efficiency dependence on multiplicity of Class-12 (high-quality) vertices for 
single tracks (left) and track pairs (right). 
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Figure C.IO: The efficiency dependence on multiplicity of CTC tracks for single tracks (left) 
and track pairs (right). 



the time dependence of the latter category is typically studied in terms of run number, 
time-integrated luminosity is a more meaningful quantity for this purpose. 



It is apparent from the results presented in Section |C.4.3| that the CTC pattern recogni- 
tion efficiencies depended strongly on the occupancy conditions inside the CTC. The results 
from Section C.4.3| also suggest that variables similar to primary- vertex and CTC-track mul- 
tiplicity be used to constrain occupancy levels in studies of run-dependent aging effects in 
the CTC. 



As outlined in Section C.2, hits were embedded into the CTC with wire efficiencies 



that were tuned from the data. The embedding techniques used in this study therefore 
did account for some time-dependent performance effects (in addition to occupancy), to the 
extent that the wire efficiencies used in the embedding were derived from data quantities (the 
numbers of used hits per superlayer) that were expected to be sensitive to time-dependent 
performance effects in the CTC. It is important to note, however, that the present study 
cannot be expected to resolve a posteriori any time-dependencies that affected the embedded 
hits directly, since the wire efficiencies were calculated using run-averaged effects and were 
administered uniformly for all the events in the sample. It should be emphasized that, 
with the exception of the embedded hits, this technique can potentially resolve all other 
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Figure C.ll: A comparison of the two-track efficiency and the product of single-track ef- 
ficiencies, as a function of primary vertex multiphcity. A tighter pt{'^^) > 0.5 GeV/c 
criterion was imposed to reduce the tracking differences between positively and negatively 
charged tracks, but the exit radius cut was kept at r^^^Q > 110 cm. The indicated uncer- 
tainties are statistical only. 
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Figure C.12: The dipion efficiency dependence on run number for all Class- 12 primary 
vertex multiplicities (left) and for one and only one Class-12 primary vertex (right). 

time-dependent effects in each of the ~36 000 events. 

Figure |C.12 illustrates the dependence of the pattern recognition efficiency on run num- 



ber. For the case where all primary vertex multiplicities are permitted, a negative slope 
in the distribution is discernible. When the efficiency is examined in events with only one 
Class-12 vertex, the slope is observed to vanish within the available statistics. The trends 



in Figure C.12 are consistent with the expectations of this embedding technique and the 



conclusion that the most pernicious effect on the pattern recognition is one of occupancy. 
C.4.5 Stereo Efficiency 



Pursuant to the discussion in Section p.3| about relative stereo and axial contributions 
to the tracking inefficiency, a stereo pattern recognition efficiency was defined, also using 



Equation |C.2| , but with the denominator populated by the subsample of events that satisfied 
the 2-dimensional axial-only matching criteria. In this expression of the stereo efficiency, 
-^found represented the number of events that, in addition to meeting the axial-only criterion, 



satisfied the 3-dimensional matching requirement. Figure p.l3 shows plots of the stereo 



pattern recognition efficiency as a function of two kinematic observables: the CTC exit 



radius, rfj^Q, and the transverse momentum, px- The stereo efficiency in Figure C.13 was 
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Figure C.13: The stereo pattern recognition efficiency for vr^ (top) and vr^ (bottom) tracks 
as a function of two kinematic observables: (left) the CTC exit radius, Tq^c, and (right) 
the transverse momentum, pT- 



observed to remain constant near 100% in r^^^Q, but to diminish appreciably in for < 
0.4 GeV/c. The dependence of the stereo pattern recognition efficiency on occupancy-related 
observables is illustrated in Figure p.l4 , which features plots of the stereo efficiency as a 
function of the primary-vertex and CTC-track multiplicities. The stereo pattern recognition 
efficiency appears to drop by ~10% with increasing primary-vertex multiplicity and by 
~5% with increasing CTC-track multiplicity. No significant differences in stereo efficiencies 
between vr^ and vr^ tracks are apparent in Figures lOl^ and |Cl4l . 



C.5 Results 



Single- and double-track pattern recognition efficiencies were calculated for each of the four 
run-number bins using Equation |C.3| with p™™ = 0.4 GeV/c. Tables and CA list 



the results of these calculations for rfj^c > HO cm with 3-dimensional and 2-dimensional 
matching criteria, respectively. Similarly, Tables |C.6| and |C.7| list the results of the calcu- 
lations for r^r^Q > 132 cm with 3-dimensional and 2-dimensional matching criteria, respec- 



tively. Tables C^ and C.8 summarize the aggregate Run 1 single- and double-track pattern 
recognition efficiencies for r^^pQ > HO cm and Tq^c > ^"^^ respectively. The totals 
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Figure C.14: The stereo pattern recognition efficiency for tt~^ (top) and vr" (bottom) tracks 
as a function of two occupancy-related observables: (left) the Class-12 VTVZ primary vertex 
multiplicity and (right) the CTC track multiplicity. 



were computed by calculating mean values weighted in time-integrated luminosity (refer to 
Table |C.1| ) and statistical uncertainty. For each aggregate efficiency, the mean statistical 
uncertainty was added in quadrature with the systematic uncertainty, which was taken to 
be equal to the maximum difference between the tracking efficiencies calculated with full 
(100%) wire hit efficiencies and those calculated with reduced wire hit efficiencies. 



C.6 Conclusions 

This study of Run lA and IB single- and double-track pattern recognition efficiencies yielded 
several qualitative and quantitative conclusions: 

1. To model accurately the CTC performance in Run 1, it was necessary to employ 
reduced wire hit efficiencies in the track embedding procedure. 

2. Throughout this study, 7r+ and vr^ tracks were treated separately. Within the statis- 
tical uncertainties alone, no significant differences between efficiencies for positively 
and negatively charged tracks were observed. 



3. A CTC track exit radius criterion of Tq^^q > 110 cm ensured that the track was in 
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3-D Efficiency 


e(vr+) 


e(vr-) 


e(7r+7r ) 


Run lA 


.952 ± .003 ± .017 
.944 ± .003 ± .016 
.921 ± .004 ± .023 
.899 ± .004 ± .030 


.956 ± .003 ± .018 
.947 ± .003 ± .014 
.923 ± .003 ± .023 
.900 ± .004 ± .026 


.924 ± .004 ± .029 
.915 ± .005 ± .023 
.869 ± .006 ± .036 
.827 ± .006 ± .043 


Run IB 


Low 

Medium 

High 



Table C.3: The single- and double-track 3-dimensional pattern recognition efficiencies mea- 
sured in the four run-number bins for tracks with px > 0.4 GeV/c and rg^c > 
The first uncertainties shown are statistical and the second are systematic. 



2-D Efficiency 




siTT-) 


e(7r"'' vr ) 


Run lA 


.957 ± .003 ± .017 
.949 ± .003 ± .016 
.933 ± .003 ± .021 
.918 ± .004 ± .025 


.960 ± .003 ± .017 
.952 ± .003 ± .013 
.931 ± .003 ± .020 
.916 ± .004 ± .019 


.931 ± .004 ± .030 
.922 ± .005 ± .023 
.883 ± .005 ± .035 
.858 ± .006 ± .035 


Run IB 


Low 

Medium 

High 



Table C.4: The single- and double-track 2-dimensional pattern recognition efficiencies mea- 
sured in the four run-number bins for tracks with jot > 0.4 GcV/c and rg^^ > 110 cm. 
The first uncertainties shown are statistical and the second are systematic. 



Aggregate Efficiency 


Matching 




e(vr-) 


e(7r+ vr ) 


Run IB 


3-D 


.917 ± .030 


.920 ± .026 


.860 ± .043 


only 


2-D 


.931 ± .025 


.930 ± .020 


.880 ± .035 


Runs 


3-D 


.927 ± .030 


.928 ± .026 


.881 ± .043 


lA + IB 


2-D 


.937 ± .025 


.936 ± .020 


.895 ± .035 



Table C.5: The total Run IB and lA -|- IB single- and double-track pattern recognition 
efficiencies, each calculated by combining the relevant run-number bins into a single mean, 
weighted by the appropriate statistical uncertainties and time-integrated luminosities. Val- 
ues are given for both the 3-dimensional and 2-dimensional matching techniques for tracks 
with pt > 0.4 GeV/c and r^c > HO cm, and the uncertainties are quadratic sums of the 
statistical and systematic components. 



3-D Efficiency 


e{n+) 




£(7r"^ TT ) 


Run lA 


.959 ± .003 ± .015 
.951 ± .003 ± .013 
.929 ± .004 ± .021 
.906 ± .004 ± .029 


.961 ± .003 ± .016 
.952 ± .003 ± .014 
.931 ± .003 ± .021 
.907 ± .004 ± .023 


.937 ± .004 ± .028 

.929 ± .005 ± .024 
.887 ± .006 ± .032 
.849 ± .007 ± .040 


Run IB 


Low 

Medium 

High 



Table C.6: The single- and double-track 3-dimensional pattern recognition efficiencies mea- 
sured in the four run-number bins for tracks with pT > 0.4 GeV/c and r^^Q > 132 cm. 
The first uncertainties shown are statistical and the second are systematic. 
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2-D Efficiency 




e(vr-) 


e(7r"'" TT ) 


Run lA 


.963 ± .003 ± .015 
.955 ± .003 ± .014 
.942 ± .003 ± .017 
.926 ± .004 ± .022 


.964 ± .003 ± .016 
.957 ± .003 ± .012 
.938 ± .003 ± .019 
.923 ± .004 ± .016 


.944 ± .004 ± .027 
.937 ± .004 ± .021 
.900 ± .005 ± .032 
.881 ± .006 ± .031 


Run IB 


Low 

Medium 

High 



Tabic C.7: The single- and double-track 2-dimensional pattern recognition efficiencies mea- 
sured in the four run- number bins for tracks with > 0.4 GeV/c and rg^c -'-^^ 
The first uncertainties shown are statistical and the second are systematic. 



Aggregate Efficiency 


Matching 


e(vr+) 


s(vr-) 


e(7r+ TT ) 


Run IB 


3-D 


.925 ± .029 


.927 ± .023 


.882 ± .040 


only 


2-D 


.939 ± .022 


.936 ± .019 


.901 ± .032 


Runs 


3-D 


.934 ± .029 


.935 ± .023 


.901 ± .040 


lA + IB 


2-D 


.944 ± .022 


.942 ± .019 


.913 ± .032 



Table C.8: The total Run IB and lA -|- IB single- and double-track pattern recognition 
efficiencies, each calculated by combining the relevant run-number bins into a single mean, 
weighted by the appropriate statistical uncertainties and time-integrated luminosities. Val- 
ues are given for both the 3-dimensional and 2-dimensional matching techniques for tracks 
with > 0.4 GeV/c and > 132 cm, and the uncertainties are quadratic sums of the 
statistical and systematic components. 
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a region of well-understood pattern recognition efficiency without undue compromise 



to the CTC and CMX fiducial acceptance (refer to Section C.4.2). 



4. The deleterious effects of CTC occupancy on the pattern recognition constituted the 
principal source of tracking inefficiency in Run 1. 

5. The pattern recognition efficiencies were studied for single and double tracks as a 
function of three occupancy-related observables, which were, in order of increasing 
correlation with the true CTC occupancy, the instantaneous luminosity, the multi- 
plicity of high-quality primary vertices, and the multiplicity of reconstructed tracks 
in the CTC (refer to Section C.4.3| ). It was recommended that investigations of non- 



occupancy time-dependent effects be performed with constraints using variables that 
were as closely correlated with the CTC occupancy as possible. 

6. The pattern recognition efficiencies of two tracks within a single event were correlated, 
a conclusion that followed from the observation that these efficiencies were driven by 
the event-by-event track- multiplicity environment conditions inside the CTC. In the 
context of primary-vertex multiplicity, the intertrack correlations were only statisti- 
cally significant at low primary-vertex multiplicities. Any double-track correlations 
became insignificant when systematic uncertainties were included in the efficiency 
measurements. 

7. CTC tracking performance in the stereo view alone appeared to contribute little to 
the observed pattern recognition inefficiencies. The stereo efficiencies possessed a 
weak dependence on occupancy and a strong dependence on pT in the region < 
0.4 GeV/c, which lay outside the practical kinematic range of this investigation (refer 
to Section pAS] ). 



8. Single- and double-track efficiencies were calculated in 2 dimensions (axial-only) and 
3 dimensions (axial and stereo) in four run-range bins spanning the Run lA and IB 
data-taking periods. Measurements were presented for tracks with > 0.4 GeV/c 
and either r^^Q > 110 cm or rg^^ > 132 cm. Aggregate Run 1 single- and double- 
track 3-dimensional and 2-dimensional pattern recognition efficiency measurements 
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were also computed for both the rfj^c ^ HO-cm and Tfj^c > 132-cm selection criteria. 



Refer to Tables U^, C.4, C.5, iJM, C.7, and C.8 in Section KlH. 



9. The total 3-dimensional Run lA + IB single-track efficiency for > 0.4 GeV/c and 
rgj,*c > 110 cm was measured to be 0.928 ± 0.020. 

10. The total 3-dimensional Run lA + IB double-track efficiency for > 0.4 GeV/c and 



rfjYQ > 110 cm was measured to be 0.881 it 0.043. 
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